Replacement Method and Enhanced Replacement Method Versus the Genetic Algorithm Approach for the Selection of Molecular Descriptors in QSPR/QSAR Theories

Development of QSARs for cysteine-containing di- and tripeptides with antioxidant activity:influence of the cysteine position

Antioxidants agents play an essential role in the food industry for improving the oxidative stability of food products. In the last years, the search for new natural antioxidants has increased due to the potential high toxicity of chemical additives. Therefore, the synthesis and evaluation of the antioxidant activity in peptides is a field of current research. In this study, we performed a Quantitative Structure Activity Relationship analysis (QSAR) of cysteine-containing 19 dipeptides and 19 tripeptides. The main objective is to bring information on the relationship between the structure of peptides and their antioxidant activity. For this purpose, 1D and 2D molecular descriptors were calculated using the PaDEL software, which provides information about the structure, shape, size, charge, polarity, solubility and other aspects of the compounds. Different QSAR model for di- and tripeptides were developed. The statistic parameters for di-peptides model (R2train = 0.947 and R2test = 0.804) and for tripeptide models (R2train = 0.923 and R2test = 0.847) indicate that the generated models have high predictive capacity. Then, the influence of the cysteine position was analyzed predicting the antioxidant activity for new di- and tripeptides, and comparing them with glutathione. In dipeptides, excepting SC, TC and VC, the activity increases when cysteine is at the N-terminal position. For tripeptides, we observed a notable increase in activity when cysteine is placed in the N-terminal position.

In Silico Antiprotozoal Evaluation of 1,4-Naphthoquinone Derivatives against Chagas and Leishmaniasis Diseases Using QSAR, Molecular Docking, and ADME Approaches

Chagas and leishmaniasis are two neglected diseases considered as public health problems worldwide, for which there is no effective, low-cost, and low-toxicity treatment for the host. Naphthoquinones are ligands with redox properties involved in oxidative biological processes with a wide variety of activities, including antiparasitic. In this work, in silico methods of quantitative structure–activity relationship (QSAR), molecular docking, and calculation of ADME (absorption, distribution, metabolism, and excretion) properties were used to evaluate naphthoquinone derivatives with unknown antiprotozoal activity. QSAR models were developed for predicting antiparasitic activity against Trypanosoma cruzi, Leishmania amazonensis, and Leishmania infatum, as well as the QSAR model for toxicity activity. Most of the evaluated ligands presented high antiparasitic activity. According to the docking results, the family of triazole derivatives presented the best affinity with the different macromolecular targets. The ADME results showed that most of the evaluated compounds present adequate conditions to be administered orally. Naphthoquinone derivatives show good biological activity results, depending on the substituents attached to the quinone ring, and perhaps the potential to be converted into drugs or starting molecules.

Exploratory and machine learning analysis of the stability constants of HgII- triazene ligands complexes

Knowing stability constants for the complexes HgII with extracting ligands is very important from environmental and therapeutic standpoints. Since the selectivity of ligands can be stated by the stability constants of cation–ligand complexes, quantitative structure–property relationship (QSPR) investigations on binding constant of HgII complexes were done. Experimental data of the stability constants in ML2 complexation of HgII and synthesized triazene ligands were used to construct and develop QSPR models. Support vector machine (SVM) and multiple linear regression (MLR) have been employed to create the QSPR models. The final model showed squared correlation coefficient of 0.917 and the standard error of calibration (SEC) value of 0.141 log K units. The proposed model presented accurate prediction with the Leave-One-Out cross validation ( Q LOO 2  = 0.756) and validated using Y-randomization and external test set. Statistical results demonstrated that the proposed models had suitable goodness of fit, predictive ability, and robustness. The results revealed the importance of charge effects and topological properties of ligand in HgII - triazene complexation.

Predicting zebrafish (Danio rerio) embryo developmental toxicity through a non-conformational QSAR approach

