SMILES-based QSAR model for arylpiperazines as high-affinity 5-HT1A receptor ligands using CORAL

In silico computational drug discovery: a Monte Carlo approach for developing a novel JAK3 inhibitors

Inhibition of Janus kinase 3 (JAK3), a member of the JAK family of tyrosine kinases, remains an essential area of research for developing treatments for autoimmune diseases, particularly cancer and rheumatoid arthritis. The recent discovery of a new JAK3 protein, PDB ID: 4Z16, offers exciting possibilities for developing inhibitors capable of forming a covalent bond with the Cys909 residue, thereby contributing to JAK3 inhibition. A powerful prediction model was constructed and validated using Monte Carlo methods, employing various internal and external techniques. This approach resulted in the prediction of eleven new molecules, which were subsequently filtered to identify six compounds exhibiting potent pIC50 values. These candidates were then subjected to ADMET analysis, molecular docking (including reversible-reversible docking with tofacitinib, an FDA-approved drug, and reversible-irreversible docking for the newly designed compounds), molecular dynamics (MD) analysis for 300 ns, and calculation of free binding energy. The results suggested that these compounds hold promise as JAK3 inhibitors. In summary, the new compounds have exhibited favorable outcomes compared to other compounds across various modeling approaches. The collective findings from these investigations provide valuable insights into the potential therapeutic applications of covalent JAK3 inhibitors, offering a promising direction for the development of novel treatments for autoimmune disorders.Communicated by Ramaswamy H. Sarma.

IFPTML Mapping of Drug Graphs with Protein and Chromosome Structural Networks vs. Pre-Clinical Assay Information for Discovery of Antimalarial Compounds

The parasite species of genus Plasmodium causes Malaria, which remains a major global health problem due to parasite resistance to available Antimalarial drugs and increasing treatment costs. Consequently, computational prediction of new Antimalarial compounds with novel targets in the proteome of Plasmodium sp. is a very important goal for the pharmaceutical industry. We can expect that the success of the pre-clinical assay depends on the conditions of assay per se, the chemical structure of the drug, the structure of the target protein to be targeted, as well as on factors governing the expression of this protein in the proteome such as genes (Deoxyribonucleic acid, DNA) sequence and/or chromosomes structure. However, there are no reports of computational models that consider all these factors simultaneously. Some of the difficulties for this kind of analysis are the dispersion of data in different datasets, the high heterogeneity of data, etc. In this work, we analyzed three databases ChEMBL (Chemical database of the European Molecular Biology Laboratory), UniProt (Universal Protein Resource), and NCBI-GDV (National Center for Biotechnology Information—Genome Data Viewer) to achieve this goal. The ChEMBL dataset contains outcomes for 17,758 unique assays of potential Antimalarial compounds including numeric descriptors (variables) for the structure of compounds as well as a huge amount of information about the conditions of assays. The NCBI-GDV and UniProt datasets include the sequence of genes, proteins, and their functions. In addition, we also created two partitions (c_assayj = c_aj and c_dataj = cd_j) of categorical variables from theChEMBL dataset. These partitions contain variables that encode information about experimental conditions of preclinical assays (c_aj) or about the nature and quality of data (c_dj). These categorical variables include information about 22 parameters of biological activity (c_a0), 28 target proteins (c_a1), and 9 organisms of assay (c_a2), etc. We also created another partition of (c_protj = c_pj) including categorical variables with biological information about the target proteins, genes, and chromosomes. These variables cover32 genes (c_p0), 10 chromosomes (c_p1), gene orientation (c_p2), and 31 protein functions (c_p3). We used a Perturbation-Theory Machine Learning Information Fusion (IFPTML) algorithm to map all this information (from three databases) into and train a predictive model. Shannon’s entropy measure Sh_k (numerical variables) was used to quantify the information about the structure of drugs, protein sequences, gene sequences, and chromosomes in the same information scale. Perturbation Theory Operators (PTOs) with the form of Moving Average (MA) operators have been used to quantify perturbations (deviations) in the structural variables with respect to their expected values for different subsets (partitions) of categorical variables. We obtained three IFPTML models using General Discriminant Analysis (GDA), Classification Tree with Univariate Splits (CTUS), and Classification Tree with Linear Combinations (CTLC). The IFPTML-CTLC presented the better performance with Sensitivity Sn(%) = 83.6/85.1, and Specificity Sp(%) = 89.8/89.7 for training/validation sets, respectively. This model could become a useful tool for the optimization of preclinical assays of new Antimalarial compounds vs. different proteins in the proteome of Plasmodium.

