In silico prediction of mosquito repellents for clothing application.

Use of protective clothing is a simple and efficient way to reduce the contacts with mosquitoes and consequently the probability of transmission of diseases spread by them. This mechanical barrier can be enhanced by the application of repellents. Unfortunately the number of available repellents is limited. As a result, there is a crucial need to find new active and safer molecules repelling mosquitoes. In this context, a structure-activity relationship (SAR) model was proposed for the design of repellents active on clothing. It was computed from a dataset of 2027 chemicals for which repellent activity on clothing was measured against Aedes aegypti. Molecules were described by means of 20 molecular descriptors encoding physicochemical properties, topological information and structural features. A three-layer perceptron was used as statistical tool. An accuracy of 87% was obtained for both the training and test sets. Most of the wrong predictions can be explained. Avenues for increasing the performances of the model have been proposed.

Topological QSAR Modelling of Carboxamides Repellent Activity to Aedes Aegypti.

Aedes aegypti vector control is of paramount importance owing to the damages induced by the various severe diseases that these insects may transmit, and the increasing risk of important outbreaks of these pathologies. Search for new chemicals efficient against Aedes aegypti, and devoid of side-effects, which have been associated to the currently most used active substance i. e. N,N-diethyl-m-toluamide (DEET), is therefore an important issue. In this context, we developed various Quantitative Structure Activity Relationship (QSAR) models to predict the repellent activity against Aedes aegypti of 43 carboxamides, by using Multiple Linear Regression (MLR) and various machine learning approaches. The easy computation of the four topological descriptors selected in this study, compared to the CODESSA descriptors used in the literature, and the predictive ability of the here proposed MLR and machine learning models developed using the software QSARINS and R, make the here proposed QSARs attractive. As demonstrated in this study, these models can be applied at the screening level, to guide the design of new alternatives to DEET.

QSAR models for predicting the toxicity of piperidine derivatives against Aedes aegypti.

QSAR models are proposed for predicting the toxicity of 33 piperidine derivatives against Aedes aegypti. From 2D topological descriptors, calculated with the PaDEL software, ordinary least squares multilinear regression (OLS-MLR) treatment from the QSARINS software and machine learning and related approaches including linear and radial support vector machine (SVM), projection pursuit regression (PPR), radial basis function neural network (RBFNN), general regression neural network (GRNN) and k-nearest neighbours (k-NN), led to four-variable models. Their robustness and predictive ability were evaluated through both internal and external validation. Determination coefficients (r2) greater than 0.85 on the training sets and 0.8 on the test sets were obtained with OLS-MLR and linear SVM. They slightly outperform PPR, radial SVM and RBFNN, whereas GRNN and k-NN showed lower performance. The easy availability of the involved structural descriptors and the simplicity of the MLR model make the corresponding model attractive at an exploratory level for proposing, from this limited dataset, guidelines in the design of new potentially active molecules.

Investigation of the influence of protein corona composition on gold nanoparticle bioactivity using machine learning approaches.

The understanding of the mechanisms and interactions that occur when nanomaterials enter biological systems is important to improve their future use. The adsorption of proteins from biological fluids in a physiological environment to form a corona on the surface of nanoparticles represents a key step that influences nanoparticle behaviour. In this study, the quantitative description of the composition of the protein corona was used to study the effect on cell association induced by 84 surface-modified gold nanoparticles of different sizes. Quantitative relationships between the protein corona and the activity of the gold nanoparticles were modelled by using several machine learning-based linear and non-linear approaches. Models based on a selection of only six serum proteins had robust and predictive results. The Projection Pursuit Regression method had the best performances (r(2) = 0.91; Q(2)loo = 0.81; r(2)ext = 0.79). The present study confirmed the utility of protein corona composition to predict the bioactivity of gold nanoparticles and identified the main proteins that act as promoters or inhibitors of cell association. In addition, the comparison of several techniques showed which strategies offer the best results in prediction and could be used to support new toxicological studies on gold-based nanomaterials.