For many years, the frequent use of synthetic chemicals in the manufacture of veterinary drugs and plague control products has raised negative effects on human health and other non-target organisms, promoting the need to employ a practical and suitable methodology for early risk identification of several thousand commercial compounds. The zebrafish (Danio rerio) embryo has been emerged as one sustainable animal model for measuring developmental toxicity, an endpoint that is included in the regulatory procedures to approve chemicals, avoiding conventional and costly toxicity assays based on animal testing. In this context, the Quantitative Structure-Activity Relationships (QSAR) theory is applied to develop a predictive model based on a well-defined zebrafish embryo developmental toxicity database reported by the ToxCast™ Phase I chemical library of the Environmental Protection Agency (U.S. EPA). By means of four freely available softwares, a set with 28,038 non-conformational descriptors that encode the largest amount of permanent structural features are readily calculated. The Replacement Method (RM) variable subset selection technique provided the best regression models. Thereby, a linear QSAR model with proper statistical quality (R_train² = 0.64, RMSE_train = 0.49) is established in agreement with the Organization for Economic Co-operation and Development principles, accomplishing each internal (loo, l15 % o, VIF and Y-randomization) and external (R_test²,R_m², Q_F1², Q_F2², Q_F3² and CCC) validation criterion. The present QSAR approach provides a useful computational tool to estimate zebrafish developmental toxicity of new, untasted or hypothetical compounds, and it can contribute to the general lack of QSAR models in the literature to predict this endpoint.

Antiprotozoal QSAR modelling for trypanosomiasis (Chagas disease) based on thiosemicarbazone and thiazole derivatives

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, remains a neglected endemic infection that affects around 8 million people worldwide and causes 12,000 premature deaths per year. Traditional chemotherapy is limited to the nitro-antiparasitic drugs Benznidazole and Nifurtimox, which present serious side effects and low long-term efficacy. Several research efforts have been made over the last decade to find new chemical structures with better effectiveness and tolerance than standard anti-Chagas drugs. Among these, new sets of thiosemicarbazone and thiazole derivatives have exhibited potent in vitro activity against T. cruzi, especially for its extracellular forms (epimastigote and trypomastigote). In this work, we have developed three antiprotozoal quantitative structure-relationship (QSAR) models for Chagas disease based on the in vitro activity data reported as IC₅₀ (μM) and CC₅₀ (μM) over the last decade, particularly by Lima-Leite's group in Brazil. The models were developed using the replacement method (RM), a technique based on Multivariable Linear Regression (MLR), and external and internal validation methodologies, like the use of a test set, Leave-one-Out (LOO) cross-validation and Y-Randomization. Two of these QSAR models were developed for trypomastigotes form of the parasite Trypanosoma cruzi, one based on IC₅₀ and the other on CC₅₀ data; while the third QSAR model was developed for its epimastigotes form based on CC₅₀ activity. Our models presented sound statistical parameters that endorses their prediction capability. Such capability was tested for a set of 13 hitherto-unknown structurally related aromatic cyclohexanone derivatives.

Prediction of toxicity of Ionic Liquids based on GC-COSMO method

In order to evaluate the toxicity of several different ionic liquids (ILs) towards the leukemia rat cell line (ICP-81), an efficient and reliable quantitative structure-activity relationships (QSAR) model is developed based on descriptors from COSMO-SAC (conductor-like screening model for segment activity coefficient) model. The distribution of screen charge density (σ-profile) of 127 ILs is calculated by GC-COSMO (group contribution based COSMO) method. Two segmentation methods toward σ-profile are used to find out the appropriate descriptors for the QSAR model. The optimal subset of descriptors is obtained by enhanced replacement method (ERM). A multiple linear regression (MLR) and multilayer perceptron technique (MLP) are used to build the linear and nonlinear models, respectively, and the applicability domain of the models is assessed by the Williams plot. It turns out that the nonlinear model based the second segmentation method (MLP-2) is the best QSAR model with an R²=0.975, MSE=0.019 for the training set and R²=0.938, MSE=0.037 for the test set. The reliability and robustness of the presented QSAR models are confirmed by Leave-One-Out (LOO) cross and external validations.

