The index of ideality of correlation: A criterion of predictability of QSAR models for skin permeability?

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.

The system of self-consistent models based on quasi-SMILES as a tool to predict the potential of nano-inhibitors of human lung carcinoma cell line A549 for different experimental conditions

The different features of the impact of nanoparticles on cells, such as the structure of the core, presence/absence of doping, quality of surface, diameter, and dose, were used to define quasi-SMILES, a line of symbols encoded the above physicochemical features of the impact of nanoparticles. The correlation weight for each code in the quasi-SMILES has been calculated by the Monte Carlo method. The descriptor, which is the sum of the correlation weights, is the basis for a one-variable model of the biological activity of nano-inhibitors of human lung carcinoma cell line A549. The system of models obtained by the above scheme was checked on the self-consistence, i.e., reproducing the statistical quality of these models observed for different distributions of available nanomaterials into the training and validation sets. The computational experiments confirm the excellent potential of the approach as a tool to predict the impact of nanomaterials under different experimental conditions. In conclusion, our model is a self-consistent model system that provides a user to assess the reliability of the statistical quality of the used approach.

Does the accounting of the local symmetry fragments in quasi-SMILES improve the predictive potential of the QSAR models of toxicity toward tadpoles?

Models of toxicity to tadpoles have been developed as single parameters based on special descriptors which are sums of correlation weights, molecular features, and experimental conditions. This information is presented by quasi-SMILES. Fragments of local symmetry (FLS) are involved in the development of the model and the use of FLS correlation weights improves their predictive potential. In addition, the index of ideality correlation (IIC) and correlation intensity index (CII) are compared. These two potential predictive criteria were tested in models built through Monte Carlo optimization. The CII was more effective than IIC for the models considered here.

Monte Carlo Optimization Method Based QSAR Modeling of Placental Barrier Permeability

PURPOSE

In order to ensure that drug administration is safe during pregnancy, it is crucial to have the possibility to predict the placental permeability of drugs in humans. The experimental method which is most widely used for the said purpose is in vitro human placental perfusion, though the approach is highly expensive and time consuming. Quantitative structure-activity relationship (QSAR) modeling represents a powerful tool for the assessment of the drug placental transfer, and can be successfully employed to be an alternative in in vitro experiments.

METHODS

The conformation-independent QSAR models covered in the present study were developed through the use of the SMILES notation descriptors and local molecular graph invariants. What is more, the Monte Carlo optimization method, was used in the test sets and the training sets as the model developer with three independent molecular splits.

RESULTS

A range of different statistical parameters was used to validate the developed QSAR model, including the standard error of estimation, mean absolute error, root-mean-square error (RMSE), correlation coefficient, cross-validated correlation coefficient, Fisher ratio, MAE-based metrics and the correlation ideality index. Once the mentioned statistical methods were employed, an excellent predictive potential and robustness of the developed QSAR model was demonstrated. In addition, the molecular fragments, which are derived from the SMILES notation descriptors accounting for the decrease or increase in the investigated activity, were revealed.

CONCLUSION

The presented QSAR modeling can be an invaluable tool for the high-throughput screening of the placental permeability of drugs.

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.

Using the Correlation Intensity Index to Build a Model of Cardiotoxicity of Piperidine Derivatives

The assessment of cardiotoxicity is a persistent problem in medicinal chemistry. Quantitative structure-activity relationships (QSAR) are one possible way to build up models for cardiotoxicity. Here, we describe the results obtained with the Monte Carlo technique to develop hybrid optimal descriptors correlated with cardiotoxicity. The predictive potential of the cardiotoxicity models (pIC50, Ki in nM) of piperidine derivatives obtained using this approach provided quite good determination coefficients for the external validation set, in the range of 0.90-0.94. The results were best when applying the so-called correlation intensity index, which improves the predictive potential of a model.

The system of self-consistent models for pesticide toxicity to Daphnia magna