Linear and non-linear modelling of the cytotoxicity of TiO2 and ZnO nanoparticles by empirical descriptors.

Titanium oxide (TiO2) and zinc oxide (ZnO) nanoparticles are among the most widely used in different applications in daily life. In this study, local regression and classification models were developed for a set of ZnO and TiO2 nanoparticles tested at different concentrations for their ability to disrupt the lipid membrane in cells. Different regression techniques were applied and compared by checking the robustness of the models and their external predictive ability. Additionally, a simple classification model was developed, which predicts the potential for disruption of the studied nanoparticles with good accuracy (overall accuracy, specificity, and sensitivity >80%) on the basis of two empirical descriptors. The present study demonstrates that empirical descriptors, such as experimentally determined size and tested concentrations, are relevant to modelling the activity of nanoparticles. This information may be useful to screen the potential for harmful effect of nanoparticles in different experimental conditions and to optimize the design of toxicological tests. Results from the present study are useful to support and refine the future application of in silico tools to nanoparticles, for research and regulatory purposes.

Structure-activity relationship (SAR) modelling of mosquito larvicides.

An attempt was made to derive structure-activity models allowing the prediction of the larvicidal activity of structurally diverse chemicals against mosquitoes. A database of 188 chemicals with their activity on Aedes aegypti larvae was constituted from analysis of original publications. The activity values were expressed in log 1/IC50 (concentration required to produce 50% inhibition of larval development, mmol). All the chemicals were encoded by means of CODESSA and autocorrelation descriptors. Partial least squares analysis, classification and regression tree, random forest and boosting regression tree analyses, Kohonen self-organizing maps, linear artificial neural networks, three-layer perceptrons, radial basis function artificial neural networks and support vector machines with linear, polynomial, radial basis function and sigmoid kernels were tested as statistical tools. Because quantitative models did not give good results, a two-class model was designed. The three-layer perceptron significantly outperformed the other statistical approaches regardless of the threshold value used to split the data into active and inactive compounds. The most interesting configuration included eight autocorrelation descriptors as input neurons and four neurons in the hidden layer. This led to more than 96% of good predictions on both the training set and external test set of 88 and 100 chemicals, respectively. From the overall simulation results, new candidate molecules were proposed which will be shortly synthesized and tested.

Structure-activity relationship study of trifluoromethylketone inhibitors of insect juvenile hormone esterase: comparison of several classification methods.

Juvenile hormone esterase (JHE) plays a key role in the development and metamorphosis of holometabolous insects. Its inhibitors could possibly be targeted for insect control. Conversely, JHE may also be involved in endocrine disruption by xenobiotics, resulting in detrimental effects in beneficial insects. There is therefore a need to know the structural characteristics of the molecules able to monitor JHE activity, and to develop SAR and QSAR studies to estimate their effectiveness. For a large diverse population of 181 trifluoromethylketones (TFKs) - the most potent JHE inhibitors known to date - we recently proposed a binary classification (active/inactive) using a support vector machine and Codessa structural descriptors. We have now examined, using the same data set and with the same descriptors, the applicability and performance of five other machine learning approaches. These have been shown able to handle high dimensional data (with descriptors possibly irrelevant or redundant) and to cope with complex mechanisms, but without delivering explicit directly exploitable models. Splitting the data into five batches (training set 80%, test set 20%) and carrying out leave-one-out cross-validation, led to good results of comparable performance, consistent with our previous support vector classifier (SVC) results. Accuracy was greater than 0.80 for all approaches. A reduced set of 15 descriptors common to all the investigated approaches showed good predictive ability (confirmed using a three-layer perceptron) and gives some clues regarding a mechanistic interpretation.

Structure-activity relationship study of trifluoromethylketones: inhibitors of insect juvenile hormone esterase.