Predicting pungency and understanding the pungency mechanism of capsaicinoids using TOPS-MODE approach

Quantitative structure-property relationship (QSPR) models were developed for predicting the pungency of a set of capsaicinoids. Multiple linear regression (MLR) coupled with topological substructural molecular descriptor (TOPS-MODE) approach was used. The best MLR model based on only five orthogonalized TOPS-MODE variables allowed us to obtain a coefficient of determination of 0.954 on the training set. The predictive power of the model was validated through a test set and several external validation parameters. This showed that the TOPS-MODE descriptors weighted by bond dipole moments, van der Waals atomic radii, and the total solute hydrogen bond basicity affected pungency. The contributions of certain bonds and fragments to pungency were used to understand the pungency mechanism of capsaicinoids. The selected model can more accurately predict pungency of capsaicinoids compared than those found in the literature, and especially bring insights into the structural features and chemical factors related to pungency.

QSAR studies of the antioxidant activity of anthocyanins

Through experimental information available from antioxidant assays of seventeen anthocyanins, and six common anthocyanidins, quantitative structure-activity relationships (QSAR) have been established in the present work. The antioxidant bioactivity has been predicted in three different lipid environments: emulsified and bulk oil (methyl linoleate) (in vitro tests) at concentrations of 50 and 250 μM, and 50 μM of the inhibitor, respectively, and in human LDL (low-density lipoprotein; "bad cholesterol") (ex vivo test) at concentrations of 2.5, 10, and 25 μM of the inhibitor. Radical scavenging activity was predicted in the assay with the 1,1-diphenyl-2-picrylhydrazyl radical (DPPH·). The QSAR models developed for each test and concentration used allowed to obtain prospective information on the constitutional and topological molecular characteristics for anthocyanin/anthocyanidin compounds. Therefore, the antioxidant activity was predicted for twenty-one compounds with unknown experimental values, leading for some of them to a favorable predicted bioactivity.

An exploratory study using QICAR models for prediction of adsorption capacity of multi-walled carbon nanotubes for heavy metal ions

The Quantitative Ion Character-Activity Relationship (QICAR) method was used for correlating metal ionic characteristics with the maximum adsorption capacity (q_max) of multi-walled carbon for heavy metals. The experimental values of q_max for 25 heavy metal ions, estimated by the Langmuir isotherm model, were used to construct a QICAR model. The genetic algorithm, enhanced replacement method and successive projection algorithm procedures were applied as variable selection algorithms to choose the optimal subsets of descriptors. The selected variables were correlated with q_max values by using partial least squares (PLS) regression. Orthogonal signal correction was applied as a pre-processing technique. Among of different variable selection methods, the enhanced replacement method displayed noticeable statistical parameters of the final model. The results of the enhancement replacement method-orthogonal correction signal-PLS model, with RMSEC = 0.733, r²_c = 0.999 and r²_p = 0.946, were excellent and dramatically better than those of other models. The developed QICAR model satisfied the internal and external validation criteria. The importance of electronegativity, ionic radius and atomic number of the heavy metal ions indicated the impact of the tendency to accept electrons and the size of ions in adsorption on carbon nanotubes.

Identification of new potential HIV-1 reverse transcriptase inhibitors by QSAR modeling and structure-based virtual screening

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) have gained a definitive place due to their unique antiviral potency, high specificity and low toxicity in antiretroviral combination therapies which are used to treat HIV. To design more specific HIV-1 inhibitors, 218 diverse non-nucleoside reverse transcriptase inhibitors with their EC₅₀ values were collected. Then, different types of molecular descriptors were calculated. Also, genetic algorithm (GA) and enhanced replacement methods (ERM) were used as the variable selection approaches to choose more relevant features. Based on selected descriptors, a classification support vector machine (SVM) model was constructed to categorize compounds into two groups of active and inactive ones. The most active compound in the set was docked and was used as the input to the Pharmit server to screen the Molport and PubChem libraries by constructing a structure-based pharmacophore model. Shape filters for the protein and ligand as well as Lipinski's rule of five have been applied to filter out the output of virtual screening from pharmacophore search. Three hundred and thirty-four compounds were finally retrieved from the virtual screening and were fed to the previously constructed SVM model. Among them, the SVM model rendered seven active compounds and they were also analyzed by docking calculations and ADME/Tox parameters.

