Ligand and structure-based discovery of phosphorus-containing compounds as potential metalloproteinase inhibitors

In this study, a methodology is proposed, combining ligand- and structure-based virtual screening tools, for the identification of phosphorus-containing compounds as inhibitors of zinc metalloproteases. First, we use Dragon molecular descriptors to develop a Linear Discriminant Analysis classification model, which is widely validated according to the OECD principles. This model is simple, robust, stable and has good discriminating power. Furthermore, it has a defined applicability domain and it is used for virtual screening of the DrugBank database. Second, docking experiments are carried out on the identified compounds that showed good binding energies to the enzyme thermolysin. Considering the potential toxicity of phosphorus-containing compounds, their toxicological profile is evaluated according to Protox II. Of the five molecules evaluated, two show carcinogenic and mutagenic potential at small LD50, not recommended as drugs, while three of them are classified as non-toxic, and could constitute a starting point for the development of new vasoactive metalloprotease inhibitor drugs. According to molecular dynamics simulation, two of them show stable interactions with the active site maintaining coordination with the metal. A high agreement is evident between QSAR, docking and molecular dynamics results, demonstrating the potentialities of the combination of these tools.

Exploring proteasome inhibition using atomic weighted vector indices and machine learning approaches

Ubiquitin-proteasome system (UPS) is a highly regulated mechanism of intracellular protein degradation and turnover. The UPS is involved in different biological activities, such as the regulation of gene transcription and cell cycle. Several researchers have applied cheminformatics and artificial intelligence methods to study the inhibition of proteasomes, including the prediction of UPP inhibitors. Following this idea, we applied a new tool for obtaining molecular descriptors (MDs) for modeling proteasome Inhibition in terms of EC₅₀ (µmol/L), in which a set of new MDs called atomic weighted vectors (AWV) and several prediction algorithms were used in cheminformatics studies. In the manuscript, a set of descriptors based on AWV are presented as datasets for training different machine learning techniques, such as linear regression, multiple linear regression (MLR), random forest (RF), K-nearest neighbors (IBK), multi-layer perceptron, best-first search, and genetic algorithm. The results suggest that these atomic descriptors allow adequate modeling of proteasome inhibitors despite artificial intelligence techniques, as a variant to build efficient models for the prediction of inhibitory activity.

A Review of Computational Approaches Targeting SARS-CoV-2 Main Protease to the Discovery of New Potential Antiviral Compounds

The new pandemic caused by the coronavirus (SARS-CoV-2) has become the biggest challenge that the world is facing today. It has been creating a devastating global crisis, causing countless deaths and great panic. The search for an effective treatment remains a global challenge owing to controversies related to available vaccines. A great research effort (clinical, experimental, and computational) has emerged in response to this pandemic, and more than 125000 research reports have been published in relation to COVID-19. The majority of them focused on the discovery of novel drug candidates or repurposing of existing drugs through computational approaches that significantly speed up drug discovery. Among the different used targets, the SARS-CoV-2 main protease (M^pro), which plays an essential role in coronavirus replication, has become the preferred target for computational studies. In this review, we examine a representative set of computational studies that use the M^pro as a target for the discovery of small-molecule inhibitors of COVID-19. They will be divided into two main groups, structure-based and ligand-based methods, and each one will be subdivided according to the strategies used in the research. From our point of view, the use of combined strategies could enhance the possibilities of success in the future, permitting to development of more rigorous computational studies in future efforts to combat current and future pandemics.

Anti-inflammatory Activity: In silico and In vivo of Sapogenins Present in Agave brittoniana subsp. brachypus (Trel.)

BACKGROUND

Agave brittoniana subsp. brachypus is an endemic plant of Cuba, which contains different steroidal sapogenins with anti-inflammatory effects. This work aims to develop computational models which allow the identification of new chemical compounds with potential anti-inflammatory activity.

METHODS

The in vivo anti-inflammatory activity was evaluated in two rat models: carrageenaninduced paw edema and cotton pellet-induced granuloma. In each study, we used 30 Sprague Dawley male rats divided into five groups containing six animals. The products isolated and administrated were fraction rich in yuccagenin and sapogenins crude.

RESULTS

The obtained model, based on a classification tree, showed an accuracy value of 86.97% for the training set. Seven compounds (saponins and sapogenins) were identified as potential antiinflammatory agents in the virtual screening. According to in vivo studies, the yuccagenin-rich fraction was the greater inhibitor of the evaluated product from Agave.

