Design and development of novel focal adhesion kinase (FAK) inhibitors using Monte Carlo method with index of ideality of correlation to validate QSAR

The coefficient of conformism of a correlative prediction (CCCP): Building up reliable nano-QSPRs/QSARs for endpoints of nanoparticles in different experimental conditions encoded via quasi-SMILES

Simulation of the physicochemical and biochemical behavior of nanomaterials has its own specifics. However, the main goal of modeling for both traditional substances and nanomaterials is the same. This is an ecologic risk assessment. The universal indicator of toxicity is the n-octanol/water partition coefficient. Mutagenicity indicates the possibility of future undesirable environmental effects, possibly greater than toxicity. Models have been proposed for the octanol/water distribution coefficient of gold nanoparticles and the mutagenicity of silver nanoparticles. Unlike the previous studies, here the models are built using an updated scheme, which includes two improvements. Firstly, the computing involves a new criterion for prediction potential, the so-called coefficient of conformism of a correlative prediction (CCCP); secondly, the Las Vegas algorithm is used to select the potentially most promising models from a group of models obtained by the Monte Carlo algorithm. Apparently, CCCP is a measure of the predictive potential (not only correlation). This can give an advantage in developing a model in comparison to using the classic determination coefficient. Likely, CCCP can be more informative than the classical determination coefficient. The Las Vegas algorithm is able to improve the model obtained by the Monte Carlo method.

Identification of potential FAK inhibitors using mol2vec molecular descriptor-based QSAR, molecular docking, ADMET study, and molecular dynamics simulation

This study aims to identify potential focal adhesion kinase (FAK) inhibitors through an integrated computational approach, combining mol2vec descriptor-based QSAR, molecular docking, ADMET study, and molecular dynamics simulation. A dataset of 437 compounds with known FAK inhibitory activities was used to develop QSAR models using machine learning algorithms combined with mol2vec descriptors. Subsequently, the most promising compounds were subjected to molecular docking against FAK to evaluate their binding affinities and key interactions. ADMET study and molecular dynamics simulation were also employed to investigate the pharmacokinetic, drug-like properties, and the stability of the protein-ligand complexes. The results showed that the mol2vec descriptor-based QSAR model established by support vector regression demonstrated good predictive performance (R2 = 0.813, RMSE = 0.453, MAE = 0.263 in case of training set, and R2 = 0.729, RMSE = 0.635, MAE = 0.477 in case of test set), indicating their reliability in identifying potent FAK inhibitors. Using this QSAR model and molecular docking, compound 21 (ZINC000004523722) was identified as the most potential compound, with predicted logIC50 value and binding energy of 2.59 and - 9.3 kcal/mol, respectively. The results of molecular dynamics simulation and ADMET study also further suggested its potential as a promising drug candidate. However, because our research was merely theoretical, additional in vitro and in vivo studies are required for the verification of these results.

Correlation intensity index-index of ideality of correlation: A hyphenated target function for furtherance of MAO-B inhibitory activity assessment

Monoamine oxidases are the enzymes involved in the management of brain homeostasis through oxidative deamination of monoamines such as neurotransmitters, tyramine etc. The excessive production of monoamine oxidase-B specifically results in numerous neurodegenerative disorders like Alzheimer's and Parkinson's diseases. Inhibitors of monoamine oxidase-B are applied in the management of these disorders. Here in this article we have developed robust hybrid descriptor based QSAR models related to 123 monoamine oxidase-B inhibitors through CORAL software by means of Monte Carlo optimization method. Three target functions were applied to prepare QSAR models and three splits were made for each target function. The most reliable, robust and better predictive QSAR models were developed with TF3 (correlation intensity index -index of ideality of correlation). Correlation intensity index showed positive effect on QSAR models. The structural features obtained from the QSAR modeling were incorporated in newly designed molecules and exhibited positive effect on their endpoint. Significant binding interactions were represented by these molecules in docking studies. Molecule B5 displayed prominent pIC50 (8.3) and binding affinity (-11.5 kcal mol-1) towards monoamine oxidase-B.

Semi-Correlations for Building Up a Simulation of Eye Irritation