Impact assessment of the rational selection of training and test sets on the predictive ability of QSAR models

This study performed an analysis of the influence of the training and test set rational selection on the quality and predictively of the quantitative structure-activity relationship (QSAR) model. The study was carried out on three different datasets of Influenza Neuraminidase (H1N1) inhibitors. The three datasets were divided into training and test sets using three rational selection methods: based on k-means, Kennard-Stone algorithm and Activity and the results were compared with Random selection. Then, a total of 31,490 mathematical models were developed and those models that presented a determination coefficient higher than: r²_train > 0.8, r²_loo > 0.7, r²_test > 0.5 and minimum standard deviation (SD) and minimum root-mean square error (RMS) were selected. The selected models were validated using the internal leave-one-out method and the predictive capacity was evaluated by the external test set. The results indicate that random selection could lead to erroneous results. In return, a rational selection allows for obtaining more reliable conclusions. The QSAR models with major predictive power were found using the k-means algorithm and selection by activity.

What do docking and QSAR tell us about the design of HIV-1 reverse transcriptase nonnucleoside inhibitors?

Despite vigorous studies, effective nonnucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) are still in demand, not only due to toxicity and detrimental side effects of currently used drugs but also because of the emergence of multidrug-resistant viral strains. In this contribution, we present results of docking of 47 inhibitors to 107 allosteric centers of HIV-1 reverse transcriptase. Based on the average binding scores, we have constructed QSAR equations to elucidate directions of further developments in the inhibitor design that come from this structural data.

An Introduction to the Basic Concepts in QSAR-Aided Drug Design

The need for the development of new drugs to combat existing and newly identified conditions is unavoidable. One of the important tools used in the advanced drug development pipeline is computer-aided drug design. Traditionally, to find a drug many ligands were synthesized and evaluated for their effectiveness using suitable bioassays and if all other drug-likeness features were met, the candidate(s) would possibly reach the market. Although this approach is still in use in advanced format, computational methods are an indispensable component of modern drug development projects. One of the methods used from very early days of rationalizing the drug design approaches is Quantitative Structure-Activity Relationship (QSAR). This chapter overviews QSAR modeling steps by introducing molecular descriptors, mathematical model development for relating biological activities to molecular structures, and model validation. At the end, several successful cases where QSAR studies were used extensively are presented.

Linear and nonlinear models for predicting fish bioconcentration factors for pesticides

This work is devoted to the applications of the multiple linear regression (MLR), multilayer perceptron neural network (MLP NN) and projection pursuit regression (PPR) to quantitative structure-property relationship analysis of bioconcentration factors (BCFs) of pesticides tested on Bluegill (Lepomis macrochirus). Molecular descriptors of a total of 107 pesticides were calculated with the DRAGON Software and selected by inverse enhanced replacement method. Based on the selected DRAGON descriptors, a linear model was built by MLR, nonlinear models were developed using MLP NN and PPR. The robustness of the obtained models was assessed by cross-validation and external validation using test set. Outliers were also examined and deleted to improve predictive power. Comparative results revealed that PPR achieved the most accurate predictions. This study offers useful models and information for BCF prediction, risk assessment, and pesticide formulation.

Application of quantitative structure-property relationship analysis to estimate the vapor pressure of pesticides

