2D Autocorrelation modeling of the activity of trihalobenzocycloheptapyridine analogues as farnesyl protein transferase inhibitors

Development and Validation of Predictive Quantitative Structure-Activity Relationship Models for Estrogenic Activities of Hydroxylated Polychlorinated Biphenyls

Disruption of the endocrine system by hydroxylated polychlorinated biphenyls (OH-PCBs) is hypothesized, among other potential mechanisms, to be mediated via nuclear receptor binding. Due to the high cost and lengthy time required to produce high-quality experimental data, empirical data to support the nuclear receptor binding hypothesis are in short supply. In the present study, two quantitative structure-activity relationship models were developed for predicting the estrogenic activities of OH-PCBs. Findings revealed that model I (for the estrogen receptor α dataset) contained five two-dimensional (2D) descriptors belonging to the classes autocorrelation, Burden modified eigenvalues, chi path, and atom type electrotopological state, whereas model II (for the estrogen receptor β dataset) contained three 2D and three 3D descriptors belonging to the classes autocorrelation, atom type electrotopological state, and Radial Distribution Function descriptors. The internal and external validation metrics reported for models I and II indicate that both models are robust, reliable, and suitable for predicting the estrogenic activities of untested OH-PCB congeners. Environ Toxicol Chem 2023;42:823-834. © 2023 SETAC.

Computational Modeling to Explain Why 5,5-Diarylpentadienamides are TRPV1 Antagonists

Several years ago, the crystallographic structures of the transient receptor potential vanilloid 1 (TRPV1) in the presence of agonists and antagonists were reported, providing structural information about its chemical activation and inactivation. TRPV1’s activation increases the transport of calcium and sodium ions, leading to the excitation of sensory neurons and the perception of pain. On the other hand, its antagonistic inactivation has been explored to design analgesic drugs. The interactions between the antagonists 5,5-diarylpentadienamides (DPDAs) and TRPV1 were studied here to explain why they inactivate TRPV1. The present work identified the structural features of TRPV1–DPDA complexes, starting with a consideration of the orientations of the ligands inside the TRPV1 binding site by using molecular docking. After this, a chemometrics analysis was performed (i) to compare the orientations of the antagonists (by using LigRMSD), (ii) to describe the recurrent interactions between the protein residues and ligand groups in the complexes (by using interaction fingerprints), and (iii) to describe the relationship between topological features of the ligands and their differential antagonistic activities (by using a quantitative structure–activity relationship (QSAR) with 2D autocorrelation descriptors). The interactions between the DPDA groups and the residues Y511, S512, T550, R557, and E570 (with a recognized role in the binding of classic ligands), and the occupancy of isoquinoline or 3-hydroxy-3,4-dihydroquinolin-2(1H)-one groups of the DPDAs in the vanilloid pocket of TRPV1 were clearly described. Based on the results, the structural features that explain why DPDAs inactivate TRPV1 were clearly exposed. These features can be considered for the design of novel TRPV1 antagonists.

Synthesis of diN-Substituted Glycyl-Phenylalanine Derivatives by Using Ugi Four Component Reaction and Their Potential as Acetylcholinesterase Inhibitors

Ugi four component reaction (Ugi-4CR) isocyanide-based multicomponent reactions were used to synthesize diN-substituted glycyl-phenylalanine (diNsGF) derivatives. All of the synthesized compounds were characterized by spectroscopic and spectrometric techniques. In order to evaluate potential biological applications, the synthesized compounds were tested in computational models that predict the bioactivity of organic molecules by using only bi-dimensional molecular information. The diNsGF derivatives were predicted as cholinesterase inhibitors. Experimentally, all of the synthesized diNsGF derivatives showed moderate inhibitory activities against acetylcholinesterase (AChE) and poor activities against butyrylcholinesterase (BuChE). Compound 7a has significant activity and selectivity against AChE, which reveals that the diNsGF scaffold could be improved to reach novel candidates by combining other chemical components of the Ugi-4CR in a high-throughput combinatorial screening experiment. Molecular docking experiments of diNsGF derivatives inside AChE suggest that these compounds placed the phenylalanine group at the peripheral site of AChE. The orientations and chemical interactions of diNsGF derivatives were analyzed, and the changeable groups were identified for future exploration of novel candidates that could lead to the improvement of diNsGF derivative inhibitory activities.