Curated Database and Preliminary AutoML QSAR Model for 5-HT1A Receptor

Introduction of a new drug to the market is a challenging and resource-consuming process. Predictive models developed with the use of artificial intelligence could be the solution to the growing need for an efficient tool which brings practical and knowledge benefits, but requires a large amount of high-quality data. The aim of our project was to develop quantitative structure–activity relationship (QSAR) model predicting serotonergic activity toward the 5-HT1A receptor on the basis of a created database. The dataset was obtained using ZINC and ChEMBL databases. It contained 9440 unique compounds, yielding the largest available database of 5-HT1A ligands with specified pKi value to date. Furthermore, the predictive model was developed using automated machine learning (AutoML) methods. According to the 10-fold cross-validation (10-CV) testing procedure, the root-mean-squared error (RMSE) was 0.5437, and the coefficient of determination (R²) was 0.74. Moreover, the Shapley Additive Explanations method (SHAP) was applied to assess a more in-depth understanding of the influence of variables on the model’s predictions. According to to the problem definition, the developed model can efficiently predict the affinity value for new molecules toward the 5-HT1A receptor on the basis of their structure encoded in the form of molecular descriptors. Usage of this model in screening processes can significantly improve the process of discovery of new drugs in the field of mental diseases and anticancer therapy.

Exploring biological efficacy of novel benzothiazole linked 2,5-disubstituted-1,3,4-oxadiazole hybrids as efficient α-amylase inhibitors: Synthesis, characterization, inhibition, molecular docking, molecular dynamics and Monte Carlo based QSAR studies

In an effort to explore a class of novel antidiabetic agents, we have made an effort to synergize the α-amylase inhibitory potential of 1,3-benzothiazole and 1,3,4-oxadiazole scaffolds by combining the two into a single structure via an ether linkage. The structure of synthesized benzothiazole clubbed oxadiazole derivatives are established by different spectral techniques. The synthesized hybrids are evaluated for their in vitro inhibitory potential against α-amylase. Compound 8f is found to be the most potent with a significant inhibition (87.5 ± 0.74% at 50 μg/mL, 82.27 ± 1.85% at 25 μg/mL and 79.94 ± 1.88% at 12.5 μg/mL) when compared to positive control acarbose (77.96 ± 2.06%, 71.17 ± 0.60%, 67.24 ± 1.16% at 50 μg/mL, 25 μg/mL and 12.5 μg/mL concentration). Molecular docking of the most potent enzyme inhibitor, 8f, shows promising interaction with the binding site of biological macromolecule Aspergillus oryzae α-amylase (PDB ID: 7TAA) and human pancreatic α-amylase (PDB ID: 3BAJ). To a step further, in-depth QSAR studies show a significant correlation between the experimental and the predicted inhibitory activities with the best R_validation²= 0.8701. The developed QSAR model can provide ample information about the structural features responsible for the increase and decrease of inhibitory activity. The mechanistic interpretation of the structure-activity relationship (SAR) is done with the help of combined computational calculations i.e. molecular docking and QSAR. Finally, molecular dynamic simulations are performed to get an insight into the binding mode of the most potent derivative with α-amylase from A. oryzae (PDB ID: 7TAA) and human pancreas (PDB ID: 3BAJ).