The juvenile hormone esterase (JHE) regulates juvenile hormone titre in insect hemolymph during its larval development. It has been suggested that JHE could be targeted for use in insect control. This enzyme can also be considered as involved in the phenomenon of endocrine disruption by xenobiotics in beneficial insects. Consequently, there is a need to know the characteristics of the molecules able to act on the JHE. Trifluoromethylketones (TFKs) are the most potent JHE inhibitors found to date and different quantitative structure-activity relationships (QSARs) have been derived for this group of chemicals. In this context, a set of 181 TFKs (118 active and 63 inactive compounds), tested on Trichoplusia ni for their JHE inhibition activity and described by physico-chemical descriptors, was split into different training and test sets to derive structure-activity relationship (SAR) models from support vector classification (SVC). A SVC model including 88 descriptors and derived from a Gaussian kernel was selected for its predictive performances. Another model computed only with 13 descriptors was also selected due to its mechanistic interpretability. This study clearly illustrates the difficulty in capturing the essential structural characteristics of the TFKs explaining their JHE inhibitory activity.

Linear and non-linear QSAR modelling of juvenile hormone esterase inhibitors.

A tight control of juvenile hormone (JH) titre is crucial during the life cycle of a holometabolous insect. JH metabolism is made through the action of enzymes, particularly the juvenile hormone esterase (JHE). Trifluoromethylketones (TFKs) are able to inhibit this enzyme to disrupt the endocrine function of the targeted insect. In this context, a set of 96 TFKs, tested on Trichoplusia ni for their JHE inhibition, was split into a training set (n = 77) and a test set (n = 19) to derive a QSAR model. TFKs were initially described by 42 CODESSA (Comprehensive Descriptors for Structural and Statistical Analysis) descriptors, but a feature selection process allowed us to consider only five descriptors encoding the structural characteristics of the TFKs and their reactivity. A classical and spline regression analysis, a three-layer perceptron, a radial basis function network and a support vector regression were experienced as statistical tools. The best results were obtained with the support vector regression (r(2) and r(test)(2) = 0.91). The model provides information on the structural features and properties responsible for the high JHE inhibition activity of TFKs.

Decision trees versus support vector machine for classification of androgen receptor ligands.

With the current concern of limiting experimental assays, increased interest now focuses on in silico models able to predict toxicity of chemicals. Endocrine disruptors cover a large number of environmental and industrial chemicals which may affect the functions of natural hormones in humans and wildlife. Structure-activity models are now increasingly used for predicting the endocrine disruption potential of chemicals. In this study, a large set of about 200 chemicals covering a broad range of structural classes was considered in order to categorize their relative binding affinity (RBA) to the androgen receptor. Classification of chemicals into four activity groups, with respect to their log RBA value, was carried out in a cascade of recursive partitioning trees, with descriptors calculated from CODESSA software and encoding topological, geometrical and quantum chemical properties. The hydrophobicity parameter (log P), Balaban index, and descriptors relying on charge distribution (maximum partial charge, nucleophilic index on oxygen atoms, charged surface area, etc.) appear to play a major role in the chemical partitioning. Separation of strongly active compounds was rather straightforward. Similarly, about 90% of the inactive compounds were identified. More intricate was the separation of active compounds into subsets of moderate and weak binders, the task requiring a more complex tree. A comparison was made with support vector machine yielding similar results.

MolDIA: XML based system of molecular diversity analysis towards virtual screening and QSPR.

In this paper a new chemoinformatics tool for Molecular Diversity Analysis (MolDIA) is introduced. The objective of this system is the analysis of molecular similarity and diversity through the treatment of structural and physicochemical information. Current needs for chemical databases include the analysis, the management and the retrieval of chemical information. The implementation of eXtended Markup Languages (XML) is proposed as a basis for representing and structuring the chemical information contained in data structures and databases. The adequate descriptor vector and related physicochemical properties have been defined and constructed. The benefits of XML in chemoinformatics are discussed, as well as, the applications of this system in a virtual screening environment.

Quantitative structure-toxicity relationships (QSTRs): a comparative study of various non linear methods. General regression neural network, radial basis function neural network and support vector machine in predicting toxicity of nitro- and cyano- aromatics to Tetrahymena pyriformis.