Imidazole-containing farnesyltransferase inhibitors: 3D quantitative structure-activity relationships and molecular docking

One of the most promising anticancer and recent antimalarial targets is the heterodimeric zinc-containing protein farnesyltransferase (FT). In this work, we studied a highly diverse series of 192 Abbott-initiated imidazole-containing compounds and their FT inhibitory activities using 3D-QSAR and docking, in order to gain understanding of the interaction of these inhibitors with FT to aid development of a rational strategy for further lead optimization. We report several highly significant and predictive CoMFA and CoMSIA models. The best model, composed of CoMFA steric and electrostatic fields combined with CoMSIA hydrophobic and H-bond acceptor fields, had r (2) = 0.878, q (2) = 0.630, and r (pred) (2) = 0.614. Docking studies on the statistical outliers revealed that some of them had a different binding mode in the FT active site based on steric bulk and available active site space, explaining why the predicted activities differed from the experimental activities.

Modeling of acetylcholinesterase inhibition by tacrine analogues using Bayesian-regularized Genetic Neural Networks and ensemble averaging

Acetylcholinesterase inhibition was modeled for a set of 136 tacrine analogues using Bayesian-regularized Genetic Neural Networks (BRGNNs). In the BRGNN approach the Bayesian-regularization avoids overtraining/overfitting and the genetic algorithm (GA) allows exploring an ample pool of 3D-descriptors. The predictive capacity of our selected model was evaluated by averaging multiple validation sets generated as members of diverse-training set neural network ensembles (NNEs). The ensemble averaging provides reliable statistics. When 40 members are assembled, the NNE provides a reliable measure of training and test set R values of 0.921 and 0.851 respectively. In other respects, the ability of the nonlinear selected GA space for differentiating the data was evidenced when the total data set was well distributed in a Kohonen Self-Organizing Map (SOM). The location of the inhibitors in the map facilitates the analysis of the connection between compounds and serves as a useful tool for qualitative predictions.

Protein radial distribution function (P-RDF) and Bayesian-Regularized Genetic Neural Networks for modeling protein conformational stability: Chymotrypsin inhibitor 2 mutants

Development of novel computational approaches for modeling protein properties is a main goal in applied Proteomics. In this work, we reported the extension of the radial distribution function (RDF) scores formalism to proteins for encoding 3D structural information with modeling purposes. Protein-RDF (P-RDF) scores measure spherical distributions on protein 3D structure of 48 amino acids/residues properties selected from the AAindex data base. P-RDF scores were tested for building predictive models of the change of thermal unfolding Gibbs free energy change (DeltaDeltaG) of chymotrypsin inhibitor 2 upon mutations. In this sense, an ensemble of Bayesian-Regularized Genetic Neural Networks (BRGNNs) yielded an optimum nonlinear model for the conformational stability. The ensemble predictor described about 84% and 70% variance of the data in training and test sets, respectively.

2D Autocorrelation, CoMFA, and CoMSIA modeling of protein tyrosine kinases' inhibition by substituted pyrido[2,3-d]pyrimidine derivatives

2D Autocorrelation, comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA) were undertaken for a series of substituted pyrido[2,3-d]pyrimidine derivatives to correlate platelet-derived growth factor receptor (PDGFR), fibroblast growth factor receptor (FGFR), and c-Src tyrosine kinases' inhibition with 2D and 3D structural properties of 22 known compounds. QSAR models with considerable internal as well as external predictive ability were obtained. The relevant 2D autocorrelation descriptors for modeling each protein tyrosine kinase (PTK) inhibitory activity were selected by genetic algorithm (GA) and multiple linear regression (MLR) approach. The 2D autocorrelation space brings different descriptors for each PTK inhibition and suggests the atomic properties relevant for the inhibitors to interact with each PTK active site. CoMFA and CoMSIA were developed with a focus on interpretative ability using coefficient contour maps. CoMSIA produced significantly better results for all correlations. The results indicate a strong correlation between the inhibitory activity of the modeled compounds and the hydrophobic and H-bond donor fields around them.