CONCLUSION

The evaluated metabolites of the Agave brittoniana subsp. Brachypus showed an interesting anti-inflammatory effect.

Ligand-based discovery of new potential acetylcholinesterase inhibitors for Alzheimer's disease treatment

The enzyme acetylcholinesterase (AChE) is currently a therapeutic target for the treatment of neurodegenerative diseases. These diseases have highly variable causes but irreversible evolutions. Although the treatments are palliative, they help relieve symptoms and allow a better quality of life, so the search for new therapeutic alternatives is the focus of many scientists worldwide. In this study, a QSAR-SVM classification model was developed by using the MATLAB numerical computation system and the molecular descriptors implemented in the Dragon software. The obtained parameters are adequate with accuracy of 88.63% for training set, 81.13% for cross-validation experiment and 81.15% for prediction set. In addition, its application domain was determined to guarantee the reliability of the predictions. Finally, the model was used to predict AChE inhibition by a group of quinazolinones and benzothiadiazine 1,1-dioxides obtained by chemical synthesis, resulting in 14 drug candidates with in silico activity comparable to acetylcholine.

Computational Modeling of Aldose Reductase Inhibitory Activity of Flavonoids Derivatives for Diabetic Complications Treatment

Background:

Diabetes mellitus is a chronic metabolic disease that constitutes a risk factor for patients infected by COVID-19. Aldose reductase (ALR2) is an enzyme that catalyzes the formation of sorbitol in the metabolism of glucose via polyols in diabetic patients and leads to a group of diabetic complications: cataracts, retinopathies, neuropathies, and nephropathies.

Introduction:

Inhibitors of this enzyme are therapeutic targets for the prophylaxis and treatment of these conditions. The aim of this work was to identify flavonoids isolated from medicinal plants, fruits, and vegetables as potential inhibitors of ALR2.

Methods:

In this study, using the MATLAB numerical computation system and the molecular descriptors implemented in the DRAGON software, a regression tree was developed, with an R2 of 0.953 and adequate parameters for its fit.

Results:

he model was validated to take into account internal and external validation procedures. Besides, the applicability domain was determined to guarantee the reliability of the predictions. Due to its good predictive power (R2 ext = 0.949), the model was used to predict the inhibition of ALR2 by flavonoids reported in dietary sources. The most promising flavonoids are Myricetin and Tricin (pIC50predicted = 7.296), which are within the application domain and meet drug-like properties for oral administration.

Conclusion:

Finally, we can conclude that the proposed tools are useful for the rapid and economical identification of flavonoid-based potential drug candidates against ALR2 in diabetic complications.

Machine learning approach to discovery of small molecules with potential inhibitory action against vasoactive metalloproteases

With the advancement of combinatorial chemistry and big data, drug repositioning has boomed. In this sense, machine learning and artificial intelligence techniques offer a priori information to identify the most promising candidates. In this study, we combine QSAR and docking methodologies to identify compounds with potential inhibitory activity of vasoactive metalloproteases for the treatment of cardiovascular diseases. To develop this study, we used a database of 191 thermolysin inhibitor compounds, which is the largest as far as we know. First, we use Dragon's molecular descriptors (0-3D) to develop classification models using Bayesian networks (Naive Bayes) and artificial neural networks (Multilayer Perceptron). The obtained models are used for virtual screening of small molecules in the international DrugBank database. Second, docking experiments are carried out for all three enzymes using the Autodock Vina program, to identify possible interactions with the active site of human metalloproteases. As a result, high-performance artificial intelligence QSAR models are obtained for training and prediction sets. These allowed the identification of 18 compounds with potential inhibitory activity and an adequate oral bioavailability profile, which were evaluated using docking. Four of them showed high binding energies for the three enzymes, and we propose them as potential dual ACE/NEP inhibitors for the control of blood pressure. In summary, the in silico strategies used here constitute an important tool for the early identification of new antihypertensive drug candidates, with substantial savings in time and money.

Computational identification of chemical compounds with potential anti-Chagas activity using a classification tree