The application of molecular descriptors in describing Quantitative Structure Property Relationships (QSPR) for the estimation of vapor pressure (VP) of pesticides is of ongoing interest. In this study, QSPR models were developed using multiple linear regression (MLR) methods to predict the vapor pressure values of 162 pesticides. Several feature selection methods, namely the replacement method (RM), genetic algorithms (GA), stepwise regression (SR) and forward selection (FS), were used to select the most relevant molecular descriptors from a pool of variables. The optimum subset of molecular descriptors was used to build a QSPR model to estimate the vapor pressures of the selected pesticides. The Replacement Method improved the predictive ability of vapor pressures and was more reliable for the feature selection of these selected pesticides. The results provided satisfactory MLR models that had a satisfactory predictive ability, and will be important for predicting vapor pressure values for compounds with unknown values. This study may open new opportunities for designing and developing new pesticide.

QSAR models for predicting octanol/water and organic carbon/water partition coefficients of polychlorinated biphenyls

Quantitative structure-property relationship modelling can be a valuable alternative method to replace or reduce experimental testing. In particular, some endpoints such as octanol-water (KOW) and organic carbon-water (KOC) partition coefficients of polychlorinated biphenyls (PCBs) are easier to predict and various models have been already developed. In this paper, two different methods, which are multiple linear regression based on the descriptors generated using Dragon software and hologram quantitative structure-activity relationships, were employed to predict suspended particulate matter (SPM) derived log KOC and generator column, shake flask and slow stirring method derived log KOW values of 209 PCBs. The predictive ability of the derived models was validated using a test set. The performances of all these models were compared with EPI Suite™ software. The results indicated that the proposed models were robust and satisfactory, and could provide feasible and promising tools for the rapid assessment of the SPM derived log KOC and generator column, shake flask and slow stirring method derived log KOW values of PCBs.

Alignment independent 3D-QSAR, quantum calculations and molecular docking of Mer specific tyrosine kinase inhibitors as anticancer drugs

Mer receptor tyrosine kinase is a promising novel cancer therapeutic target in many human cancers, because abnormal activation of Mer has been implicated in survival signaling and chemoresistance. 3D-QSAR analyses based on alignment independent descriptors were performed on a series of 81 Mer specific tyrosine kinase inhibitors. The fractional factorial design (FFD) and the enhanced replacement method (ERM) were applied and tested as variable selection algorithms for the selection of optimal subsets of molecular descriptors from a much greater pool of such regression variables. The data set was split into 65 molecules as the training set and 16 compounds as the test set. All descriptors were generated by using the GRid INdependent descriptors (GRIND) approach. After variable selection, GRIND were correlated with activity values (pIC50) by PLS regression. Of the two applied variable selection methods, ERM had a noticeable improvement on the statistical parameters of PLS model, and yielded a q (2) value of 0.77, an [Formula: see text] of 0.94, and a low RMSEP value of 0.25. The GRIND information contents influencing the affinity on Mer specific tyrosine kinase were also confirmed by docking studies. In a quantum calculation study, the energy difference between HOMO and LUMO (gap) implied the high interaction of the most active molecule in the active site of the protein. In addition, the molecular electrostatic potential energy at DFT level confirmed results obtained from the molecular docking. The identified key features obtained from the molecular modeling, enabled us to design novel kinase inhibitors.

Quantitative Structure-Antifungal Activity Relationships for cinnamate derivatives

Quantitative Structure-Activity Relationships (QSAR) are established with the aim of analyzing the fungicidal activities of a set of 27 active cinnamate derivatives. The exploration of more than a thousand of constitutional, topological, geometrical and electronic molecular descriptors, which are calculated with Dragon software, leads to predictions of the growth inhibition on Pythium sp and Corticium rolfsii fungi species, in close agreement to the experimental values extracted from the literature. A set containing 21 new structurally related cinnamate compounds is prepared. The developed QSAR models are applied to predict the unknown fungicidal activity of this set, showing that cinnamates like 38, 28 and 42 are expected to be highly active for Pythium sp, while this is also predicted for 28 and 34 in C. rolfsii.