The OECD recognizes that data on a compound's ability to treat eye irritation are essential for the assessment of new compounds on the market. In silico models are frequently used to provide information when experimental data are lacking. Semi-correlations, as they are called, can be useful to build up categorical models for eye irritation. Semi-correlations are latent regressions that can be used when the endpoint is expressed by two values: 1 for an active molecule and 0 for an inactive molecule. The regression line is based on the descriptor values which serve to distribute the data into four classes: true positive, true negative, false positive, and false negative. These values are applied to calculate the corresponding statistical criterion for assessing the predictive potential of the categorical model. In our model, the descriptor is the sum of what are termed correlation weights. These are defined by optimization using the Monte Carlo method. The target function of the optimization is related to the determination coefficient and the mean absolute error for the training set. Our model gives results that are better than those previously reported for the same endpoint.

Quasi-SMILES as a basis to build up models of endpoints for nanomaterials

Simplified molecular input-line entry system (SMILES) is a format for representing of the molecular structure. Quasi-SMILES is an extended format for representing molecular structure data and some eclectic data, which in principle could be applied to improve a model's predictive potential. Nano-quantitative structure-property relationships (nano-QSPRs) for energy gap (Eg, eV) of the metals oxide nanoparticles based on the quasi-SMILES give a predictive model for Eg, characterized by the following statistical quality for external validation set n = 22, R2 = 0.83, RMSE = 0.267.

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.

Quantitative structure activity relationship studies of androgen receptor binding affinity of endocrine disruptor chemicals with index of ideality of correlation, their molecular docking, molecular dynamics and ADME studies

Endocrine disrupter chemicals (EDCs) are both natural and man-made chemicals that mimic, block or interfere with human hormonal system. In the present manuscript, QSAR modeling was performed for the androgen disruptors that interfere with biosynthesis, metabolism or action of androgens that causes adverse effects on male reproductive system. A set of 96 EDCs that exhibited affinity towards androgen receptors (Log RBA) in rats were employed for carrying out QSAR studies using Hybrid descriptors (combination of HFG and SMILES) through Monte Carlo Optimization. Using index of ideality of correlation (TF2), five splits were formed and predictability of five models resulting from these splits was assessed by various validation parameters. Models resulted from first split was the top most one with R2validation = 0.7878. Structural attributes responsible for change in endpoint were studied by employing correlation weights of structural attributes. In order to further validate the model, new EDCs were designed using these attributes. In silico molecular modelling studies were performed to assess the detailed interactions with the receptor. The binding energies of all the designed compounds were observed to be better than lead and are in the range of -10.46 to -14.80. Molecular dynamics simulation of 100 ns was performed for ED01 and NED05. The results revealed that the protein-ligand complex bearing NED05 was more stable than lead ED01 exhibiting better interactions with the receptor. Further, in an attempt to assess their metabolism, ADME studies were evaluated using SwissADME. The developed model enables to predict the characteristics of designed compounds in an authentic way.Communicated by Ramaswamy H. Sarma.

Development of reliable quantitative structure-toxicity relationship models for toxicity prediction of benzene derivatives using semiempirical descriptors

The Health and environmental hazards of benzene and nitrobenzene (NB) derivatives have remained a topic of interest of researchers. In silico methods for prediction of toxicity of chemicals have proved their worth in accurate forecast of environmental as well as health toxicity and are strongly recommended by regulatory authorities. Two quantitative structure-toxicity relationship (QSTR) models explaining Scenedesmus obliquus toxicity trends among 39 benzene derivatives and Tetrahymena pyriformis toxicity of 103 NB and 392 benzene derivatives are developed using semiempirical quantum chemical parameters. The best constructed QSTR models have good fitting ability (R² = 0.8053, 0.7591, and 0.8283) and robustness (Q²_LOO = 0.7507, 0.7227, and 0.8194; Q²_LMO = 0.7338, 0.7153, and 0.8172). The external predictivity of all the models are quite good (R²_EXT = 0.8256, 0.9349, and 0.8698). Electronegativity, Cosmo volume, total energy, and molecular weight are responsible for the increase and decrease of toxicity of benzene derivatives against S. obliquus while electronegativity, electrophilicity index, the heat of formation, total energy, hydrophobicity, and cosmo volume are responsible for modulation of toxicity of NB and benzene derivatives toward T. pyriformis. These models fulfill the requirements of all the five OECD principles.

QSAR study of tetrahydropteridin derivatives as polo-like kinase 1(PLK1) Inhibitors with molecular docking and dynamics study