Quantitative structure activity relationship studies of novel hydrazone derivatives as α-amylase inhibitors with index of ideality of correlation

The present manuscript describes the synthesis, α-amylase inhibition, in silico studies and in-depth quantitative structure-activity relationship (QSAR) of a library of aroyl hydrazones based on benzothiazole skeleton. All the compounds of the developed library are characterized by various spectral techniques. α-Amylase inhibitory potential of all compounds has been explored, where compound 7n exhibits remarkable α-amylase inhibition of 87.5% at 50 µg/mL. Robust QSAR models are made by using the balance of correlation method in CORAL software. The chemical structures at different concentration with optimal descriptors are represented by SMILES. A data set of 66 SMILES of 22 hydrazones at three distinct concentrations are prepared. The significance of the index of ideality of correlation (IIC) with applicability domain (AD) is also studied at depth. A QSAR model with best Rvalidation2 = 0.8587 for split 1 is considered as a leading model. The outliers and promoters of increase and decrease of endpoint are also extracted. The binding modes of the most active compound, that is, 7n in the active site of Aspergillus oryzae α-amylase (PDB ID: 7TAA) are also explored by in silico molecular docking studies. Compound 7n displays high resemblance in binding mode and pose with the standard drug acarbose. Molecular dynamics simulations performed on protein-ligand complex for 100 ns, the protein gets stabilised after 20 ns and remained below 2 Å for the remaining simulation. Moreover, the deviation observed in RMSF during simulation for each amino acid residue with respect to Cα carbon atom is insignificant.

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 Monte Carlo Method as a Tool to Build up Predictive QSPR/QSAR

Background:

The Monte Carlo method has a wide application in various scientific researches. For the development of predictive models in a form of the quantitative structure-property / activity relationships (QSPRs/QSARs), the Monte Carlo approach also can be useful. The CORAL software provides the Monte Carlo calculations aimed to build up QSPR/QSAR models for different endpoints.

Methods:

Molecular descriptors are a mathematical function of so-called correlation weights of various molecular features. The numerical values of the correlation weights give the maximal value of a target function. The target function leads to a correlation between endpoint and optimal descriptor for the visible training set. The predictive potential of the model is estimated with the validation set, i.e. compounds that are not involved in the process of building up the model.

Results:

The approach gave quite good models for a large number of various physicochemical, biochemical, ecological, and medicinal endpoints. Bibliography and basic statistical characteristics of several CORAL models are collected in the present review. In addition, the extended version of the approach for more complex systems (nanomaterials and peptides), where behaviour of systems is defined by a group of conditions besides the molecular structure is demonstrated.

Conclusion:

The Monte Carlo technique available via the CORAL software can be a useful and convenient tool for the QSPR/QSAR analysis.

The using of the Index of Ideality of Correlation (IIC) to improve predictive potential of models of water solubility for pesticides

Models for water solubility of pesticides suggested in this manuscript are important data from point of view of ecologic engineering. The Index of Ideality of Correlation (IIC) of groups of quantitative structure-property relationships (QSPRs) for water solubility of pesticides related to the calibration sets was used to identify good in silico models. This comparison confirmed the high IIC set provides better statistical quality of the model for the validation set. Though there are large databases on solubility, the reliable prediction of the endpoint for new substances which are potential pesticides is an important ecologic task. Unfortunately, predictive models for various endpoints suffer overtraining, and the IIC serves to avoid or at least reduce this. Thus, the approach suggested has both theoretical and economic effects for ecology.

Nucleobase sequence based building up of reliable QSAR models with the index of ideality correlation using Monte Carlo method