Quantitative structure-property/activity relationships (QSPRs/QSARs) are a tool of modern theoretical and computational chemistry. The self-consistent model system is both a method to build up a group of QSPR/QSAR models and an approach to checking the reliability of these models. Here, a group of models of pesticide toxicity toward Daphnia magna for different distributions into training and test sub-sets is compared. This comparison is the basis for formulating the system of self-consistent models. The so-called index of the ideality of correlation (IIC) has been used to improve the above models' predictive potential of pesticide toxicity. The predictive potential of the suggested models should be classified as high since the average value of the determination coefficient for the validation sets is 0.841, and the dispersion is 0.033 (on all five models). The best model (number 4) has an average determination coefficient of 0.89 for the external validation sets (related to all five splits).

The validation of predictive potential via the system of self-consistent models: the simulation of blood-brain barrier permeation of organic compounds

CONTEXT

To apply the quantitative relationships "structure-endpoint" approach, the reliability of prediction is necessary but sometimes challenging to achieve. In this work, an attempt is made to accomplish the reliability of forecasts by creating a set of random partitions of data into training and validation sets, followed by constructing random models. A system of random models for a helpful approach should be self-consistent, giving a similar or at least comparable statistical quality of the predictions for models obtained using different splits of available data into training and validation sets.

METHOD

The carried out computer experiments aimed at obtaining blood-brain barrier permeation models showed that, in principle, can be used such an approach (the Monte Carlo optimization of the correlation weights for different molecular features) for the above purpose taking advantage of specific algorithms to optimize the modelling steps with applying of new statistical criteria such as the index of ideality of correlation (IIC) and the correlation intensity index (CII). The results so obtained are good and better than what was reported previously. The suggested approach to validation of models is non-identic to traditionally applied manners of the checking up models. The concept of validation can be used for arbitrary models (not only for models of the blood-brain barrier).

The enhancement scheme for the predictive ability of QSAR: A case of mutagenicity

Mutagenicity is one of the most dangerous properties from the point of view of medicine and ecology. Experimental determination of mutagenicity remains a costly process, which makes it attractive to identify new hazardous compounds based on available experimental data through in silico methods or quantitative structure-activity relationships (QSAR). A system for constructing groups of random models is proposed for comparing various molecular features extracted from SMILES and graphs. For mutagenicity (mutagenicity values were expressed by the logarithm of the number of revertants per nanomole assayed by Salmonella typhimurium TA98-S9 microsomal preparation) models, the Morgan connectivity values are more informative than the comparison of quality for different rings in molecules. The resulting models were tested with the previously proposed model self-consistency system. The average value of the determination coefficient for the validation set is 0.8737 ± 0.0312.

CORAL: Model of Ecological Impact of Heavy Metals on Soils via the Study of Modification of Concentration of Biomolecules in Earthworms (Eisenia fetida)

The traditional application for quantitative structure-property/activity relationships (QSPRs/QSARs) in the fields of thermodynamics, toxicology or drug design is predicting the impact of molecular features using data on the measurable characteristics of substances. However, it is often necessary to evaluate the influence of various exposure conditions and environmental factors, besides the molecular structure. Different enzyme-driven processes lead to the accumulation of metal ions by the worms. Heavy metals are sequestered in these organisms without being released back into the soil. In this study, we propose a novel approach for modeling the absorption of heavy metals, such as mercury and cobalt by worms. The models are based on optimal descriptors calculated for the so-called quasi-SMILES, which incorporate strings of codes reflecting experimental conditions. We modeled the impact on the levels of proteins, hydrocarbons, and lipids in an earthworm's body caused by different combinations of concentrations of heavy metals and exposure time observed over two months of exposure with a measurement interval of 15 days.

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 Novel Therapeutics for Schizophrenia Treatment Based on a Selective Positive Allosteric Modulation of α1-Containing GABAARs-In Silico Approach