PLK1 is the key target for dealing with different cancer because it plays an important role in cell proliferation. According to the regulation of OECD, a QSAR model was developed from a dataset of 68 tetrahydropteridin derivatives. Three descriptors (maxHaaCH, ATSC7i, AATS7m) were considered for the development of the QSAR model. The reliability and predictability of the developed QSAR model were evaluated by various statistical parameters (r² = 0.8213, r²_ext = 0.8771 and CCC_ext = 0.9364). The maxHaaCH descriptor is positively correlated to pIC₅₀ whereas, the ATSC7i and AATS7m are negatively correlated with pIC₅₀. The QSAR model explains all the structural features and shows a good correlation with the activity. Based on molecular modelling techniques, five compounds (D1-D5) were designed. Molecular docking and dynamics studies of the most active compound were performed with PDB ID: 2RKU. The results of the present investigation may be employed to identify and develop effective inhibitors for the treatment of PLK1-related pathophysiological disorders.

Quantitative structure-activity relationship modeling for predication of inhibition potencies of imatinib derivatives using SMILES attributes

Chronic myelogenous leukemia (CML) which is resulted from the BCR-ABL tyrosine kinase (TK) chimeric oncoprotein, is a malignant clonal disorder of hematopoietic stem cells. Imatinib is used as an inhibitor of BCR-ABL TK in the treatment of CML patients. The main object of the present manuscript is focused on constructing quantitative activity relationships (QSARs) models for the prediction of inhibition potencies of a large series of imatinib derivatives against BCR-ABL TK. Herren, the inbuilt Monte Carlo algorithm of CORAL software is employed to develop QSAR models. The SMILES notations of chemical structures are used to compute the descriptor of correlation weights (CWs). QSAR models are established using the balance of correlation method with the index of ideality of correlation (IIC). The data set of 306 molecules is randomly divided into three splits. In QSAR modeling, the numerical value of R², Q², and IIC for the validation set of splits 1 to 3 are in the range of 0.7180-0.7755, 0.6891-0.7561, and 0.4431-0.8611 respectively. The numerical result of [Formula: see text] > 0.5 for all three constructed models in the Y-randomization test validate the reliability of established models. The promoters of increase/decrease for pIC₅₀ are recognized and used for the mechanistic interpretation of structural attributes.

CORAL: Development of a hybrid descriptor based QSTR model to predict the toxicity of dioxins and dioxin-like compounds with correlation intensity index and consensus modelling