QSAR Study on the Relaxant Agents from Some Mexican Medicinal Plants and Synthetic Related Organic Compounds

Quantitative Structure-Activity Relationship studies were performed to describe and predict the antispasmodic activity of some molecules isolated from Mexican Medicinal Flora as well as for some synthetic ones based on stilbenoid bioisosteres. The relaxant activity of these molecules was taken from experiments on rat and guinea-pig ileum tissues. Given that there is some evidence of species-specific on the relaxant effects, two data sets were proposed, one for rat ileum and the other for guinea-pig ileum. These data were statistically treated in order to find a Quantitative Structure-Activity Relationship model that could describe the corresponding biological models. The goodness of prediction for the best models was measured in terms of the Leave-One-Out Cross-Validation R(2) (LOO q(2)) and the correlation coefficients of regressions through the origin (RTO R(2)0). Results show that papaverine activity could not be used as reference in rat ileum tests; however, this molecule can be used as a good reference molecule in guinea-pig ileum tests. Our study shows that MATS5p and R8m+ descriptors are the most important descriptors in predicting the rat ileum activity and that atomic polarizability is the main atomic property. On the other hand, the R3u GETAWAY descriptor turns out to be important in predicting the guinea-pig ileum activity where the influence/distance of substituents on these molecules could describe the observed activity.

Proteometric study of ghrelin receptor function variations upon mutations using amino acid sequence autocorrelation vectors and genetic algorithm-based least square support vector machines

Functional variations on the human ghrelin receptor upon mutations have been associated with a syndrome of short stature and obesity, of which the obesity appears to develop around puberty. In this work, we reported a proteometrics analysis of the constitutive and ghrelin-induced activities of wild-type and mutant ghrelin receptors using amino acid sequence autocorrelation (AASA) approach for protein structural information encoding. AASA vectors were calculated by measuring the autocorrelations at sequence lags ranging from 1 to 15 on the protein primary structure of 48 amino acid/residue properties selected from the AAindex database. Genetic algorithm-based multilinear regression analysis (GA-MRA) and genetic algorithm-based least square support vector machines (GA-LSSVM) were used for building linear and non-linear models of the receptor activity. A genetic optimized radial basis function (RBF) kernel yielded the optimum GA-LSSVM models describing 88% and 95% of the cross-validation variance for the constitutive and ghrelin-induced activities, respectively. AASA vectors in the optimum models mainly appeared weighted by hydrophobicity-related properties. However, differently to the constitutive activity, the ghrelin-induced activity was also highly dependent of the steric features of the receptor.

QSAR modeling of matrix metalloproteinase inhibition by N-hydroxy-alpha-phenylsulfonylacetamide derivatives

The main molecular features which determine the selectivity of a set of 80 N-hydroxy-alpha-phenylsulfonylacetamide derivatives (HPSAs) in the inhibition of three matrix metalloproteinases (MMP-1, MMP-9, and MMP-13) have been identified by using linear and nonlinear predictive models. The molecular information has been encoded in 2D autocorrelation descriptors, obtained from different weighting schemes. The linear models were built by multiple linear regression (MLR) combined with genetic algorithm (GA), and a robust QSAR mapping paradigm. The Bayesian-regularized genetic neural network (BRGNN) was employed for nonlinear modeling. In such approaches each model could have its own set of input variables. All models were predictive according to internal and external validation experiments; but the best results correspond to nonlinear ones. The 2D autocorrelation space brings different descriptors for each MMP inhibition, and suggests the atomic properties relevant for the inhibitors to interact with each MMP active site. On the basis of the current results, the reported models have the potential to discover new potent and selective inhibitors and bring useful molecular information about the ligand specificity for MMP S(1)(') and S(2)(') subsites.