This study describes in silico designing of aptamers against the influenza virus using Monte Carlo method. Aptamers are short, single-stranded oligonucleotides and these bind to an ample range of biologically important proteins which are related to many disease conditions. The affinities and specificities of aptamers are comparable to antibodies. In the medicinal chemistry, quantitative structure-activity relationship (QSAR) is an important skill which is used for drug design and development. To study the inhibitory activity of aptamers, we have developed QSAR models based on Monte Carlo method. The nucleobase sequence descriptors B_k, BB_k and BBB_k are used to generate the QSAR models. A number of statistical benchmarks together with index of ideality of correlation (IIC) is considered to validate the build QSAR models. Data set of 98 aptamers is divided into four random splits. The statistical criteria R² = 0.8711 and CCC = 0.9207 of the validation set of split 3 are best, so the build QSAR model of split 3 is the paramount model. The aptamer fragment responsible for the promotors of endpoint increase and decrease are also determined. These fragments are applied to design new nine aptamers from the lead aptamer APT01.Communicated by Ramaswamy H. Sarma.

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.

Prediction of Biochemical Endpoints by the CORAL Software: Prejudices, Paradoxes, and Results

Quantitative structure-activity relationships (QSARs) for prediction of toxicological endpoints built up with the CORAL software are discussed. Prejudices related to these QSAR models are listed. Possible ways to improve the software are discussed.

Heme Oxygenase Database (HemeOxDB) and QSAR Analysis of Isoform 1 Inhibitors

Due to increasing interest in the field of heme oxygenases (HOs), we built a ligand database called HemeOxDB that includes the entire set of known HO-1 and HO-2 inhibitors, resulting in more than 400 compounds. The HemeOxDB is available online at http://www.researchdsf.unict.it/hemeoxdb/, and having a robust search engine allows end users to build complex queries, sort tabulated results, and generate color-coded two- and three-dimensional graphs. This database will grow to be a tool for the design of potent and selective HO-1 or HO-2 inhibitors. We were also interested in virtually searching for alternative inhibitors, and, for the first time in the field of HOs, a quantitative structure-activity relationship (QSAR) model was built using half-maximal inhibitory concentration (IC₅₀ ) values of the whole set of known HO-1 inhibitors, taken from the HemeOxDB and employing the Monte Carlo technique. The statistical quality suggested that the model is robust and possesses desirable predictive potential. The screening of US Food and Drug Administration (FDA)-approved drugs, external to our dataset, suggested new predicted inhibitors, opening the way for replacing imidazole groups. The HemeOxDB and the QSAR model reported herein may help in prospectively identifying or repurposing new drugs with optimal structural attributes for HO enzyme inhibition.

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.

Quasi-SMILES for Nano-QSAR Prediction of Toxic Effect of Al2O3 Nanoparticles

The level of malondialdehyde (MDA) in wet tissue of different organs is utilized as a measure of toxic effect. The numerical data on the concentration of MDA in wet tissue of liver, kidneys, brain, and heart of rat is examined as the endpoint which are impacted by different dose (mg/kg), exposure time (3 and 14 days) and single oral treatment of aluminium nano-oxide (Al2O3) with 30 nm or 40 nm. An attempt to develop predictive model for this endpoint has been carried out in this work. SMILES is a traditional tool to represent molecular structure for QSPRs/QSARs. In contrast to traditional SMILES, so-called quasi-SMILES can be a tool to build up quantitative features – property / activity relationships (QFPRs/QFARs) for endpoints which are not defined by solely molecular structure, but by a group of physicochemical and/or biochemical conditions. The quasi-SMILES is the representation of the above eclectic conditions whereas the QFPR/QFAR are models of endpoints which are modified under impacts of these eclectic conditions.

Use of the Monte Carlo Method for OECD Principles-Guided QSAR Modeling of SIRT1 Inhibitors

SIRT1 inhibitors offer therapeutic potential for the treatment of a number of diseases including cancer and human immunodeficiency virus infection. A diverse series of 45 compounds with reported SIRT1 inhibitory activity has been employed for the development of quantitative structure-activity relationship (QSAR) models using the Monte Carlo optimization method. This method makes use of simplified molecular input line entry system notation of the molecular structure. The QSAR models were built up according to OECD principles. Three subsets of three splits were examined and validated by respective external sets. All the three described models have good statistical quality. The best model has the following statistical characteristics: R²  = 0.8350, Q²_test  = 0.7491 for the test set and R²  = 0.9655, Q²_ext  = 0.9261 for the validation set. In the mechanistic interpretation, structural attributes responsible for the endpoint increase and decrease are defined. Further, the design of some prospective SIRT1 inhibitors is also presented on the basis of these structural attributes.

