CORAL: Quantitative structure-activity relationship models for estimating toxicity of organic compounds in rats

Development of a Quasi-Quantitative Structure-Activity Relationship Model for Prediction of the Immobilization Response of Daphnia magna Exposed to Metal-Based Nanomaterials

The conventional Hill equation model is suitable to fit dose-response data obtained from performing (eco)toxicity assays. Models based on quasi-quantitative structure-activity relationships (QSARs) to estimate the Hill coefficient ( n H ) ${n}_{{\rm{H}}})$ were developed with the aim of predicting the response of the invertebrate species Daphnia magna to exposure to metal-based nanomaterials. Descriptors representing the pristine properties of nanoparticles and media conditions were coded to a quasi-simplified molecular input line entry system and correlated to experimentally derived values of n H ${n}_{{\rm{H}}}$ . Monte Carlo optimization was used to model the set of n H ${n}_{{\rm{H}}}$  values, and the model was trained on the basis of reported dose-response relationships of 60 data sets (n = 367 individual response observations) of 11 metal-based nanomaterials as obtained from 20 literature reports. The model simulates the training data well, with only 2.3% deviation between experimental and modeled response data. The technique was employed to predict the dose-response relationships of 15 additional data sets (n = 72 individual observations) not included in model development of seven metal-based nanomaterials from 10 literature reports, with an average error of 3.5%. Combining the model output with either the median effective concentration value or any other known effect level as obtained from experimental data allows the prediction of full dose-response curves of D. magna immobilization. This model is an accurate screening tool that allows the determination of the shape and slope of dose-response curves, thereby greatly reducing experimental effort in case of novel advanced metal-based nanomaterials or the prediction of responses in altered exposure media. This screening model is compliant with the 3Rs (replacement, reduction, and refinement) principle, which is embraced by the scientific and regulatory communities dealing with nano-safety. Environ Toxicol Chem 2022;41:1439-1450. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Synergistic Behavior of Plant Proteins and Biobased Latexes in Bioplastic Food Packaging Materials: Experimental and Machine Learning Study

Plant-based proteins are attractive components which may serve as sustainable alternatives to current petrochemical products. Both soy protein and major corn protein, zein, are of interest in food packaging applications due to their sustainability, biodegradation properties, and inherent physicochemical properties. This study discusses the development of bioplastic materials, where it explores the effects of combining zein, soy protein, and plasticizing latexes derived from plant oil-based monomers (POBMs) on properties of resulting bioplastic films. By looking for synergistic effects of soy protein's inherent film formation ability and zein's higher strength, we prepare strong yet flexible soy-zein films as materials, called proteoposites. Incorporation of natural additive POBM-latexes helps to plasticize and hydrophobize the bioplastic films and thus to improve mechanical and barrier properties. Variation of the POBM-latexes' particle size further aims to enhance the performance of resulting bioplastic films. As a result, modified soy-zein proteoposite films with improved moisture resistance, enhanced mechanical behavior, and greater barrier properties were developed. Machine learning-based computational models were utilized in order to find main structural factors affecting the bioplastic's properties and develop a quantitative structure-property relationship model between the physicochemical properties of the film components and the resulted bioplastics' properties and performance. The developed model effectively predicts experimental outcomes with >85% (R²: 0.85) accuracy. The newly synthesized proteoposites confirmed the machine learning model predictions. As a result, proteoposite films made of two plant proteins and modified with POBM-latexes can be considered as an attractive and viable replacement for petrochemical food packaging products.

Correlation intensity index: Building up models for mutagenicity of silver nanoparticles