Prediction of toxicity of 203 nitro- and cyano-aromatic chemicals to Tetrahymena pyriformis was carried out by radial basis function neural network, general regression neural network and support vector machine, in non-linear response surface methodology. Toxicity was predicted from hydrophobicity parameter (log Kow) and maximum superdelocalizability (Amax). Special attention was drawn to prediction ability and robustness of the models, investigated both in a leave-one-out and 10-fold cross validation (CV) processes. The influence that the corresponding changes in the learning sets during these CV processes could have on a common external test set including 41 compounds was also examined. This allowed us to establish the stability of the models. The non linear results slightly outperform (as expected) multilinear relationships (MLR) and also favourably compete with various other non linear approaches recently proposed by Ren (J. Chem. Inf. Comput. Sci., 43 1679 (2003)).

QSAR and classification study of 1,4-dihydropyridine calcium channel antagonists based on least squares support vector machines.

The least squares support vector machine (LSSVM), as a novel machine learning algorithm, was used to develop quantitative and classification models as a potential screening mechanism for a novel series of 1,4-dihydropyridine calcium channel antagonists for the first time. Each compound was represented by calculated structural descriptors that encode constitutional, topological, geometrical, electrostatic, quantum-chemical features. The heuristic method was then used to search the descriptor space and select the descriptors responsible for activity. Quantitative modeling results in a nonlinear, seven-descriptor model based on LSSVM with mean-square errors 0.2593, a predicted correlation coefficient (R(2)) 0.8696, and a cross-validated correlation coefficient (R(cv)(2)) 0.8167. The best classification results are found using LSSVM: the percentage (%) of correct prediction based on leave one out cross-validation was 91.1%. This paper provides a new and effective method for drug design and screening.

CISOC-PSCT: a predictive system for carcinogenic toxicity.

A SAR based carcinogenic toxicity prediction system, CISOC-PSCT, was developed. It consisted of two principal phases: the construction of relationships between structural descriptors and carcinogenic toxicity indices, and prediction of the toxicity from the SAR model. The training set included 2738 carcinogenic and 4130 non-carcinogenic compounds. Three predefined topological types of substructures termed Star, Path and Ring were used to generate the descriptors for each structure in the training set. In this system, the defined carcinogenic toxicity index (CTI) was obtained from the probability of a structural descriptor to either belong to the carcinogenic or non-carcinogenic compounds. Based on these structural descriptors and their CTI, a SAR model was derived. Then the carcinogenic possibility (CP) and the carcinogenic impossibility (CIP) of compounds were predicted. The model was tested from a testing set of 304 carcinogenic compounds (MDL toxicity database), 460 non-carcinogenic compounds (CMC database) and 94 compounds extracted from two traditional Chinese medicine herbs.

Comparative study of QSAR/QSPR correlations using support vector machines, radial basis function neural networks, and multiple linear regression.

Support vector machines (SVMs) were used to develop QSAR models that correlate molecular structures to their toxicity and bioactivities. The performance and predictive ability of SVM are investigated and compared with other methods such as multiple linear regression and radial basis function neural network methods. In the present study, two different data sets were evaluated. The first one involves an application of SVM to the development of a QSAR model for the prediction of toxicities of 153 phenols, and the second investigation deals with the QSAR model between the structures and the activities of a set of 85 cyclooxygenase 2 (COX-2) inhibitors. For each application, the molecular structures were described using either the physicochemical parameters or molecular descriptors. In both studied cases, the predictive ability of the SVM model is comparable or superior to those obtained by MLR and RBFNN. The results indicate that SVM can be used as an alternative powerful modeling tool for QSAR studies.

CoMFA/CoMSIA/HQSAR and Docking Study of the Binding Mode of Selective Cyclooxygenase (COX-2) Inhibitors.