Global and local QSPR models to predict supercooled vapour pressure for organic compounds

In this study, a quantitative structure-property relationship (QSPR) approach was used for estimation of logarithmic values of supercooled liquid vapour pressure (log P_L) of a large set of structurally diverse organic compounds. This set includes 12 local sets of aromatic and aliphatic hydrocarbons, polychlorinated biphenyls, ethers, polychlorinated and brominated diphenylethers, polychlorinated naphthalenes and alcohols. Some simple models based on the linear relationship between log P_L and VolSurf descriptors were developed as global models, and a general equation as a simple way to calculate the supercooled liquid vapour pressure of organic chemicals was provided. A descriptor representing the hydrophilic regions (WO₁) of organic chemicals showed the highest correlation with log P_L and resulted in a one-parameter global model characterized by satisfactory statistical performance; calibration (r²_c) and prediction (r²_p) correlation coefficient of 0.84 and 0.85, respectively. Moreover, local QSPR models were also developed for each subset of organic compounds and, as expected, the statistical results obtained from these models were better than the global one. From the descriptors involved in the models, it is concluded that the hydrophilic and hydrophobic regions at different energy levels and polarizability usually determine the variation of supercooled liquid vapour pressure of organic compounds.

An Introduction to the Basic Concepts in QSAR-Aided Drug Design

The need for the development of new drugs to combat existing and newly identified conditions is unavoidable. One of the important tools used in the advanced drug development pipeline is computer-aided drug design. Traditionally, to find a drug many ligands were synthesized and evaluated for their effectiveness using suitable bioassays and if all other drug-likeness features were met, the candidate(s) would possibly reach the market. Although this approach is still in use in advanced format, computational methods are an indispensable component of modern drug development projects. One of the methods used from very early days of rationalizing the drug design approaches is Quantitative Structure-Activity Relationship (QSAR). This chapter overviews QSAR modeling steps by introducing molecular descriptors, mathematical model development for relating biological activities to molecular structures, and model validation. At the end, several successful cases where QSAR studies were used extensively are presented.

QSAR models for thiophene and imidazopyridine derivatives inhibitors of the Polo-Like Kinase 1

The inhibitory activity of 103 thiophene and 33 imidazopyridine derivatives against Polo-Like Kinase 1 (PLK1) expressed as pIC50 (-logIC50) was predicted by QSAR modeling. Multivariate linear regression (MLR) was employed to model the relationship between 0D and 3D molecular descriptors and biological activities of molecules using the replacement method (MR) as variable selection tool. The 136 compounds were separated into several training and test sets. Two splitting approaches, distribution of biological data and structural diversity, and the statistical experimental design procedure D-optimal distance were applied to the dataset. The significance of the training set models was confirmed by statistically higher values of the internal leave one out cross-validated coefficient of determination (Q2) and external predictive coefficient of determination for the test set (Rtest2). The model developed from a training set, obtained with the D-optimal distance protocol and using 3D descriptor space along with activity values, separated chemical features that allowed to distinguish high and low pIC50 values reasonably well. Then, we verified that such model was sufficient to reliably and accurately predict the activity of external diverse structures. The model robustness was properly characterized by means of standard procedures and their applicability domain (AD) was analyzed by leverage method.

QSAR analysis on tacrine-related acetylcholinesterase inhibitors

Background

The evaluation of the clinical effects of Tacrine has shown efficacy in delaying the deterioration of the symptoms of Alzheimer’s disease, while confirming the adverse events consisting mainly in the elevated liver transaminase levels. The study of tacrine analogs presents a continuous interest, and for this reason we establish Quantitative Structure-Activity Relationships on their Acetylcholinesterase inhibitory activity.

Results

Ten groups of new developed Tacrine-related inhibitors are explored, which have been experimentally measured in different biochemical conditions and AChE sources. The number of included descriptors in the structure-activity relationship is characterized by ‘Rule of Thumb’. The 1502 applied molecular descriptors could provide the best linear models for the selected Alzheimer’s data base and the best QSAR model is reported for the considered data sets.