Nanomaterials become significant component of economics. Consequently, nanomaterials become object of environmental sciences. There is a traditional list of endpoints which are indicators of the ecological risk. Mutagenicity is one of important component in this list. The quasi-SMILES approach, that in contrast to majority of work dedicated to modelling behaviour of nanomaterials gives possibility to consider experimental conditions as well as other circumstances which can impact the behaviour of nanomaterials is suggested. This is carried out via so-called quasi-SMILES. The quasi-SMILES is a line on of codes that contains all the above available eclectic data. Modelling process aimed to build up a model involves Correlation Intensity Index (CII) that is a new criterion of predictive potential of models. The scheme of calculation of CII is described in this work in the first time. The applying of CII together with Index of Ideality Correlation (IIC) in modelling of mutagenicity of silver nanoparticles by the Monte Carlo method using the CORAL software (http://www.insilico.eu/coral) indicates that application of the CII improves the predictive potential of these models for three random splits into the training set (75%) and validation set (25%).

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 R² 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.

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.

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

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

QSAR Development for Plasma Protein Binding: Influence of the Ionization State

Purpose

This study explored several strategies to improve the performance of literature QSAR models for plasma protein binding (PPB), such as a suitable endpoint transformation, a correct representation of chemicals, more consistency in the dataset, and a reliable definition of the applicability domain.

Methods

We retrieved human fraction unbound (Fu) data for 670 compounds from the literature and carefully checked them for consistency. Descriptors were calculated taking account of the ionization state of molecules at physiological pH (7.4), in order to better estimate the affinity of molecules to blood proteins. We used different algorithms and chemical descriptors to explore the most suitable strategy for modeling the endpoint. SMILES (simplified molecular input line entry system)-based string descriptors were also tested with the CORAL software (CORelation And Logic). We did an outlier analysis to establish the models to use (or not to use) in case of well recognized families.

Results

Internal validation of the selected models returned Q² values close to 0.60. External validation also gave r² values always greater than 0.60. The CORAL descriptor based model for √fu was the best, with r² 0.74 in external validation.

Conclusions

Performance in prediction confirmed the robustness of all the derived models and their suitability for real-life purposes, i.e. screening chemicals for their ADMET profiling. Optimization of descriptors can be useful in order to obtain the correct results with a ionized molecule.

Electronic supplementary material

The online version of this article (10.1007/s11095-018-2561-8) contains supplementary material, which is available to authorized users.

Finding the structural requirements of diverse HIV-1 protease inhibitors using multiple QSAR modelling for lead identification

Multiple Quantitative Structure-Activity Relationship (QSAR) analysis is widely used in drug discovery for lead identification. Human Immunodeficiency Virus (HIV) protease is one of the key targets for the treatment of Acquired Immunodeficiency Syndrome (AIDS). One of the major challenges for the design of HIV-1 protease inhibitors (HIV PRIs) is to increase the inhibitory activities against the enzyme to a level where the problem associated to drug resistance may be considerably delayed. Herein, chemometric analyses were performed with 346 structurally diverse HIV PRIs with experimental bioactivities against a sub-type B mutant to develop highly predictable QSAR models and also to identify the effective structural determinants for higher affinity against HIV PR. The QSAR models were developed using OCHEM-based machine learning tools (ASNN, FSMLR, KNN, RF, MANN and XGBoost), with descriptors calculated by eight different software packages. Simultaneously, a Monte Carlo optimization-based QSAR modelling was performed using SMILES and graph-based descriptors to understand fragment and topochemical contributions. To validate the actual predictability of all these models, an additional set of 104 compounds (also with known experimental activities) with slightly different chemical space were employed. This ligand-based study serves as a crucial benchmark for further development of the HIV protease inhibitors with improved activities.

Use of Simplified Molecular Input Line Entry System and molecular graph based descriptors in prediction and design of pancreatic lipase inhibitors

Aim: The inhibition of pancreatic lipase (PL) enzyme is the most explored strategy for the treatment of obesity. The present study describes the development of quantitative structure–activity relationship (QSAR) models for a diverse set of 293 PL inhibitors by means of the Monte Carlo optimization technique. Methodology & results: The hybrid optimal descriptors were used to build QSAR models with three subsets of three splits. The developed QSAR models were further validated with corresponding external sets. The best QSAR model has the following statistical particulars: R2 = 0.752, Q LOO 2 = 0 . 736 for the test set and R2 = 0.768, Q F 1 2 = 0 . 628 , Q F 2 2 = 0 . 621 for the validation set. Conclusion: The developed QSAR models were robust, stable and predictive and led to the design of novel PL inhibitors.

CORAL: Building up QSAR models for the chromosome aberration test

A high level of chromosomal aberrations in peripheral blood lymphocytes may be an early marker of cancer risk, but data on risk of specific cancers and types of chromosomal aberrations are limited. Consequently, the development of predictive models for chromosomal aberrations test is important task. Majority of models for chromosomal aberrations test are so-called knowledge-based rules system. The CORAL software (http://www.insilico.eu/coral, abbreviation of “CORrelation And Logic”) is an alternative for knowledge-based rules system. In contrast to knowledge-based rules system, the CORAL software gives possibility to estimate the influence upon the predictive potential of a model of different molecular alerts as well as different splits into the training set and validation set. This possibility is not available for the approaches based on the knowledge-based rules system. Quantitative Structure–Activity Relationships (QSAR) for chromosome aberration test are established for five random splits into the training, calibration, and validation sets. The QSAR approach is based on representation of the molecular structure by simplified molecular input-line entry system (SMILES) without data on physicochemical and/or biochemical parameters. In spite of this limitation, the statistical quality of these models is quite good.

CORAL: QSAR models for carcinogenicity of organic compounds for male and female rats

Quantitative structure - activity relationships (QSARs) for carcinogenicity (rats, TD₅₀) have been built up using the CORAL software. Different molecular features, which are extracted from simplified molecular input-line entry system (SMILES) serve as the basis for building up a model. Correlation weights for the molecular features are calculated by means of the Monte Carlo optimization. Using the numerical data on the correlation weights, one can calculate a model of carcinogenicity as a mathematical function of descriptors, which are sum of the corresponding correlation weights. In other words, the correlation weights provide the maximal correlation coefficient between the descriptor and carcinogenicity, for the training set. This correlation was assessed via external validation set. Finally, lists of molecular alerts in aspects of carcinogenicity for male rats and for female rats were compared and their differences were discussed.

CORAL: Binary classifications (active/inactive) for drug-induced liver injury

INTRODUCTION

The data on human hepatotoxcity (drug-induced liver injury) is extremely important information from point of view of drug discovery. Experimental clinical data on this endpoint is scarce. Experimental way to extend databases on this endpoint is extremely difficult. Quantitative structure - activity relationships (QSAR) is attractive alternative of the experimental approach.

METHODS

Predictive models for human hepatotoxicity (drug-induced liver injury) have been built up by the Monte Carlo method with using of the CORAL software (http://www.insilico.eu/coral). These models are the binary classifications into active class and inactive class. These models are calculated with so-called "semi correlations" described in this work. The Mattews correlation coefficient of these models for external validation sets ranged from 0.52 to 0.62.

RESULTS DISCUSSION

The approach has been checked up with a group of random splits into the training and validation sets. These stochastic experiments have shown the stability of results: predictability of the models for various splits. Thus, the attempt to build up the classification QSAR model by means of the Monte Carlo technique, based on representation of the molecular structure via simplified molecular input line entry systems (SMILES) and hydrogen suppressed graph (HSG) using the CORAL software (http://www.insilico.eu/coral) has shown ability of this approach to provide quite good prediction of the examined endpoint (drug-induced liver injury).

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.

Monte Carlo-based quantitative structure-activity relationship models for toxicity of organic chemicals to Daphnia magna

Quantitative structure-activity relationships (QSARs) for toxicity of a large set of 758 organic compounds to Daphnia magna were built up. The simplified molecular input-line entry system (SMILES) was used to represent the molecular structure. The Correlation and Logic (CORAL) software was utilized as a tool to develop the QSAR models. These models are built up using the Monte Carlo method and according to the principle "QSAR is a random event" if one checks a group of random distributions in the visible training set and the invisible validation set. Three distributions of the data into the visible training, calibration, and invisible validation sets are examined. The predictive potentials (i.e., statistical characteristics for the invisible validation set of the best model) are as follows: n = 87, r²  = 0.8377, root mean square error = 0.564. The mechanistic interpretations and the domain of applicability of built models are suggested and discussed. Environ Toxicol Chem 2016;35:2691-2697. © 2016 SETAC.

Towards predicting the solubility of CO2 and N2 in different polymers using a quasi-SMILES based QSPR approach

The solubility of gases in various polymers plays an important role for the design of new polymeric materials. Quantitative structure-property relationship (QSPR) models were designed to predict the solubility of gases such as CO2 and N2 in polyethylene (PE), polypropylene (PP), polystyrene (PS), polyvinyl acetate (PVA) and poly (butylene succinate) (PBS) at different temperatures and pressures by using quasi-SMILES codes. The dataset of 315 systems was split randomly into training, calibration and validation sets; random split 1 led to 214 training (r(2) = 0.870 and RMSE = 0.019), 51 calibration (r(2) = 0.858 and RMSE = 0.020) and 50 validation (r(2) = 0.869 and RMSE = 0.017) sets. The suggested approach based on the quasi-SMILES, which are analogues of the traditional SMILES gives reasonable good predictions for solubility of CO2 and N2 in different polymers. The described methodology is universal for situations where the aim is to predict the response of an eclectic system upon a variety of physicochemical and/or biochemical conditions.

Building up a QSAR model for toxicity toward Tetrahymena pyriformis by the Monte Carlo method: A case of benzene derivatives

Data on toxicity toward Tetrahymena pyriformis is indicator of applicability of a substance in ecologic and pharmaceutical aspects. Quantitative structure-activity relationships (QSARs) between the molecular structure of benzene derivatives and toxicity toward T. pyriformis (expressed as the negative logarithms of the population growth inhibition dose, mmol/L) are established. The available data were randomly distributed three times into the visible training and calibration sets, and invisible validation sets. The statistical characteristics for the validation set are the following: r(2)=0.8179 and s=0.338 (first distribution); r(2)=0.8682 and s=0.341 (second distribution); r(2)=0.8435 and s=0.323 (third distribution). These models are built up using only information on the molecular structure: no data on physicochemical parameters, 3D features of the molecular structure and quantum mechanics descriptors are involved in the modeling process.

Nano-QSAR: Model of mutagenicity of fullerene as a mathematical function of different conditions

The experimental data on the bacterial reverse mutation test (under various conditions) on C60 nanoparticles for the cases (i) TA100, and (ii) WP2uvrA/pkM101 are examined as endpoints. By means of the optimal descriptors calculated with the Monte Carlo method a mathematical model of these endpoints has been built up. The models are a mathematical function of eclectic data such as (i) dose (g/plate); (ii) metabolic activation (i.e. with mix S9 or without mix S9); and (iii) illumination (i.e. darkness or irradiation). The eclectic data on different conditions were represented by so-called quasi-SMILES. In contrast to the traditional SMILES which are representation of molecular structure, the quasi-SMILES are representation of conditions by sequence of symbols. The calculations were carried out with the CORAL software, available on the Internet at http://www.insilico.eu/coral. The main idea of the suggested descriptors is the accumulation of all available eclectic information in the role of logical and digital basis for building up a model. The computational experiments have shown that the described approach can be a tool to build up models of mutagenicity of fullerene under different conditions.

Self-Assembled Binary Nanoscale Systems: Multioutput Model with LFER-Covariance Perturbation Theory and an Experimental-Computational Study of NaGDC-DDAB Micelles

Studies of the self-aggregation of binary systems are of both theoretical and practical importance. They provide an opportunity to investigate the influence of the molecular structure of the hydrophobe on the nonideality of mixing. On the other hand, linear free energy relationship (LFER) models, such as Hansch's equations, may be used to predict the properties of chemical compounds such as drugs or surfactants. However, the task becomes more difficult once we want to predict simultaneaously the effect over multiple output properties of binary systems of perturbations under multiple input experimental boundary conditions (b(j)). As a consequence, we need computational chemistry or chemoinformatics models that may help us to predict different properties of the autoaggregation process of mixed surfactants under multiple conditions. In this work, we have developed the first model that combines perturbation theory (PT) and LFER ideas. The model uses as input covariance PT operators (CPTOs). CPTOs are calculated as the difference between covariance ΔCov((i)μ(k)) functions before and after multiple perturbations in the binary system. In turn, covariances calculated as the product of two Box-Jenkins operators (BJO) operators. BJOs are used to measure the deviation of the structure of different chemical compounds from a set of molecules measured under a given subset of experimental conditions. The best CPT-LFER model found predicted the effects of 25,000 perturbations over 9 different properties of binary systems. We also reported experimental studies of different experimental properties of the binary system formed by sodium glycodeoxycholate and didodecyldimethylammonium bromide (NaGDC-DDAB). Last, we used our CPT-LFER model to carry out a 1000 data point simulation of the properties of the NaGDC-DDAB system under different conditions not studied experimentally.

Evaluation of QSAR models for predicting the partition coefficient (log P) of chemicals under the REACH regulation

The partition coefficient (log P) is a physicochemical parameter widely used in environmental and health sciences and is important in REACH and CLP regulations. In this regulatory context, the number of existing experimental data on log P is negligible compared to the number of chemicals for which it is necessary. There are many models to predict log P and we have selected a number of free programs to examine how they predict the log P of chemicals registered for REACH and to evaluate wheter they can be used in place of experimental data. Some results are good, especially if the information on the applicability domain of the models is considered, with R(2) values from 0.7 to 0.8 and root mean square error (RMSE) from 0.8 to 1.5.

CORAL: model for no observed adverse effect level (NOAEL)

The in vivo repeated dose toxicity (RDT) test is intended to provide information on the possible risk caused by repeated exposure to a substance over a limited period of time. The measure of the RDT is the no observed adverse effect level (NOAEL) that is the dose at which no effects are observed, i.e., this endpoint indicates the safety level for a substance. The need to replace in vivo tests, as required by some European Regulations (registration, evaluation authorization and restriction of chemicals) is leading to the searching for reliable alternative methods such as quantitative structure-activity relationships (QSAR). Considering the complexity of the RDT endpoint, for which data quality is limited and depends anyway on the study design, the development of QSAR for this endpoint is an attractive task. Starting from a dataset of 140 organic compounds with NOAEL values related to oral short term toxicity in rats, we developed a QSAR model based on optimal descriptors calculated with simplified molecular input-line entry systems and the graph of atomic orbitals by the Monte Carlo method, using CORAL software. Three different splits into the training, calibration, and validation sets are studied. The mechanistic interpretation of these models in terms of molecular fragment with positive or negative contributions to the endpoint is discussed. The probabilistic definition for the domain of applicability is suggested.

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.

Simplified molecular input line entry system-based: QSAR modelling for MAP kinase-interacting protein kinase (MNK1)

Quantitative structure-activity relationship (QSAR) models were built for the prediction of inhibition (pIC50, i.e. negative logarithm of the 50% effective concentration) of MAP kinase-interacting protein kinase (MNK1) by 43 potent inhibitors. The pIC50 values were modelled with five random splits, with the representations of the molecular structures by simplified molecular input line entry system (SMILES). QSAR model building was performed by the Monte Carlo optimisation using three methods: classic scheme; balance of correlations; and balance correlation with ideal slopes. The robustness of these models were checked by parameters as rm(2), r(*)m(2), [Formula: see text] and randomisation technique. The best QSAR model based on single optimal descriptors was applied to study in vitro structure-activity relationships of 6-(4-(2-(piperidin-1-yl) ethoxy) phenyl)-3-(pyridin-4-yl) pyrazolo [1,5-a] pyrimidine derivatives as a screening tool for the development of novel potent MNK1 inhibitors. The effects of alkyl group, -OH, -NO2, F, Cl, Br, I, etc. on the IC50 values towards the inhibition of MNK1 were also reported.

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.

Discovery of Potential, Non-Toxic Influenza Virus Inhibitor by Computational Techniques

Influenza infection continues to be a major problem in many parts of the world. Rimantadine is a first-line drug used to treat the influenza infection by targeting M2 proton channel. However, S31N mutation in M2 proton channel diminishes the efficiency of rimantadine and creates resistance. To address this issue, the present study was aimed to screen the effective lead candidate against drug resistance strain of influenza from DrugBank database. Initially, the lead molecules were filtered using Lipinski rule of five and the drug likeliness property. Subsequently, the data reduction was carried out by employing molecular docking study. Finally, molecular dynamics simulations techniques were performed to validate the lead compound. Most importantly, the -p LD50 of the screened lead molecule was calculated using CORAL software to estimate the Rat oral toxicity. Accordingly, memantine may possibly become a promising lead compound of rimantadine-resistant influenza virus strain.

Optimal descriptor as a translator of eclectic data into endpoint prediction: mutagenicity of fullerene as a mathematical function of conditions

The experimental data on the bacterial reverse mutation test on C60 nanoparticles (TA100) is examined as an endpoint. By means of the optimal descriptors calculated with the Monte Carlo method a mathematical model of the endpoint has been built up. The model is the mathematical function of (i) dose (g/plate); (ii) metabolic activation (i.e. with S9 mix or without S9 mix); and (iii) illumination (i.e. dark or irradiation). The statistical quality of the model is the following: n=10, r(2)=0.7549, q(2)=0.5709, s=7.67, F=25 (Training set); n=5, r(2)=0.8987, s=18.4 (Calibration set); and n=5, r(2)=0.6968, s=10.9 (Validation set).

QSPR modelling of dielectric constants of π-conjugated organic compounds by means of the CORAL software

Simplified molecular input line entry system (SMILES) notations of 116 π-conjugated organic compounds have been used in three random splits to develop single optimal descriptor based quantitative structure-property relationships (QSPR) models for the prediction of dielectric constants by CORAL (CORrelation And Logic). Four kinds of optimal descriptors were obtained based on SMILES, hydrogen suppressed graph (HSG), graph of atomic orbitals (GAO) and hybrid descriptors. The Monte Carlo optimization was carried out for each random split by three different methods: (i) classic scheme; (ii) balance of correlations; and (iii) balance of correlations with ideal slopes. The QSPR models gave reliable and accurate values of dielectric constant for all the π-conjugated organic compounds. SMILES and the hybrid-based QSPR model provided the best accuracy for the prediction of dielectric constants. Statistical characteristics of the QSPR model-1 based on classic scheme method are n = 110, r(2) = 0.860, Q(2) = 0.860, s = 1.84, MAE = 1.30 and F = 696 (training set), n = 6, r(2) = 0.947, Q(2) = 0.876, s = 0.955, MAE = 0.647 and F = 71 (test set). These QSPR models are further validated by an external validation set of 25 molecules and the robustness is checked by parameters like k, kk, rm(2), r(*)m(2), average rm(2), ∆rm(2)and randomization technique ([Formula: see text]).