Quantitative Structure-Activity Relationship for High Affinity 5-HT1A Receptor Ligands Based on Norm Indexes

Arylpiperazine derivatives are promising 5-hydroxytryptamine (5-HT) receptor ligands which can inhibit serotonin reuptake effectively. In this work, some norm index descriptors were proposed and further utilized to develop a model for predicting 5-HT1A receptor affinity (pKi) of 88 arylpiperazine derivatives. Results showed that this new model could provide satisfactory predictions with the square of the correction coefficient (R(2)) of 0.8891 and the squared correlation coefficient of cross-validation (Q(2)) of 0.8082, respectively. In addition, the applicability domain of this model was validated by using the leverage approach and results which suggested potential large scale for further utilization of this model. The results of statistical values and validation tests demonstrated that our proposed norm index based model could be successfully applied for predicting the affinity 5-HT1A receptor ligands of arylpiperazine derivatives.

QSAR as a random event: a case of NOAEL

Quantitative structure-activity relationships (QSAR) for no observed adverse effect levels (NOAEL, mmol/kg/day, in logarithmic units) are suggested. Simplified molecular input line entry systems (SMILES) were used for molecular structure representation. Monte Carlo method was used for one-variable models building up for three different splits into the "visible" training set and "invisible" validation. The statistical quality of the models for three random splits are the following: split 1 n = 180, r (2) = 0.718, q (2) = 0.712, s = 0.403, F = 454 (training set); n = 17, r (2) = 0.544, s = 0.367 (calibration set); n = 21, r (2) = 0.61, s = 0.44, r m (2) = 0.61 (validation set); split 2 n = 169, r (2) = 0.711, q (2) = 0.705, s = 0.409, F = 411 (training set); n = 27, r (2) = 0.512, s = 0.461 (calibration set); n = 22, r (2) = 0.669, s = 0.360, r m (2) = 0.63 (validation set); split 3 n = 172, r (2) = 0.679, q (2) = 0.672, s = 0.420, F = 360 (training set); n = 19, r (2) = 0.617, s = 0.582 (calibration set); n = 21, r (2) = 0.627, s = 0.367, r m (2) = 0.54 (validation set). All models are built according to OCED principles.

Monte Carlo QSAR models for predicting organophosphate inhibition of acetycholinesterase

A series of 278 organophosphate compounds acting as acetylcholinesterase inhibitors has been studied. The Monte Carlo method was used as a tool for building up one-variable quantitative structure-activity relationship (QSAR) models for acetylcholinesterase inhibition activity based on the principle that the target endpoint is treated as a random event. As an activity, bimolecular rate constants were used. The QSAR models were based on optimal descriptors obtained from Simplified Molecular Input-Line Entry System (SMILES) used for the representation of molecular structure. Two modelling approaches were examined: (1) 'classic' training-test system where the QSAR model was built with one random split into a training, test and validation set; and (2) the correlation balance based QSAR models were built with two random splits into a sub-training, calibration, test and validation set. The DModX method was used for defining the applicability domain. The obtained results suggest that studied activity can be determined with the application of QSAR models calculated with the Monte Carlo method since the statistical quality of all build models was very good. Finally, structural indicators for the increase and the decrease of the bimolecular rate constant are defined. The possibility of using these results for the computer-aided design of new organophosphate compounds is presented.

Quasi-QSAR for mutagenic potential of multi-walled carbon-nanotubes