For the development of atypical antipsychotics, the selective positive allosteric modulation of the ionotropic GABAA receptor (GABAAR) has emerged as a promising approach. In the presented research, two unrelated methods were used for the development of QSAR models for selective positive allosteric modulation of 1-containing GABAARs with derivatives of imidazo [1,2-a]-pyridine. The development of conformation-independent QSAR models, based on descriptors derived from local molecular graph invariants and SMILES notation, was achieved with the Monte Carlo optimization method. From the vast pool of 0D, 1D, and 2D molecule descriptors, the GA-MLR method developed additional QSAR models. Various statistical methods were utilised for the determination of the developed models’ robustness, predictability, and overall quality, and according to the obtained results, all QSAR models are considered good. The molecular fragments that have a positive or negative impact on the studied activity were obtained from the studied molecules’ SMILES notations, and according to the obtained results, nine novel compounds were designed. The binding affinities to GABAAR of designed compounds were assessed with the application of molecular docking studies and the obtained results showed a high correlation with results obtained from QSAR modeling. To assess all designed molecules’ “drug-likeness”, their physicochemical descriptors were computed and utilised for the prediction of medicinal chemistry friendliness, pharmacokinetic properties, ADME parameters, and druglike nature.

A regression-based QSAR-model to predict acute toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica): Calibration, validation, and future developments to support risk assessment of chemicals in amphibians

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

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.

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.

Development of novel antipsychotic agents by inhibiting dopamine transporter – in silico approach

Various in silico methods were employed for the development of antipsychotic agents by dopamine transporter inhibition.

Monte Carlo optimization method based QSAR modeling of postmortem redistribution of structurally diverse drugs

The Monte Carlo optimization method was employed for the development of the QSAR model for the prediction for postmortem redistribution of structurally diverse drugs.

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.

The Monte Carlo approach to model and predict the melting point of imidazolium ionic liquids using hybrid optimal descriptors

Ionic liquids (ILs) have captured intensive attention owing to their unique properties such as high thermal stability, negligible vapour pressure, high dissolution capacity and high ionic conductivity as well as their wide applications in various scientific fields including organic synthesis, catalysis, and industrial extraction processes. Many applications of ionic liquids (ILs) rely on the melting point (T_m). Therefore, in the present manuscript, the melting points of imidazolium ILs are studied employing a quantitative structure–property relationship (QSPR) approach to develop a model for predicting the melting points of a data set of imidazolium ILs. The Monte Carlo algorithm of CORAL software is applied to build up a robust QSPR model to calculate the values T_m of 353 imidazolium ILs. Using a combination of SMILES and hydrogen-suppressed molecular graphs (HSGs), the hybrid optimal descriptor is computed and used to generate the QSPR models. Internal and external validation parameters are also employed to evaluate the predictability and reliability of the QSPR model. Four splits are prepared from the dataset and each split is randomly distributed into four sets i.e. training set (≈33%), invisible training set (≈31%), calibration set (≈16%) and validation set (≈20%). In QSPR modelling, the numerical values of various statistical features of the validation sets such as R_Validation², Q_Validation², and IIC_Validation are found to be in the range of 0.7846–0.8535, 0.7687–0.8423 and 0.7424–0.8982, respectively. For mechanistic interpretation, the structural attributes which are responsible for the increase/decrease of T_m are also extracted.

The melting points of imidazolium ILs are studied employing a quantitative structure–property relationship (QSPR) approach to develop a model for predicting the melting points of a data set of imidazolium ILs.

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

Can the Monte Carlo method predict the toxicity of binary mixtures?

Risk assessment of toxicants mainly is a result of experiments with single substances. However, toxicity in natural ecosystems typically does not result from single toxicant exposure but is rather a result of exposure to mixtures of toxicants. It is not surprising a mixture of toxicity is a subject of eco-toxicological interest for several decades. A quantitative structure-activity relationships (QSAR)-based approach is an attractive approach to assessing the joint effects in the binary mixtures. The validity of the proposed approach was demonstrated by comparing the predicted values against the experimentally determined values. Simplified molecular input-line entry system (SMILES) is used for the representation of the molecular structures of components of two-component mixtures to build up QSAR. The SMILES-based models are improving if the Monte Carlo optimization aimed to define 2D-optimal descriptors apply the so-called index of ideality of correlation (IIC), which is a mathematical function of both the correlation coefficient and mean absolute error calculated for the positive and negative difference between observed and calculated values of toxicity. The average statistical quality of these models (for the validation set) is n=25, R²=0.95, and RMSE=0.375.

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

Inverse Mechanistic Modeling of Transdermal Drug Delivery for Fast Identification of Optimal Model Parameters