Chagas disease is endemic to 21 Latin American countries and is a great public health problem in that region. Current chemotherapy remains unsatisfactory; consequently the need to search for new drugs persists. Here we present a new approach to identify novel compounds with potential anti-chagasic action. A large dataset of 584 compounds, obtained from the Drugs for Neglected Diseases initiative, was selected to develop the computational model. Dragon software was used to calculate the molecular descriptors and WEKA software to obtain the classification tree. The best model shows accuracy greater than 93.4% for the training set; the tree was also validated using a 10-fold cross-validation procedure and through a test set, achieving accuracy values over 90.5% and 92.2%, correspondingly. The values of sensitivity and specificity were around 90% in all series; also the false alarm rate values were under 10.5% for all sets. In addition, a simulated ligand-based virtual screening for several compounds recently reported as promising anti-chagasic agents was carried out, yielding good agreement between predictions and experimental results. Finally, the present work constitutes an example of how this rational computer-based method can help reduce the cost and increase the rate in which novel compounds are developed against Chagas disease.

Thorough evaluation of OECD principles in modelling of 1-[(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine derivatives using QSARINS

The human immunodeficiency virus is a lethal pathology considered as a worldwide problem. The search for new strategies for the treatment of this disease continues to be a great challenge in the scientific community. In this study, a series of 107 derivatives of 1-[(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine, previously evaluated experimentally against HIV-I reverse transcriptase, was used to model antiretroviral activity. A model of linear regression, implemented in the QSARINS software, was developed with a genetic algorithm for variable selection. The fit of its parameters was good and exhaustive validation, according to the OECD regulatory principles, was performed. Also, the applicability domain was established. In addition, its robustness (r ² = 0.84), stability (Q ² _LOO = 0.81; Q ² _LMO = 0.80) and good predictive power (r ² _EXT = 0.85) is proved. So, it was used to predict the antiretroviral activity of eight compounds obtained by rational drug design. Finally, it can be affirmed that the proposed tools allow the rapid and economic identification of potential antiretroviral drugs.

In Silico Assessment of ADME Properties: Advances in Caco-2 Cell Monolayer Permeability Modeling

One of the main goals of in silico Caco-2 cell permeability models is to identify those drug substances with high intestinal absorption in human (HIA). For more than a decade, several in silico Caco-2 models have been made, applying a wide range of modeling techniques; nevertheless, their capacity for intestinal absorption extrapolation is still doubtful. There are three main problems related to the modest capacity of obtained models, including the existence of inter- and/or intra-laboratory variability of recollected data, the influence of the metabolism mechanism, and the inconsistent in vitro-in vivo correlation (IVIVC) of Caco-2 cell permeability. This review paper intends to sum up the recent advances and limitations of current modeling approaches, and revealed some possible solutions to improve the applicability of in silico Caco-2 permeability models for absorption property profiling, taking into account the above-mentioned issues.

Discrete Fourier Transform-Based Multivariate Image Analysis: Application to Modeling of Aromatase Inhibitory Activity

We recently generalized the formerly alignment-dependent multivariate image analysis applied to quantitative structure-activity relationships (MIA-QSAR) method through the application of the discrete Fourier transform (DFT), allowing for its application to noncongruent and structurally diverse chemical compound data sets. Here we report the first practical application of this method in the screening of molecular entities of therapeutic interest, with human aromatase inhibitory activity as the case study. We developed an ensemble classification model based on the two-dimensional (2D) DFT MIA-QSAR descriptors, with which we screened the NCI Diversity Set V (1593 compounds) and obtained 34 chemical compounds with possible aromatase inhibitory activity. These compounds were docked into the aromatase active site, and the 10 most promising compounds were selected for in vitro experimental validation. Of these compounds, 7419 (nonsteroidal) and 89 201 (steroidal) demonstrated satisfactory antiproliferative and aromatase inhibitory activities. The obtained results suggest that the 2D-DFT MIA-QSAR method may be useful in ligand-based virtual screening of new molecular entities of therapeutic utility.

Prediction of aquatic toxicity of benzene derivatives using molecular descriptor from atomic weighted vectors

Several descriptors from atom weighted vectors are used in the prediction of aquatic toxicity of set of organic compounds of 392 benzene derivatives to the protozoo ciliate Tetrahymena pyriformis (log(IGC50)^-1). These descriptors are calculated using the MD-LOVIs software and various Aggregation Operators are examined with the aim comparing their performances in predicting aquatic toxicity. Variability analysis is used to quantify the information content of these molecular descriptors by means of an information theory-based algorithm. Multiple Linear Regression with Genetic Algorithms is used to obtain models of the structure-toxicity relationships; the best model shows values of Q²=0.830 and R²=0.837 using six variables. Our models compare favorably with other previously published models that use the same data set. The obtained results suggest that these descriptors provide an effective alternative for determining aquatic toxicity of benzene derivatives.

State of the Art Review and Report of New Tool for Drug Discovery

BACKGROUND

There are a great number of tools that can be used in QSAR/QSPR studies; they are implemented in several programs that are reviewed in this report. The usefulness of new tools can be proved through comparison, with previously published approaches. In order to perform the comparison, the most usual is the use of several benchmark datasets such as DRAGON and Sutherland's datasets.

METHODS

Here, an exploratory study of Atomic Weighted Vectors (AWVs), a new tool useful for drug discovery using different datasets, is presented. In order to evaluate the performance of the new tool, several statistics and QSAR/QSPR experiments are performed. Variability analyses are used to quantify the information content of the AWVs obtained from the tool, by means of an information theory-based algorithm.

RESULTS

Principal components analysis is used to analyze the orthogonality of these descriptors, for which the new MDs from AWVs provide different information from those codified by DRAGON descriptors (0-2D). The QSAR models are obtained for every Sutherland's dataset, according to the original division into training/test sets, by means of the multiple linear regression with genetic algorithm (MLR-GA). These models have been validated and compared favorably to several previously published approaches, using the same benchmark datasets.

CONCLUSION

The obtained results show that this tool should be a useful strategy for the QSAR/QSPR studies, despite its simplicity.

Machine learning-based models to predict modes of toxic action of phenols to Tetrahymena pyriformis

The phenols are structurally heterogeneous pollutants and they present a variety of modes of toxic action (MOA), including polar narcotics, weak acid respiratory uncouplers, pro-electrophiles, and soft electrophiles. Because it is often difficult to determine correctly the mechanism of action of a compound, quantitative structure-activity relationship (QSAR) methods, which have proved their interest in toxicity prediction, can be used. In this work, several QSAR models for the prediction of MOA of 221 phenols to the ciliated protozoan Tetrahymena pyriformis, using Chemistry Development Kit descriptors, are reported. Four machine learning techniques (ML), k-nearest neighbours, support vector machine, classification trees, and artificial neural networks, have been used to develop several models with higher accuracies and predictive capabilities for distinguishing between four MOAs. They showed global accuracy values between 95.9% and 97.7% and area under Receiver Operator Curve values between 0.978 and 0.998; additionally, false alarm rate values were below 8.2% for training set. In order to validate our models, cross-validation (10-folds-out) and external test-set were performed with good behaviour in all cases. These models, obtained with ML techniques, were compared with others previously reported by other researchers, and the improvement was significant.

Larvicidal activity prediction against Aedes aegypti mosquito using computational tools

BACKGROUND & OBJECTIVES

Aedes aegypti is an important vector for transmission of dengue, yellow fever, chikun- gunya, arthritis, and Zika fever. According to the World Health Organization, it is estimated that Ae. aegypti causes 50 million infections and 25,000 deaths per year. Use of larvicidal agents is one of the recommendations of health organizations to control mosquito populations and limit their distribution. The aim of present study was to deduce a mathematical model to predict the larvicidal action of chemical compounds, based on their structure.

METHODS

A series of different compounds with experimental evidence of larvicidal activity were selected to develop a predictive model, using multiple linear regression and a genetic algorithm for the selection of variables, implemented in the QSARINS software. The model was assessed and validated using the OECDs principles.

RESULTS

The best model showed good value for the determination coefficient (R2 = 0.752), and others parameters were appropriate for fitting (s = 0.278 and RMSEtr = 0.261). The validation results confirmed that the model hasgood robustness (Q2LOO = 0.682) and stability (R2-Q2LOO = 0.070) with low correlation between the descriptors (KXX = 0.241), an excellent predictive power (R2 ext = 0.834) and was product of a non-random correlation R2 Y-scr = 0.100).

INTERPRETATION & CONCLUSION

The present model shows better parameters than the models reported earlier in the literature, using the same dataset, indicating that the proposed computational tools are more efficient in identifying novel larvicidal compounds against Ae. aegypti.

Prediction of acute toxicity of phenol derivatives using multiple linear regression approach for Tetrahymena pyriformis contaminant identification in a median-size database

In this article, the modeling of inhibitory grown activity against Tetrahymena pyriformis is described. The 0-2D Dragon descriptors based on structural aspects to gain some knowledge of factors influencing aquatic toxicity are mainly used. Besides, it is done by some enlarged data of phenol derivatives described for the first time and composed of 358 chemicals. It overcomes the previous datasets with about one hundred compounds. Moreover, the results of the model evaluation by the parameters in the training, prediction and validation give adequate results comparable with those of the previous works. The more influential descriptors included in the model are: X3A, MWC02, MWC10 and piPC03 with positive contributions to the dependent variable; and MWC09, piPC02 and TPC with negative contributions. In a next step, a median-size database of nearly 8000 phenolic compounds extracted from ChEMBL was evaluated with the quantitative-structure toxicity relationship (QSTR) model developed providing some clues (SARs) for identification of ecotoxicological compounds. The outcome of this report is very useful to screen chemical databases for finding the compounds responsible of aquatic contamination in the biomarker used in the current work.

Tyrosinase enzyme: 1. An overview on a pharmacological target

The tyrosinase enzyme (EC 1.14.18.1) is an oxidoreductase inside the general enzyme classification and is involved in the oxidation and reduction process in the epidermis. These chemical reactions that the enzyme catalyzes are of principal importance in the melanogenesis process. This process of melanogenesis is related to the melanin formation, a heteropolymer of indolic nature that provides the different tonalities in the skin and helps to the protection from the ultraviolet radiation. However, a pigment overproduction, come up by the action of the tyrosinase, can cause different disorders in the skin related to the hyperpigmentation. Several studies mainly focused on the characteristics of the enzyme have been reported. In this work, an approximation to general aspects related to this enzyme is made. Besides, it is treated the researches that have been published in the part of the biochemical anatomy dealing with diseases associated with this protein (melanogenesis), its active place and its physiological states, the molecular mechanism, the methods carried out to detect the inhibitory activity, and the used substrates.

Comparative study to predict toxic modes of action of phenols from molecular structures

Quantitative structure-activity relationship models for the prediction of mode of toxic action (MOA) of 221 phenols to the ciliated protozoan Tetrahymena pyriformis using atom-based quadratic indices are reported. The phenols represent a variety of MOAs including polar narcotics, weak acid respiratory uncouplers, pro-electrophiles and soft electrophiles. Linear discriminant analysis (LDA), and four machine learning techniques (ML), namely k-nearest neighbours (k-NN), support vector machine (SVM), classification trees (CTs) and artificial neural networks (ANNs), have been used to develop several models with higher accuracies and predictive capabilities for distinguishing between four MOAs. Most of them showed global accuracy of over 90%, and false alarm rate values were below 2.9% for the training set. Cross-validation, complementary subsets and external test set were performed, with good behaviour in all cases. Our models compare favourably with other previously published models, and in general the models obtained with ML techniques show better results than those developed with linear techniques. We developed unsupervised and supervised consensus, and these results were better than our ML models, the results of rule-based approach and other ensemble models previously published. This investigation highlights the merits of ML-based techniques as an alternative to other more traditional methods for modelling MOA.

A novel approach to predict aquatic toxicity from molecular structure

The main aim of the study was to develop quantitative structure-activity relationship (QSAR) models for the prediction of aquatic toxicity using atom-based non-stochastic and stochastic linear indices. The used dataset consist of 392 benzene derivatives, separated into training and test sets, for which toxicity data to the ciliate Tetrahymena pyriformis were available. Using multiple linear regression, two statistically significant QSAR models were obtained with non-stochastic (R2=0.791 and s=0.344) and stochastic (R2=0.799 and s=0.343) linear indices. A leave-one-out (LOO) cross-validation procedure was carried out achieving values of q2=0.781 (scv=0.348) and q2=0.786 (scv=0.350), respectively. In addition, a validation through an external test set was performed, which yields significant values of Rpred2 of 0.762 and 0.797. A brief study of the influence of the statistical outliers in QSAR's model development was also carried out. Finally, our method was compared with other approaches implemented in the Dragon software achieving better results. The non-stochastic and stochastic linear indices appear to provide an interesting alternative to costly and time-consuming experiments for determining toxicity.

Estimation of ADME Properties in Drug Discovery: Predicting Caco-2 Cell Permeability Using Atom-Based Stochastic and Non-stochastic Linear Indices

The in vitro determination of the permeability through cultured Caco-2 cells is the most often-used in vitro model for drug absorption. In this report, we use the largest data set of measured P(Caco-2), consisting of 157 structurally diverse compounds. Linear discriminant analysis (LDA) was used to obtain quantitative models that discriminate higher absorption compounds from those with moderate-poorer absorption. The best LDA model has an accuracy of 90.58% and 84.21% for training and test set. The percentage of good correlation, in the virtual screening of 241 drugs with the reported values of the percentage of human intestinal absorption (HIA), was greater than 81%. In addition, multiple linear regression models were developed to predict Caco-2 permeability with determination coefficients of 0.71 and 0.72. Our method compares favorably with other approaches implemented in the Dragon software, as well as other methods from the international literature. These results suggest that the proposed method is a good tool for studying the oral absorption of drug candidates.

Atom-based stochastic and non-stochastic 3D-chiral bilinear indices and their applications to central chirality codification

Non-stochastic and stochastic 2D bilinear indices have been generalized to codify chemical structure information for chiral drugs, making use of a trigonometric 3D-chirality correction factor. In order to evaluate the effectiveness of this novel approach in drug design we have modeled the angiotensin-converting enzyme inhibitory activity of perindoprilate's sigma-stereoisomers combinatorial library. Two linear discriminant analysis models, using non-stochastic and stochastic linear indices, were obtained. The models had shown an accuracy of 95.65% for the training set and 100% for the external prediction set. Next the prediction of the sigma-receptor antagonists of chiral 3-(3-hydroxyphenyl)piperidines by multiple linear regression analysis was carried out. Two statistically significant QSAR models were obtained when non-stochastic (R(2)=0.953 and s=0.238) and stochastic (R(2)=0.961 and s=0.219) 3D-chiral bilinear indices were used. These models showed adequate predictive power (assessed by the leave-one-out cross-validation experiment) yielding values of q(2)=0.935 (s(cv)=0.259) and q(2)=0.946 (s(cv)=0.235), respectively. Finally, the prediction of the corticosteroid-binding globulin binding affinity of steroids set was performed. The obtained results are rather similar to most of the 3D-QSAR approaches reported so far. The validation of this method was achieved by comparison with previous reports applied to the same data set. The non-stochastic and stochastic 3D-chiral linear indices appear to provide a very interesting alternative to other more common 3D-QSAR descriptors.

Predicting antitrichomonal activity: A computational screening using atom-based bilinear indices and experimental proofs

Existing Trichomonas vaginalis therapies are out of reach for most trichomoniasis people in developing countries and, where available, they are limited by their toxicity (mainly in pregnant women) and their cost. New antitrichomonal agents are needed to combat emerging metronidazole-resistant trichomoniasis and reduce the side effects associated with currently available drugs. Toward this end, atom-based bilinear indices, a new TOMOCOMD-CARDD molecular descriptor, and linear discriminant analysis (LDA) were used to discover novel, potent, and non-toxic lead trichomonacidal chemicals. Two discriminant functions were obtained with the use of non-stochastic and stochastic atom-type bilinear indices for heteroatoms and H-bonding of heteroatoms. These atomic-level molecular descriptors were calculated using a weighting scheme that includes four atomic labels, namely atomic masses, van der Waals volumes, atomic polarizabilities, and atomic electronegativities in Pauling scale. The obtained LDA-based QSAR models, using non-stochastic and stochastic indices, were able to classify correctly 94.51% (90.63%) and 93.41% (93.75%) of the chemicals in training (test) sets, respectively. They showed large Matthews' correlation coefficients (C); 0.89 (0.79) and 0.87 (0.85), for the training (test) sets, correspondingly. The result of predictions on the 15% full-out cross-validation test also evidenced the robustness and predictive power of the obtained models. In addition, canonical regression analyses corroborated the statistical quality of these models (R(can) of 0.749 and of 0.845, correspondingly); they were also used to compute biological activity canonical scores for each compound. On the other hand, a close inspection of the molecular descriptors included in both equations showed that several of these molecular fingerprints are strongly interrelated with each other. Therefore, these models were orthogonalized using the Randić orthogonalization procedure. These classification functions were then applied to find new lead antitrichomonal agents and six compounds were selected as possible active compounds by computational screening. The designed compounds were synthesized and tested for in vitro activity against T. vaginalis. Out of the six compounds that were designed, and synthesized, three molecules (chemicals VA5-5a, VA5-5c, and VA5-12b) showed high to moderate cytocidal activity at the concentration of 10 microg/ml, other two compounds (VA5-8pre and VA5-8) showed high cytocidal and cytostatic activity at the concentration of 100 microg/ml and 10 microg/ml, correspondingly, and the remaining chemical (compound VA5-5e) was inactive at these assayed concentrations. Nonetheless, these compounds possess structural features not seen in known trichomonacidal compounds and thus can serve as excellent leads for further optimization of antitrichomonal activity. The LDA-based QSAR models presented here can be considered as a computer-assisted system that could potentially significantly reduce the number of synthesized and tested compounds and increase the chance of finding new chemical entities with antitrichomonal activity.

Atom-based 3D-chiral quadratic indices. Part 2: Prediction of the corticosteroid-binding globulinbinding affinity of the 31 benchmark steroids data set

A quantitative structure-activity relationship (QSAR) study to predict the relative affinities of the steroid 'benchmark' data set to the corticosteroid-binding globulin (CBG) is described. It is shown that the 3D-chiral quadratic indices closely correlate with the measured CBG affinity values for the 31 steroids. The calculated descriptors were correlated with biological data through multiple linear regressions. Two statistically significant models were obtained when non-stochastic (R = 0.924 and s = 0.46) as well as stochastic (R = 0.929 and s = 0.46) 3D-chiral quadratic indices were used. A leave-one-out (LOO) approach to model validation is used here; the best results obtained in the cross-validation procedure with non-stochastic (q2 = 0.781) and stochastic (q2 = 0.735) 3D-chiral quadratic indices are better or similar to most of the 3D-QSAR approaches reported so far. These results support the idea that the 3D-chiral quadratic indices may be helpful in prediction of the corticosteroid-binding affinity for new compounds.

3D-chiral Atom, Atom-type, and Total Non-stochastic and Stochastic Molecular Linear Indices and their Applications to Central Chirality Codification

Non-stochastic and stochastic 2D linear indices have been generalized to codify chemical structure information for chiral drugs, making use of a trigonometric 3D-chirality correction factor. These descriptors circumvent the inability of conventional 2D non-stochastic [Y. Marrero-Ponce. J. Chem. Inf. Comp., Sci. l 44 (2004) 2010] and stochastic [Y. Marrero-Ponce, et al. Bioorg. Med. Chem., 13 (2005) 1293] linear indices to distinguish sigma-stereoisomers. In order to test the potential of this novel approach in drug design we have modelled the angiotensin-converting enzyme inhibitory activity of perindoprilate's sigma-stereoisomers combinatorial library. Two linear discriminant analysis models, using non-stochastic and stochastic linear indices, were obtained. The models showed an accuracy of 100% and 96.65% for the training set; and 88.88% and 100% in the external test set, respectively. Canonical regression analysis corroborated the statistical quality of these models (R(can) of 0.78 and of 0.77) and was also used to compute biology activity canonical scores for each compound. After that, the prediction of the sigma-receptor antagonists of chiral 3-(3-hydroxyphenyl)piperidines by linear multiple regression analysis was carried out. Two statistically significant QSAR models were obtained when non-stochastic (R2 = 0.982 and s = 0.157) and stochastic (R2 = 0.941 and s = 0.267) 3D-chiral linear indices were used. The predictive power was assessed by the leave-one-out cross-validation experiment, yielding values of q2 = 0.982 (s(cv) = 0.186) and q2 = 0.90 (s(cv) = 0.319), respectively. Finally, the prediction of the corticosteroid-binding globulin binding affinity of steroids set was performed. The best results obtained in the cross-validation procedure with non-stochastic (q2 = 0.904) and stochastic (q2 = 0.88) 3D-chiral linear indices are rather similar to most of the 3D-QSAR approaches reported so far. The validation of this method was achieved by comparison with previous reports applied to the same data set. The non-stochastic and stochastic 3D-chiral linear indices appear to provide an interesting alternative to other more common 3D-QSAR descriptors.

Protein linear indices of the ‘macromolecular pseudograph α-carbon atom adjacency matrix’ in bioinformatics. Part 1: Prediction of protein stability effects of a complete set of alanine substitutions in Arc repressor

A novel approach to bio-macromolecular design from a linear algebra point of view is introduced. A protein's total (whole protein) and local (one or more amino acid) linear indices are a new set of bio-macromolecular descriptors of relevance to protein QSAR/QSPR studies. These amino-acid level biochemical descriptors are based on the calculation of linear maps on Rn[f k(xmi):Rn-->Rn] in canonical basis. These bio-macromolecular indices are calculated from the kth power of the macromolecular pseudograph alpha-carbon atom adjacency matrix. Total linear indices are linear functional on Rn. That is, the kth total linear indices are linear maps from Rn to the scalar R[f k(xm):Rn-->R]. Thus, the kth total linear indices are calculated by summing the amino-acid linear indices of all amino acids in the protein molecule. A study of the protein stability effects for a complete set of alanine substitutions in the Arc repressor illustrates this approach. A quantitative model that discriminates near wild-type stability alanine mutants from the reduced-stability ones in a training series was obtained. This model permitted the correct classification of 97.56% (40/41) and 91.67% (11/12) of proteins in the training and test set, respectively. It shows a high Matthews correlation coefficient (MCC=0.952) for the training set and an MCC=0.837 for the external prediction set. Additionally, canonical regression analysis corroborated the statistical quality of the classification model (Rcanc=0.824). This analysis was also used to compute biological stability canonical scores for each Arc alanine mutant. On the other hand, the linear piecewise regression model compared favorably with respect to the linear regression one on predicting the melting temperature (tm) of the Arc alanine mutants. The linear model explains almost 81% of the variance of the experimental tm (R=0.90 and s=4.29) and the LOO press statistics evidenced its predictive ability (q2=0.72 and scv=4.79). Moreover, the TOMOCOMD-CAMPS method produced a linear piecewise regression (R=0.97) between protein backbone descriptors and tm values for alanine mutants of the Arc repressor. A break-point value of 51.87 degrees C characterized two mutant clusters and coincided perfectly with the experimental scale. For this reason, we can use the linear discriminant analysis and piecewise models in combination to classify and predict the stability of the mutant Arc homodimers. These models also permitted the interpretation of the driving forces of such folding process, indicating that topologic/topographic protein backbone interactions control the stability profile of wild-type Arc and its alanine mutants.

Tomocomd-Cardd, a novel approach for computer-aided ? rational? drug design: I. Theoretical and experimental assessment of a promising method for computational screening and in silico design of new anthelmintic compounds

In this work, the TOMOCOMD-CARDD approach has been applied to estimate the anthelmintic activity. Total and local (both atom and atom-type) quadratic indices and linear discriminant analysis were used to obtain a quantitative model that discriminates between anthelmintic and non-anthelmintic drug-like compounds. The obtained model correctly classified 90.37% of compounds in the training set. External validation processes to assess the robustness and predictive power of the obtained model were carried out. The QSAR model correctly classified 88.18% of compounds in this external prediction set. A second model was performed to outline some conclusions about the possible modes of action of anthelmintic drugs. This model permits the correct classification of 94.52% of compounds in the training set, and 80.00% of good global classification in the external prediction set. After that, the developed model was used in virtual in silico screening and several compounds from the Merck Index, Negwer's handbook and Goodman and Gilman were identified by models as anthelmintic. Finally, the experimental assay of one organic chemical (G-1) by an in vivo test coincides fairly well (100%) with model predictions. These results suggest that the proposed method will be a good tool for studying the biological properties of drug candidates during the early state of the drug-development process.

Atom, atom-type and total molecular linear indices as a promising approach for bioorganic and medicinal chemistry: theoretical and experimental assessment of a novel method for virtual screening and rational design of new lead anthelmintic

Helminth infections are a medical problem in the world nowadays. In this paper a novel atom-level chemical descriptor has been applied to estimate the anthelmintic activity. Total and local linear indices and linear discriminant analysis were used to obtain a quantitative model that discriminates between anthelmintic and non-anthelmintic drug-like compounds. The discriminant model has an accuracy of 90.11% in the training set, with a high Matthews' correlation coefficient (MCC=0.80). To assess the robustness and predictive power of the obtained model, internal (leave-n-out) and external validation process was performed. The QSAR model correctly classified 88.55% of compounds in this external prediction set, yielding a MCC of 0.77. Another LDA model was carried out to outline some conclusions about the possible modes of action of anthelmintic drugs. It has an accuracy of 93.50% in the training set, and 80.00% in the external prediction set. After that, the developed model was used in the virtual--in silico--screening and several compounds from the Merck Index, Negwer's Handbook and Goodman and Gilman were identified by the model as anthelmintic. Finally, the experimental assay of an organic chemical (a furylethylene derivative) by an in vivo test permits us to carry out an assessment of the model. An accuracy of 100% with the theoretical predictions was observed. These results suggest that the proposed method will be a good tool for studying the biological properties of drug candidates during the early state of the drug-development process.