Amino acid sequence autocorrelation vectors and bayesian-regularized genetic neural networks for modeling protein conformational stability: Gene V protein mutants

Development of novel computational approaches for modeling protein properties from their primary structure is the main goal in applied proteomics. In this work, we reported the extension of the autocorrelation vector formalism to amino acid sequences for encoding protein structural information with modeling purposes. Amino acid sequence autocorrelation (AASA) vectors were calculated by measuring the autocorrelations at sequence lags ranging from 1 to 15 on the protein primary structure of 48 amino acid/residue properties selected from the AAindex data base. A total of 720 AASA descriptors were tested for building predictive models of the change of thermal unfolding Gibbs free energy change (delta deltaG) of gene V protein upon mutation. In this sense, ensembles of Bayesian-regularized genetic neural networks (BRGNNs) were used for obtaining an optimum nonlinear model for the conformational stability. The ensemble predictor described about 88% and 66% variance of the data in training and test sets respectively. Furthermore, the optimum AASA vector subset not only helped to successfully model unfolding stability but also well distributed wild-type and gene V protein mutants on a stability self-organized map (SOM), when used for unsupervised training of competitive neurons.

Ensembles of Bayesian-regularized genetic neural networks for modeling of acetylcholinesterase inhibition by huprines

Acetylcholinesterase inhibition was modeled for a set of huprines using ensembles of Bayesian-regularized Genetic Neural Networks. In the Bayesian-regularized Genetic Neural Network approach the Bayesian regularization avoids overfitted regressions and the genetic algorithm allows exploring a wide pool of three-dimensional descriptors. The predictive capacity of our selected model was evaluated by averaging multiple validation sets generated as members of neural network ensembles. When 60 members are assembled, the neural network ensemble provides a reliable measure of training and test set R(2)-values of 0.945 and 0.850 respectively. In other respects, the ability of the nonlinear selected genetic algorithm space for differentiate the data were evidenced when total data set was well distributed in a Kohonen self-organizing map. The analysis of the self-organizing map zones allows establishing the main structural features differentiated by our vectorial space.

Quantitative Structure?Activity Relationship Modeling of Growth Hormone Secretagogues Agonist Activity of some Tetrahydroisoquinoline 1-Carboxamides

Growth hormone secretagogue agonist activities for a data set of 45 tetrahydroisoquinoline 1-carboxamides were modeled using several kinds of molecular descriptors from dragon software. A linear model with six variables selected from a large pool of two-dimensional descriptors described 80% of cross-validation data variance. Similar results were found for a model obtained from a pool of three-dimensional descriptors. Size and hydrophilicity-related atomic properties such as mass, polarizability, and van der Waals volume were determined to be the most relevant for the differential growth hormone secretagogue agonist activities of the compounds studied. In addition, Artificial Neural Networks were trained using optimum variables from the linear models; however, they were found to overfit the data and resulted in similar or lower predictive power.

Bayesian-regularized genetic neural networks applied to the modeling of non-peptide antagonists for the human luteinizing hormone-releasing hormone receptor

Bayesian-regularized genetic neural networks (BRGNNs) were used to model the binding affinity (IC(50)) for 128 non-peptide antagonists for the human luteinizing hormone-releasing hormone (LHRH) receptor using 2D spatial autocorrelation vectors. As a preliminary step, a linear dependence was established by multiple linear regression (MLR) approach, selecting the relevant descriptors by genetic algorithm (GA) feature selection. The linear model showed to fit the training set (N=102) with R(2)=0.746, meanwhile BRGNN exhibited a higher value of R(2)=0.871. Beyond the improvement of training set fitting, the BRGNN model overcame the linear one by being able to describe 85% of test set (N=26) variance in comparison with 73% the MLR model. Our non-linear QSAR model illustrates the importance of an adequate distribution of atomic properties represented in topological frames and reveals the electronegativities, masses and polarizabilities as the most influencing atomic properties in the structures of the heterocycles under analysis for having an appropriate LHRH antagonistic activity. Furthermore, the ability of the non-linear selected variables for differentiating the data was evidenced when total data set was well distributed in a Kohonen self-organizing map (SOM).