In the present study, ninety-five halogenated dioxins and related chemicals (dibenzo-p-dioxins, dibenzofurans, biphenyls, and naphthalene) with endpoint pEC₅₀ were used to develop twelve quantitative structure toxicity relationship (QSTR) models using inbuilt Monte Carlo algorithm of CORAL software. The hybrid optimal descriptor of correlation weights (DCW) using a combination of SMILES and HSG (hydrogen suppressed graph) was employed to generate QSTR models. Three target functions i.e. TF₁ (W_IIC=W_CII=0), TF₂ (W_IIC= 0.3 & W_CII=0) and TF₃ (W_IIC= 0.0 &W_CII=0.3) were employed to develop robust QSTR models and the statistical outcomes of each target function were compared with each other. The correlation intensity index (CII) was found a reliable benchmark of the predictive potential for QSTR models. The numerical value of the determination coefficient of the validation set of split 1 computed by TF₃ was found highest (R_Valid²=0.8438). The fragments responsible for the toxicity of dioxins and related chemicals were also identified in terms of the promoter of increase/decrease for pEC₅₀. Three random splits (Split 1, Split 2 and Split 4) were selected for the extraction of the promoter of increase/decrease for pEC₅₀. In the last, consensus modelling was performed using the intelligent consensus tool of DTC lab (https://dtclab.webs.com/software-tools). The original consensus model, which was created by combining four distinct models employing the split 4 arrangement, was more predictive for the validation set and the numerical value of the determination coefficient of the test set (validation set) was increased from 0.8133 to 0.9725. For the validation set of split 4, the mean absolute error (MAE 100%) was also lowered from 0.513 to 0.2739.

Paradox of 'ideal correlations': improved model for air half-life of persistent organic pollutants

The persistence of organic pollutants is an important environmental property due to the extended possibility to have an impact of corresponding substances. In many cases, the experimental values of the thousands of contaminants are missing. The object of the study is novel computational modelling for air pollutions. Quantitative structure-property relationship (QSPR) for air half-life has been built using the Monte Carlo method with applying the index of ideality of correlation (IIC). The basis of the predictive model of air half-life is the representation of the molecular structure by simplifying molecular input-line entry system (SMILES) and numerical data on the above endpoint (expressed by hours) converted to a decimal logarithm. The statistical quality of the model has been checked up with different validation metrics and is quite good. Paradoxically, the improvement of the statistical quality via the IIC for the validation set is done in detriment to the training set. The new model has performed better than those obtained previously on the same set of compounds, for the prediction of new compounds in the validation set. Some semi-quantitative indicators for the mechanistic interpretation of the model are suggested.

Carcinogenicity prediction using the index of ideality of correlation

Carcinogenicity testing is necessary to protect human health and comply with regulations, but testing it with the traditionally used two-year rodent studies is time-consuming and expensive. In certain cases, such as for impurities, alternative methods may be convenient. Thus there is an urgent need for alternative approaches for reliable and robust assessments of carcinogenicity. The Monte Carlo technique with CORAL software is a tool to tackle this task for unknown compounds using available experimental data for a representative set of compounds. The models can be constructed with the simplified molecular input line entry system without additional physicochemical descriptors. We describe here a model based on a data set of 1167 substances. Matthew's correlation coefficient values for calibration and validation sets are 0.747 and 0.577, respectively. Double bonds between carbon atoms and double bonds of oxygen atoms are the molecular features that indicate the carcinogenic potential of a compound.

Monte Carlo Method and GA-MLR-Based QSAR Modeling of NS5A Inhibitors against the Hepatitis C Virus

Hepatitis C virus (HCV) is a serious disease that threatens human health. Despite consistent efforts to inhibit the virus, it has infected more than 58 million people, with 300,000 deaths per year. The HCV nonstructural protein NS5A plays a critical role in the viral life cycle, as it is a major contributor to the viral replication and assembly processes. Therefore, its importance is evident in all currently approved HCV combination treatments. The present study identifies new potential compounds for possible medical use against HCV using the quantitative structure–activity relationship (QSAR). In this context, a set of 36 NS5A inhibitors was used to build QSAR models using genetic algorithm multiple linear regression (GA-MLR) and Monte Carlo optimization and were implemented in the software CORAL. The Monte Carlo method was used to build QSAR models using SMILES-based optimal descriptors. Four splits were performed and 24 QSAR models were developed and verified through internal and external validation. The model created for split 3 produced a higher value of the determination coefficients using the validation set (R2 = 0.991 and Q2 = 0.943). In addition, this model provides interesting information about the structural features responsible for the increase and decrease of inhibitory activity, which were used to develop eight novel NS5A inhibitors. The constructed GA-MLR model with satisfactory statistical parameters (R2 = 0.915 and Q2 = 0.941) confirmed the predicted inhibitory activity for these compounds. The Absorption, Distribution, Metabolism, Elimination, and Toxicity (ADMET) predictions showed that the newly designed compounds were nontoxic and exhibited acceptable pharmacological properties. These results could accelerate the process of discovering new drugs against HCV.

The searching for agents for Alzheimer's disease treatment via the system of self-consistent models

Robust quantitative structure-activity relationships (QSARs) for hBACE-1 inhibitors (pIC50) for a large database (n = 1706) are established. New statistical criteria of the predictive potential of models are suggested and tested. These criteria are the index of ideality of correlation (IIC) and the correlation intensity index (CII). The system of self-consistent models is a new approach to validate the predictive potential of QSAR-models. The statistical quality of models obtained using the CORAL software (http://www.insilico.eu/coral) for the validation sets is characterized by the average determination coefficient R²_v= 0.923, and RMSE =0.345. Three new promising molecular structures which can become inhibitors hBACE-1 are suggested.

CORAL: Monte Carlo based global QSAR modelling of Bruton tyrosine kinase inhibitors using hybrid descriptors

Global QSAR modelling was performed to predict the pIC₅₀ values of 233 diverse heterocyclic compounds as BTK inhibitors with the Monte Carlo algorithm of CORAL software using the DCW hybrid descriptors extracted from SMILES notations of molecules. The dataset of 233 BTK inhibitors was randomly split into training, invisible training, calibration and validation sets. The index of ideality of correlation was also applied to build and judge the predictability of the QSAR models. Eight global QSAR models based on the hybrid optimal descriptor using two target functions, i.e. TF₁ (W_IIC = 0) and TF₂ (W_IIC = 0.2) have been constructed. The statistical parameters of QSAR models computed by TF₂ are more reliable and robust and were used to predict the pIC₅₀ values. The model constructed for split 4 via TF₂ is regarded as the best model and the numerical values of r²_Train, r²_Valid, Q²_Train and Q²_Valid are equal to 0.7981, 0.7429, 0.7898 and 0.6784, respectively. By internal and external validation techniques, the predictability and reliability of the designed models have been assessed. The structural attributes responsible for the increase and decrease of pIC₅₀ of BTK inhibitors were also identified.

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).

Prediction of power conversion efficiency of phenothiazine-based dye-sensitized solar cells using Monte Carlo method with index of ideality of correlation

Simplified molecular-input line-entry system (SMILES) notation and inbuilt Monte Carlo algorithm of CORAL software were employed to construct generative and prediction QSPR models for the analysis of the power conversion efficiency (PCE) of 215 phenothiazine derivatives. The dataset was divided into four splits and each split was further divided into four sets. A hybrid descriptor, a combination of SMILES and hydrogen suppressed graph (HSG), was employed to build reliable and robust QSPR models. The role of the index of ideality of correlation (IIC) was also studied in depth. We performed a comparative study to predict PCE using two target functions (TF₁ without IIC and TF₂ with IIC). Eight QSPR models were developed and the models developed with TF₂ was shown robust and reliable. The QSPR model generated from split 4 was considered a leading model. The different statistical benchmarks were computed for the lead model and these were rtraining set2=0.7784; rinvisible training set2=0.7955; rcalibration set2=0.7738; rvalidation set2=0.7506; Qtraining set2=0.7691; Qinvisible training set2=0.7850; Qcalibration set2=0.7501;  Qvalidation set2=0.7085; IIC_{training set} = 0.8590; IIC_{invisible training set} = 0.8297; IIC_{calibration set} = 0.8796; IIC_{validation set} = 0.8293, etc. The promoters of increase and decrease of endpoint PCE were also extracted.

The Monte Carlo method to build up models of the hydrolysis half-lives of organic compounds

The hydrolysis of organic chemicals such as pesticides, pollutants, or drugs can affect the fate and behaviour of environmental contaminants, so it is of interest to evaluate the stability of substances in water for various purposes. For the registration of organic compounds in Europe, information on hydrolysis must be presented. However, the experimental measurements of all chemicals would require enormous resources, and computational models may become attractive. Applying the CORAL software (http://www.insilico.eu/coral) quantitative structure-property relationships (QSPRs) were built up to model hydrolysis. The 2D-optimal descriptor is calculated with so-called correlation weights for attributes of simplified molecular input-line entry systems (SMILES). The correlation weights are obtained as results of the special Monte Carlo optimization. The nature of (five- and six-member) rings is an important component of this approach. Another important component is the atom pair proportions for nitrogen, oxygen, and sulphur. The statistical quality of the best model is: n = 44, r² = 0.74 (training set); n = 14, r² = 0.75 (calibration set); and n = 12, r² = 0.80 (validation set).

The index of ideality of correlation improves the predictive potential of models of the antioxidant activity of tripeptides from frog skin (Litoria rubella)

It is usually held that good-quality models for the biological activity of peptides must take into account their 3D architecture and descriptors of quantum mechanics. However, the present study shows that it is possible to build up models without these complex calculations. The structure of tripeptides represented by sequences of one-symbol abbreviations of the corresponding amino acids serves to build up quantitative structure-activity relationships for the antioxidant activity of tripeptides from frog skin. The statistical quality of the best model for the validation set is n = 27, r² = 0.93, RMSE = 0.15.

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.

A Monte Carlo method based QSPR model for prediction of reaction rate constants of hydrated electrons with organic contaminants

The Monte Carlo algorithm was applied to formulate a robust quantitative structure-property relationship (QSPR) model to compute the reactions rate constants of hydrated electron values for a data set of 309 water contaminants containing 125 aliphatic and 184 phenyl-based chemicals. The QSPR models were computed with the hybrid optimal descriptors which were procured by combining the SMILES and hydrogen-suppressed molecular graph for both classes of compounds. Approximately 75% of the total experimental data set was randomly divided into training and invisible training sets, while approximately 25% was divided into calibration and validation sets. The authenticity and robustness of the developed QSPR models were also judged by the Index of Ideality of Correlation. In QSPR modelling of aliphatic compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.852-0.905, 0.815-0.894, 0.839-0.897 and 0.737-0.867, respectively. Whereas, in the QSPR modelling of phenyl-based compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.867-0.896, 0.852-0.865, 0.816-0.850 and 0.760-0.762, respectively. The structural attributes, which are promoters of l o g K e a q - increase/decrease are also extracted from the SMILES notation for mechanistic interpretation. These QSPR models can also be applied to compute the reaction rate constants of organic contaminants.

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%).

QSAR modelling of larvicidal phytocompounds against Aedes aegypti using index of ideality of correlation

Aedes aegypti is the primary vector of several infectious viruses that cause yellow, dengue, chikungunya, and Zika fevers. Recently, plant-derived products have been tested as safe and eco-friendly larvicides against Ae. aegypti. The present study aimed to improve QSAR models for 62 larvicidal phytocompounds against Ae. aegypti via the Monte Carlo method based on the index of the ideality of correlation (IIC) criterion. The representation of structures was done with SMILES. Three splits were prepared randomly and three QSAR models were constructed using IIC target function. The molecular descriptors were selected from SMILES descriptors and the hydrogen-filled molecular graphs. The predictability of three models was evaluated on the validation sets, the r ² of which was 0.9770, 0.8660, and 0.8565 for models 1 to 3, respectively. The statistical results of three randomized splits indicated that robust, simple, predictive, and reliable models were obtained for different sets. From the modelling results, important descriptors were identified to enhance and reduce the larvicidal activity of compounds. Based on the identified important descriptors, some new structures of larvicidal compounds were proposed. The larvicidal activity of novel molecules designed further was supported by docking studies. Using the simple QSAR model, one can predict pLC₅₀ of new similarity larvicidal phytocompounds.

In silico enhancement of azo dye adsorption affinity for cellulose fibre through mechanistic interpretation under guidance of QSPR models using Monte Carlo method with index of ideality correlation

Azo dyes are a group of chemical moieties joined by azo (-N=N-) group with potential usefulness in different industrial applications. But these dyes are not devoid of hazardous consequence because of poor affinity for the fibre and discharge into the water stream. The chemical aspects of 72 azo dyes towards cellulose fibre in terms of their affinity by QSPR have been explored in the present work. We have employed two approaches, namely balance of correlation without IIC (TF₁) and balance of correlation with IIC (TF₂), to generate 16 QSAR models from 8 splits. The determination coefficient of calibration and validation set was found higher when the QSPR models were developed using the index of ideality correlation (IIC) parameter (TF₂). The model developed with TF₂ for split 3 was considered as a prominent model because the determination coefficient of the validation set was maximum (r ² = 0.9468). The applicability domain (AD) was also analysed based on 'statistical defect', d(A) for a SMILES attribute. The mechanistic interpretation was done by identifying the SMILES attributes responsible for the promoter of endpoint increase and promoter of endpoint decrease. These SMILES attributes were applied to design 15 new dyes with higher affinity for cellulose fibre.

'Ideal correlations' for the predictive toxicity to Tetrahymena pyriformis

OBJECTIVES

Predictive models for toxicity to Tetrahymena pyriformis are an important component of natural sciences. The present study aims to build up a predictive model for the endpoint using the so-called index of ideality of correlation (IIC). Besides, the comparison of the predictive potential of these models with the predictive potential of models suggested in the literature is the task of the present study.

METHODS

The Monte Carlo technique is a tool to build up the predictive model applied in this study. The molecular structure is represented via a simplified molecular input-line entry system (SMILES). The IIC is a statistical characteristic sensitive to both the correlation coefficient and mean absolute error. Applying of the IIC to build up quantitative structure-activity relationships (QSARs) for the toxicity to Tetrahymena pyriformis improves the predictive potential of those models for random splits into the training set and the validation set. The calculation was carried out with CORAL software (http://www.insilico.eu/coral).

RESULTS

The statistical quality of the suggested models is incredibly good for the external validation set, but the statistical quality of the models for the training set is modest. This is the paradox of ideal correlation, which is obtained with applying the IIC.

CONCLUSIONS

The Monte Carlo technique is a convenient and reliable way to build up a predictive model for toxicity to Tetrahymena pyriformis. The IIC is a useful statistical criterion for building up predictive models as well as for the assessment of their statistical quality.

Identification of structural fingerprints for ABCG2 inhibition by using Monte Carlo optimization, Bayesian classification, and structural and physicochemical interpretation (SPCI) analysis

The human breast cancer resistance protein (BCRP), one of the members of the large ATP binding cassette (ABC) transporter superfamily, is crucial for resistance against chemotherapeutic agents. Currently, it has been emerged as one of the best biological targets for the designing of small molecule drugs capable of eliminating multidrug resistance in breast cancer. In order to gain insights into the relationship between the molecular structure of compounds and the ABCG2 inhibition, a multi-QSAR approach using different methods was performed on a dataset of 294 ABCG2 inhibitors with diverse scaffolds. The best models obtained by different chemometric methods have the following statistical characteristics: Monte Carlo Optimization-based QSAR (sensitivity = 0.905, specificity = 0.6255, accuracy = 0.756, and MCC = 0.545), Bayesian classification model (sensitivity = 0.735, specificity = 0.775, and concordance = 0.757); structural and physicochemical interpretation analysis-random forest method (balance accuracy = 0.750, sensitivity = 0.810, and specificity = 0.700). Additionally, structural fingerprints modulating the ABCG2 inhibitory properties were identified from the best models of each method and also validated with each other. The current modelling study is an attempt to get a deep insight into the different important structural fingerprints modulating ABCG2 inhibition.

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.

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.

Predictive QSAR modeling for the antioxidant activity of natural compounds derivatives based on Monte Carlo method

In this research, QSAR modeling was carried out through SMILES of compounds and on the basis of the Monte Carlo method to predict the antioxidant activity of 79 derivatives of pulvinic acid, 23 of coumarine, as well as nine structurally non-related compounds against three radiation sources of Fenton, gamma, and UV. QSAR model was designed through CORAL software, as well as a newer optimizing method well known as the index of ideality correlation. The full set of antioxidant compounds were randomly distributed into four sets, including training, invisible training, validation, and calibration; this division was repeated three times randomly. The optimal descriptors were picked up from a hybrid model by the combination of the hydrogen-suppressed graph and SMILES descriptors based on the objective function. These models' predictability was assessed on the sets of validation. The results of three randomized sets showed that simple, robust, reliable, and predictive models were achieved for training, invisible training, validation, and calibration sets of all three models. The central decrease/increase descriptors were identified. This simple QSAR can be useful to predict antioxidant activity of numerous antioxidants.

The index of ideality of correlation: models of the flash points of ternary mixtures

Reliable information related to the flash point of ternary mixtures assists in the rational classification of different ternary mixtures of liquids.

Whether the Validation of the Predictive Potential of Toxicity Models is a Solved Task?

Different kinds of biological activities are defined by complex biochemical interactions, which are termed as a "mathematical function" not only of the molecular structure but also for some additional circumstances, such as physicochemical conditions, interactions via energy and information effects between a substance and organisms, organs, cells. These circumstances lead to the great complexity of prediction for biochemical endpoints, since all "details" of corresponding phenomena are practically unavailable for the accurate registration and analysis. Researchers have not a possibility to carry out and analyse all possible ways of the biochemical interactions, which define toxicological or therapeutically attractive effects via direct experiment. Consequently, a compromise, i.e. the development of predictive models of the above phenomena, becomes necessary. However, the estimation of the predictive potential of these models remains a task that is solved only partially. This mini-review presents a collection of attempts to be used for the above-mentioned task, two special statistical indices are proposed, which may be a measure of the predictive potential of models. These indices are (i) Index of Ideality of Correlation; and (ii) Correlation Contradiction Index.

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.

The Correlation Contradictions Index (CCI): Building up reliable models of mutagenic potential of silver nanoparticles under different conditions using quasi-SMILES

The interpretation of the mutagenic potential of silver nanoparticles as a mathematical function of (i) dose; (ii) coating; and (iii) type of mutagenicity (TA98 and TA100) gives quantitative models with good statistical quality. So-called quasi-SMILES are used to represent examined objects (silver nanoparticles under different conditions) for building up models. Simplified molecular input-line entry systems (SMILES) is a well-known sequence of symbols for representation of the molecular structure. Quasi-SMILES is a similar sequence of symbols for representation of experimental conditions. The Correlation Contradiction Index (CCI) calculated with data on the calibration set gives possibility to predict quality of correlation of "experimental vs. calculated values of endpoint" for external validation set.

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.