Conclusion

The QSAR models developed in this work have a satisfactory predictive ability, and are obtained by selecting the most representative molecular descriptors of the chemical structure, represented through more than a thousand of constitutional, topological, geometrical, quantum-mechanical and electronic descriptor types.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0084-0) contains supplementary material, which is available to authorized users.

3D-QSAR studies on the toxicity of substituted benzenes to Tetrahymena pyriformis: CoMFA, CoMSIA and VolSurf approaches

Three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis were performed on the toxicity of a large set of substituted benzenes toward ciliate Tetrahymena pyriformis. The 3D-QSAR studies were carried out using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and VolSurf techniques. The optimal CoMFA and CoMSIA models obtained from the training set were all statistically significant with correlation coefficients (R(2)) greater than 0.79 and absolute error less than 0.33 in log units. The predictive ability of the models was externally evaluated through the prediction of a test set (20 percent of the whole data set) that were not included in the training set. A simple and fairly good predictive linear model based on VolSurf descriptors was also developed that showed an adequate prediction power of the toxicity (pIGC50) of substituted benzenes. Validation, reliability and robustness of models were also evaluated by leave-one-out, leave-four-out, bootstrapping and progressive scrambling approaches. The results confirmed that in addition to hydrophobic effects, electrostatic and H-bonding interactions also play important roles in the toxicity of substituted benzenes. The information obtained from CoMFA and CoMSIA 3-D contour maps could be useful to explain the toxicity mechanism of substituted benzenes.

QSAR study of liver specificity of carcinogenicity of N-nitroso compounds

The quantitative structure-activity relationship (QSAR) of N-nitroso compounds (NOCs) for rat liver was developed by a topological sub-structural molecular-descriptors (TOPS-MODE) approach to predict non-liver-carcinogenic and liver-carcinogenic N-nitroso compounds based on a data set of 108 NOCs. Three descriptors calculated solely from the molecular structures of the compounds were selected by enhanced replacement method (ERM) and were weighted, respectively, with atomic weight, bond dipole moments and Abraham solute descriptor partition between water and aqueous solvent systems to indicate the importance of their roles in liver specificity. A detailed discussion on these three descriptors was carried out, and the contributions of different fragments to rat-liver specificity and the interactions among fragments were analyzed. Such results can offer some useful theoretical references for understanding the chemical structural and biological factors related to the liver-specific biological activity of NOCs.

Design of skin penetration enhancers using replacement methods for the selection of the molecular descriptors

Transdermal delivery of certain drugs is challenging because of skin barrier resistance. This study focuses on the implementation of feature-selection algorithms to design chemical penetration enhancers. A database, consisting of 145 polar and nonpolar chemicals, was chosen for the investigation. Replacement, enhanced replacement and stepwise algorithms were applied to identify relevant structural properties of these compounds. The descriptors were calculated using Molecular Modeling Pro™ Plus. Based on the coefficient of determination, the replacement methods outperformed the stepwise approach in selecting the features that best correlated with the flux enhancement ratio. An artificial neural network model was built to map a subset of descriptors from sixty-one nonpolar enhancers onto the output vector. The R² value improved from 0.68, for a linear model, to 0.74, which shows that the improved framework might be effective in the design of compounds with user-defined properties.

Quantitative structure-activity relationships of antimicrobial fatty acids and derivatives against Staphylococcus aureus

Fatty acids and derivatives (FADs) are resources for natural antimicrobials. In order to screen for additional potent antimicrobial agents, the antimicrobial activities of FADs against Staphylococcus aureus were examined using a microplate assay. Monoglycerides of fatty acids were the most potent class of fatty acids, among which monotridecanoin possessed the most potent antimicrobial activity. The conventional quantitative structure-activity relationship (QSAR) and comparative molecular field analysis (CoMFA) were performed to establish two statistically reliable models (conventional QSAR: R(2)=0.942, Q(2)(LOO)=0.910; CoMFA: R(2)=0.979, Q(2)=0.588, respectively). Improved forecasting can be achieved by the combination of these two models that provide a good insight into the structure-activity relationships of the FADs and that may be useful to design new FADs as antimicrobial agents.