Transdermal drug delivery systems are a key technology to administer drugs with a high first-pass effect in a non-invasive and controlled way. Physics-based modeling and simulation are on their way to become a cornerstone in the engineering of these healthcare devices since it provides a unique complementarity to experimental data and additional insights. Simulations enable to virtually probe the drug transport inside the skin at each point in time and space. However, the tedious experimental or numerical determination of material properties currently forms a bottleneck in the modeling workflow. We show that multiparameter inverse modeling to determine the drug diffusion and partition coefficients is a fast and reliable alternative. We demonstrate this strategy for transdermal delivery of fentanyl. We found that inverse modeling reduced the normalized root mean square deviation of the measured drug uptake flux from 26 to 9%, when compared to the experimental measurement of all skin properties. We found that this improved agreement with experiments was only possible if the diffusion in the reservoir holding the drug was smaller than the experimentally measured diffusion coefficients suggested. For indirect inverse modeling, which systematically explores the entire parametric space, 30,000 simulations were required. By relying on direct inverse modeling, we reduced the number of simulations to be performed to only 300, so a factor 100 difference. The modeling approach's added value is that it can be calibrated once in-silico for all model parameters simultaneously by solely relying on a single measurement of the drug uptake flux evolution over time. We showed that this calibrated model could accurately be used to simulate transdermal patches with other drug doses. We showed that inverse modeling is a fast way to build up an accurate mechanistic model for drug delivery. This strategy opens the door to clinically ready therapy that is tailored to patients.

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.

The unreliability of the reliability criteria in the estimation of QSAR for skin sensitivity: A pun or a reliable law?

Some new products, which include common personal-care products, drugs, household items, can be hazardous in aspect personal care products/cosmetics and their ingredients (i.e. the above can effect human skin). International organizations (e.g. the Organisation for Economic Co-operation and Development-OECD) recommend evaluating individual ingredients when assessing the safety of personal care or cosmetic products. Thus, checking up that "popular at the market" substances are non-toxic, do not penetrate into or through normal or compromised human skin, and therefore, pose no risk to human health is an essential element of modern toxicology. The development of reliable models of toxicological endpoints is a tool to carry out the above checking up via quantitative structure-activity relationships (QSARs). The reliability of the QSAR is the current task of mathematical statistics. Recently, the index of ideality of correlation (IIC) and correlation intensity index (CII) were suggested as criteria of predictive potential (i.e. reliability) of QSAR-models. Here, the abilities of these criteria were studied for the case of building up models for skin sensitivity (LLNA, local lymph node assay). Computational experiments have confirmed that the IIC demonstrates an obvious ability to improve the predictive potential of models of skin sensitization. The applying of the CII for the case of skin sensitization also improves the quality of the model. However, the best models for skin sensitization were observed if the above-mentioned criteria are applied jointly (n = 268; R² = 0.60; RMSE = 0.63).

Monte-Carlo method-based QSAR model to discover phytochemical urease inhibitors using SMILES and GRAPH descriptors

Urease inhibitors are known to play a vital role in the field of medicine as well as agriculture. Special attention is attributed to the development of novel urease inhibitors with a view to treat the Helicobacter pylori infection. Amongst a number of urease inhibitors, a large number of molecules fail in vivo and in clinical trials due to their hydrolytic instability and toxicity profile. The search for potential inhibitors may require screening of large and diverse databases of small molecules and to design novel molecules. We developed a Monte-Carlo method-based QSAR model to predict urease inhibiting potency of molecules using SMILES and GRAPH descriptors on an existing diverse database of urease inhibitors. The QSAR model satisfies all the statistical parameters required for acceptance as a good model. The model is applied to identify urease inhibitors among the wide range of compounds in the phytochemical database, NPACT, as a test case. We combine the ligand-based and structure-based drug discovery methods to improve the accuracy of the prediction. The method predicts pIC₅₀ and estimates docking score of compounds in the database. The method may be applied to any other database or compounds designed in silico to discover novel drugs targeting urease.Communicated by Ramaswamy H. Sarma.

Development of potential therapeutics for pain treatment by inducing Sigma 1 receptor antagonism – in silico approach

QSAR modeling with computer-aided drug design were used for the in silico development of novel therapeutics for pain treatment.

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.

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

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

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.

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

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

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.