Linear and nonlinear QSAR study of N-hydroxy-2-[(phenylsulfonyl)amino]acetamide derivatives as matrix metalloproteinase inhibitors

The inhibitory activity (IC50) toward matrix metalloproteinases (MMP-1, MMP-2, MMP-3, MMP-9, and MMP-13) of N-hydroxy-2-[(phenylsulfonyl)amino]acetamide derivatives (HPSAAs) has been successfully modeled using 2D autocorrelation descriptors. The relevant molecular descriptors were selected by linear and nonlinear genetic algorithm (GA) feature selection using multiple linear regression (MLR) and Bayesian-regularized neural network (BRANN) approaches, respectively. The quality of the models was evaluated by means of cross-validation experiments and the best results correspond to nonlinear ones (Q2>0.7 for all models). Despite the high correlation between the studied compound IC50 values, the 2D autocorrelation space brings different descriptors for each MMP inhibition. On the basis of these results, these models contain useful molecular information about the ligand specificity for MMP S'1, S1, and S'2 pockets.

2D Autocorrelation modeling of the negative inotropic activity of calcium entry blockers using Bayesian-regularized genetic neural networks

Negative inotropic potency of 60 benzothiazepine-like calcium entry blockers (CEBs), Diltiazem analogs, was successfully modeled using Bayesian-regularized genetic neural networks (BRGNNs) and 2D autocorrelation vectors. This approach yielded reliable and robust models whilst by means of a linear genetic algorithm (GA) search routine no multilinear regression model was found describing more than 50% of the training set. On the contrary, the optimum neural network predictor with five inputs described about 84% and 65% variances of 50 randomly selected training and test sets. Autocorrelation vectors in the nonlinear model contained information regarding 2D spatial distributions on the CEB structure of van der Waals volumes, electronegativities, and polarizabilities. However, a sensitivity analysis of the network inputs pointed out to the electronegativity and polarizability 2D topological distributions at substructural fragments of sizes 3 and 4 as the most relevant features governing the nonlinear modeling of the negative inotropic potency.

Genetic neural network modeling of the selective inhibition of the intermediate-conductance Ca2+-activated K+ channel by some triarylmethanes using topological charge indexes descriptors

Selective inhibition of the intermediate-conductance Ca(2+)-activated K(+ )channel (IK (Ca)) by some clotrimazole analogs has been successfully modeled using topological charge indexes (TCI) and genetic neural networks (GNNs). A neural network monitoring scheme evidenced a highly non-linear dependence between the IK (Ca) blocking activity and TCI descriptors. Suitable subsets of descriptors were selected by means of genetic algorithm. Bayesian regularization was implemented in the network training function with the aim of assuring good generalization qualities to the predictors. GNNs were able to yield a reliable predictor that explained about 97% data variance with good predictive ability. On the contrary, the best multivariate linear equation with descriptors selected by linear genetic search, only explained about 60%. In spite of when using the descriptors from the linear equations to train neural networks yielded higher fitted models, such networks were very unstable and had relative low predictive ability. However, the best GNN BRANN 2 had a Q ( 2 ) of LOO of cross-validation equal to 0.901 and at the same time exhibited outstanding stability when calculating 80 randomly constructed training/test sets partitions. Our model suggested that structural fragments of size three and seven have relevant influence on the inhibitory potency of the studied IK (Ca) channel blockers. Furthermore, inhibitors were well distributed regarding its activity levels in a Kohonen self-organizing map (KSOM) built using the inputs of the best neural network predictor.