Improvement of the Prediction Power of the CoMFA and CoMSIA Models on Histamine H3 Antagonists by Different Variable Selection Methods

The aim of this study is to enhance the predictivity power of CoMFA and CoMSIA models by means of different variable selection algorithms. The genetic algorithm (GA), successive projection algorithm (SPA), stepwise multiple linear regression (SW-MLR), and the enhanced replacement method (ERM) were used and tested as variable selection algorithms. Then, the selected variables were used to generate a simple and predictive model by the multilinear regression algorithm. A set of 74 histamine H₃ antagonists were split into 40 compounds as a training set, and 17 compounds as a test set, by the Kennard-Stone algorithm. Before splitting the data, 17 compounds were randomly selected from the pool of the whole data set as an evaluation set without any supervision, pretreatment, or visual inspection. Among applied variable selection algorithms, ERM had noticeable improvement on the statistical parameters. The r² values of training, test, and evaluation sets for the ERM-MLR model using CoMFA fields were 0.9560, 0.8630, and 0.8460 and using the CoMSIA fields were 0.9800, 0.8521, and 0.9080, respectively. In this study, the principles of organization for economic cooperation and development (OECD) for regulatory acceptability of QSARs are considered.

Quantitative structure-activity relationships of mosquito larvicidal chalcone derivatives

The mosquito larvicidal activities of a series of chalcones and some derivatives were subjected to a quantitative structure-activity relationship (QSAR) study, using more than a thousand constitutional, topological, geometrical, and electronic molecular descriptors calculated with Dragon software. The larvicidal activity values for 28 active compounds of the series were predicted, showing in general a good approximation to the experimental values found in the literature. Chalcones having one or both electron-rich rings showed high toxicity. However, the activity of chalcones was reduced by electron-withdrawing groups, and this was roughly diminished by derivatization of the carbonyl group. A set of six chalcones being structurally similar to some of the active ones, with a still unknown larvicidal activity, were prepared. Their activity values were predicted by applying the developed QSAR models, showing that two chalcones of such set, both 32 and 34, were expected to be highly active.

Predicting carcinogenicity and understanding the carcinogenic mechanism of N-nitroso compounds using a TOPS-MODE approach

A linear discriminant analysis (LDA) coupled with an enhanced replacement method (ERM) was used as an alternative method to predict the carcinogenicity of N-nitroso compounds (NOCs) in rats. This presented LDA based on the topological substructural molecular descriptors (TOPS-MODE) approach was developed to predict the carcinogenic and noncarcinogenic activity on a data set of 111 NOCs with a good classification value of 90.1%. The predictive power of the LDA model was validated through an external validation set (37 compounds) with a prediction accuracy of 94.6% and a leave-one-out cross-validation procedure (LOOCV) with a good prediction of 86.5%. This methodology showed that the TOPS-MODE descriptors weighted, respectively, by bond dipole moment and Abraham solute descriptor dipolarity/polarizability affected the NOC carcinogenicity. The contributions of certain bonds and fragments to carcinogenicity were used to assess biotransformation and carcinogenic mechanisms. The positive contribution of the carbon-nitrogen single bond (between the N-nitroso group and α-carbon to the N-nitroso group) indicated that the α-hydroxylation reaction could occur at the α-carbon or otherwise not occur. Similarly, the contributions from the molecular fragment could be applied to indicate whether the fragments generated an alkylating agent. These results suggested that this approach could discriminate between carcinogenic and noncarcinogenic NOCs, thereby providing insight into the structural features and chemical factors related to NOC carcinogenicity.