Available on the Internet, the CORAL software (http://www.insilico.eu/coral) has been used to build up quasi-quantitative structure-activity relationships (quasi-QSAR) for prediction of mutagenic potential of multi-walled carbon-nanotubes (MWCNTs). In contrast with the previous models built up by CORAL which were based on representation of the molecular structure by simplified molecular input-line entry system (SMILES) the quasi-QSARs based on the representation of conditions (not on the molecular structure) such as concentration, presence (absence) S9 mix, the using (or without the using) of preincubation were encoded by so-called quasi-SMILES. The statistical characteristics of these models (quasi-QSARs) for three random splits into the visible training set and test set and invisible validation set are the following: (i) split 1: n=13, r(2)=0.8037, q(2)=0.7260, s=0.033, F=45 (training set); n=5, r(2)=0.9102, s=0.071 (test set); n=6, r(2)=0.7627, s=0.044 (validation set); (ii) split 2: n=13, r(2)=0.6446, q(2)=0.4733, s=0.045, F=20 (training set); n=5, r(2)=0.6785, s=0.054 (test set); n=6, r(2)=0.9593, s=0.032 (validation set); and (iii) n=14, r(2)=0.8087, q(2)=0.6975, s=0.026, F=51 (training set); n=5, r(2)=0.9453, s=0.074 (test set); n=5, r(2)=0.8951, s=0.052 (validation set).

Optimal descriptor as a translator of eclectic data into prediction of cytotoxicity for metal oxide nanoparticles under different conditions

The Monte Carlo technique has been used to build up quantitative structure-activity relationships (QSARs) for prediction of dark cytotoxicity and photo-induced cytotoxicity of metal oxide nanoparticles to bacteria Escherichia coli (minus logarithm of lethal concentration for 50% bacteria pLC50, LC50 in mol/L). The representation of nanoparticles include (i) in the case of the dark cytotoxicity a simplified molecular input-line entry system (SMILES), and (ii) in the case of photo-induced cytotoxicity a SMILES plus symbol '^'. The predictability of the approach is checked up with six random distributions of available data into the visible training and calibration sets, and invisible validation set. The statistical characteristics of these models are correlation coefficient 0.90-0.94 (training set) and 0.73-0.98 (validation set).

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.

Optimal nano-descriptors as translators of eclectic data into prediction of the cell membrane damage by means of nano metal-oxides

Systematization of knowledge on nanomaterials has become a necessity with the fast growth of applications of these species. Building up predictive models that describe properties (both beneficial and hazardous) of nanomaterials is vital for computational sciences. Classic quantitative structure-property/activity relationships (QSPR/QSAR) are not suitable for investigating nanomaterials because of the complexity of their molecular architecture. However, some characteristics such as size, concentration, and exposure time can influence endpoints (beneficial or hazardous) related to nanoparticles and they can therefore be involved in building a model. Application of the optimal descriptors calculated with the so-called correlation weights of various concentrations and different exposure times are suggested in order to build up a predictive model for cell membrane damage caused by a series of nano metal-oxides. The numerical data on correlation weights are calculated by the Monte Carlo method. The obtained results are in good agreement with the experimental data.

Optimal descriptor as a translator of eclectic information into the prediction of membrane damage: the case of a group of ZnO and TiO2 nanoparticles

The development of quantitative structure-activity relationships for nanomaterials needs representation of molecular structure of extremely complex molecular systems. Obviously, various characteristics of nanomaterial could impact associated biochemical endpoints. Following features of TiO2 and ZnO nanoparticles (n=42) are considered here: (i) engineered size (nm); (ii) size in water suspension (nm); (iii) size in phosphate buffered saline (PBS, nm); (iv) concentration (mg/L); and (v) zeta potential (mV). The damage to cellular membranes (units/L) is selected as an endpoint. Quantitative features-activity relationships (QFARs) are calculated by the Monte Carlo technique for three distributions of data representing values associated with membrane damage into the training and validation sets. The obtained models are characterized by the following average statistics: 0.78<r(2)training<0.92 and 0.67<r(2)validation<0.83.