The intermolecular interaction between four types of anti-inflammatory inhibitors (oxazoles, pyrazoles, pyrroles and imidazoles) and COX-2 receptor was studied. The results of docking suggest that they have similar interaction mechanism. The most active compounds of these four types of inhibitors could both form several hydrogen bonds with residues His90, Arg513, Leu352 and Arg120, and develop hydrophobic interaction with residues Phe518, Leu352 and Leu359. This is consistent with the investigation reported by R. G. Kurumbail et al. (Nature. 1996, 384, 644-648). A common 3D-QSAR model could be constructed with these four categories of COX-2 inhibitors using the method of docking- guided conformer selection. The cross-validated q2 values are found as 0.741 and 0.632 for CoMFA and CoMSIA respectively. And the non-cross-validated r2 values are 0.887 and 0.885. 54 inhibitors constitute the test set used to validate the model. The results show that this model possesses good predictive ability for diverse COX-2 inhibitors. Furthermore, a HQSAR model was used to evaluate the influence of substituents on anti-inflammatory activity. Compared with the results of previous works, our model possesses significantly better prediction ability. It could help us to well understand the interaction mechanism between inhibitors and COX-2 receptor, and to make quantitative prediction of their inhibitory activities.

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking.

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.

Comparative study of non nucleoside inhibitors with HIV-1 reverse transcriptase based on 3D-QSAR and docking.

The intermolecular interaction between two types of non nucleoside reverse transcriptase inhibitors (NNRTIs), HEPT and TIBO, and HIV reverse transcriptase receptor (HIVRT) was investigated. The result of docking study showed that two types of NNRTIs presented similar interaction mechanism with HIVRT. The most active compound of every type of inhibitors could form one hydrogen bond with the residue Lys101 and has hydrophobic interaction with residues Tyr181, Tyr188 and Tyr318, etc. Three 3D-QSAR models including two partial correlation models (one for each family of HEPT and TIBO) and a mixed model gathering two families were constructed. Comparative study of these models indicated that the mixed model offered the strongest prediction ability. For this model, the cross-validated q2 values were 0.720 and 0.675, non-cross-validated r2 values were 0.940 and 0.920 for CoMFA and CoMSIA, respectively. It has been validated by using a test set of 27 inhibitors. Compared with previously reported works, our model showed better prediction ability. It could help us to insight the interaction between NNRTIs and HIVRT, and to design new anti-HIV NNRTIs inhibitors.

Theoretical study of fast repair of DNA damage by cistanoside C and analogs: mechanism and docking.

Experiments show that the natural substances phenylpropanoid glycosides (PPGs) extracted from pelicularis spicata are capable of repairing DNA damaged by oxygen radicals. Based on kinetic measurements and experiments on tumor cells, a theoretical study of the interaction between PPG molecules and isolated DNA bases, as well as a DNA fragment has been performed. An interaction mechanism reported early has been refined. The docking calculations performed using junction minimization of nucleic acids (JUMNA) software showed that the PPG molecules can be docked into the minor groove of DNA and form complexes with the geometry suitable for an electron transfer between guanine radical and the ligand. Such complexes can be formed without major distortions of DNA structure and are further stabilized by the interaction with the rhamnosyl side-groups.

Quantitative prediction of liquid chromatography retention of N-benzylideneanilines based on quantum chemical parameters and radial basis function neural network.

Based on quantum chemical parameters and a simple numerical coding, the liquid chromatography retention of bifunctionally substituted N-benzylideneaniles (NBA) has been predicted using a radial basis function neural network (RBFNN) model. The quantum chemical parameters involved in the model are dipole moment (m), energies of the highest occupied and lowest unoccupied molecular orbitals (E(homo,) E(lumo)), net charge of the most negative atom (Q(min)), sum of absolute values of the charges of all atoms in two given functional groups (Delta), total energy of the molecule (E(T)), weight of the molecule (W), and numerical coding (N). N was used to indicate the different positions of two substituents. The predictive values are consistent with the experimental results. The mean relative error of the testing set is 1.6%, and the maximum relative error is less than 5.0%. In this work the success of the whole modeling process only depends on the optimization of the spread parameter in network.