Three-dimensional quantitative structure-activity relationship analyses of piperidine-based CCR5 receptor antagonists

Prediction of PKCθ inhibitory activity using the Random Forest Algorithm

This work is devoted to the prediction of a series of 208 structurally diverse PKCθ inhibitors using the Random Forest (RF) based on the Mold(2) molecular descriptors. The RF model was established and identified as a robust predictor of the experimental pIC(50) values, producing good external R(2) (pred) of 0.72, a standard error of prediction (SEP) of 0.45, for an external prediction set of 51 inhibitors which were not used in the development of QSAR models. By using the RF built-in measure of the relative importance of the descriptors, an important predictor-the number of group donor atoms for H-bonds (with N and O)-has been identified to play a crucial role in PKCθ inhibitory activity. We hope that the developed RF model will be helpful in the screening and prediction of novel unknown PKCθ inhibitory activity.

Tracking chemical processing pathways in combinatorial polymer libraries via data mining

Changes in the molecular structure and composition of interpenetrating polymer networks (IPNs) can be used to tailor their properties. While the properties of IPNs are typically different than polymer blends, a clear understanding of the impact of changing polymerization sequence on the physical properties and the corresponding molecular bonding is needed. To address this issue, a data mining approach is used to identify the change with polymerization sequence of tensile and rheological properties of acrylate-epoxy IPNs. The experimental approach used to study the molecular structure is high throughput Fourier transform infrared (FTIR) spectroscopy. Analysis of the FTIR spectra of IPNs synthesized with different polymerization sequences leads to an understanding of the molecular bonding responsible for the tensile and rheological properties. From the interpretation of the wavenumber bands and associated molecular bonds, we find that the polymerization sequence most affects hydrogen bonding and aromatic ring bond energies. This work defines the relationships between chemistry, structure, processing, and properties of the IPN samples.

Quantitative structure-activity relationships studies of CCR5 inhibitors and toxicity of aromatic compounds using gene expression programming

Quantitative structure-activity relationship (QSAR) study of chemokine receptor 5 (CCR5) binding affinity of substituted 1-(3,3-diphenylpropyl)-piperidinyl amides and ureas and toxicity of aromatic compounds have been performed. The gene expression programming (GEP) was used to select variables and produce nonlinear QSAR models simultaneously using the selected variables. In our GEP implementation, a simple and convenient method was proposed to infer the K-expression from the number of arguments of the function in a gene, without building the expression tree. The results were compared to those obtained by artificial neural network (ANN) and support vector machine (SVM). It has been demonstrated that the GEP is a useful tool for QSAR modeling.

Comparative docking study of anibamine as the first natural product CCR5 antagonist in CCR5 homology models

Anibamine, a novel pyridine quaternary alkaloid recently isolated from Aniba sp., has been found to effectively bind to the chemokine receptor CCR5 with an IC(50) at 1 microM in competition with (125)I-gp120, an HIV viral envelope protein binding to CCR5 with high affinity. Since CCR5, a G-protein-coupled receptor, is an essential coreceptor for the human immunodeficiency virus type I (HIV-1) entry to host cells, a CCR5 antagonist that inhibits the cellular entry of HIV-1 provides a new therapy choice for the treatment of HIV. Anibamine provides a novel structural skeleton that is remarkably different from all lead compounds previously identified as CCR5 antagonists. Here, we report comparative docking studies of anibamine with several other known CCR5 antagonists in two CCR5 homology models built based on the crystal structures of bovine rhodopsin and human beta(2)-adrenergic receptor. The binding pocket of anibamine has some common features shared with other high affinity CCR5 antagonists, suggesting that they may bind in similar binding sites and/or modes. At the same time, several unique binding features of anibamine were identified, and it will likely prove beneficial in future molecular design of novel CCR5 antagonists based on the anibamine scaffold.

3D-QSAR studies of substituted 1-(3, 3-diphenylpropyl)-piperidinyl amides and ureas as CCR5 receptor antagonists

3D-QSAR studies on the derivatives of 1-(3,3-diphenylpropyl)-piperidinyl amide and urea as CCR5 receptor antagonists were performed by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA) methods to rationalize the structural requirements responsible for the inhibitory activity of these compounds. The global minimum energy conformer of the template molecule, the most active and pharmacokinetically stable molecule of the series, was obtained by systematic search and used to build structures of the molecules in the dataset. The best predictions for the CCR5-receptor were obtained with the CoMFA standard model (q (2) = 0.787, r (2) = 0.962) and CoMSIA model combined steric, electrostatic and hydrophobic fields (q (2) = 0.809, r (2) = 0.951). The predictive ability of CoMFA and CoMSIA were determined using a test set of 12 compounds giving predictive correlation coefficients of 0.855 and 0.83, respectively, indicating good predictive power. Further, the robustness of the model was verified by bootstrapping analysis. The contour maps produced by the CoMFA and CoMSIA models were used to identify the structural features relevant to the biological activity in this series. Based on the CoMFA and CoMSIA analysis, we have identified some key features in the series that are responsible for CCR5 antagonistic activity which may be used to design more potent 1-(3,3-diphenylpropyl)-piperidinyl derivatives and predict their activity prior to synthesis.

3D-QSAR study of hallucinogenic phenylalkylamines by using CoMFA approach

The three-dimensional quantitative structure-activity relationship (3D-QSAR) has been studied on 90 hallucinogenic phenylalkylamines by the comparative molecular field analysis (CoMFA). Two conformations were compared during the modeling. Conformation I referred to the amino group close to ring position 6 and conformation II related to the amino group trans to the phenyl ring. Satisfactory results were obtained by using both conformations. There were still differences between the two models. The model based on conformation I got better statistical results than the one about conformation II. And this may suggest that conformation I be preponderant when the hallucinogenic phenylalkylamines interact with the receptor. To further confirm the predictive capability of the CoMFA model, 18 compounds with conformation I were randomly selected as a test set and the remaining ones as training set. The best CoMFA model based on the training set had a cross-validation coefficient q (2) of 0.549 at five components and non cross-validation coefficient R (2) of 0.835, the standard error of estimation was 0.219. The model showed good predictive ability in the external test with a coefficient R (pre) (2) of 0.611. The CoMFA coefficient contour maps suggested that both steric and electrostatic interactions play an important role. The contributions from the steric and electrostatic fields were 0.450 and 0.550, respectively.

CCR5 receptor antagonists: Discovery and SAR study of guanylhydrazone derivatives

High throughput screening (HTS) led to the identification of the guanylhydrazone of 2-(4-chlorobenzyloxy)-5-bromobenzaldehyde as a CCR5 receptor antagonist. Initial modifications of the guanylhydrazone series indicated that substitution of the benzyl group at the para-position was well tolerated. Substitution at the 5-position of the central phenyl ring was critical for potency. Replacement of the guanylhydrazone group led to the discovery of a novel series of CCR5 antagonists.

Investigation of substituent effect of 1-(3,3-diphenylpropyl)-piperidinyl phenylacetamides on CCR5 binding affinity using QSAR and virtual screening techniques

A linear quantitative-structure activity relationship model is developed in this work using Multiple Linear Regression Analysis as applied to a series of 51 1-(3,3-diphenylpropyl)-piperidinyl phenylacetamides derivatives with CCR5 binding affinity. For the selection of the best variables the Elimination Selection-Stepwise Regression Method (ES-SWR) is utilized. The predictive ability of the model is evaluated against a set of 13 compounds. Based on the produced QSAR model and an analysis on the domain of its applicability, the effects of various structural modifications on biological activity are investigated. The study leads to a number of guanidine derivatives with significantly improved predicted activities.

QSAR analyses of 3-(4-benzylpiperidin-1-yl)-N-phenylpropylamine derivatives as potent CCR5 antagonists

CCR5 receptor binding affinity of a series of 3-(4-benzylpiperidin-1-yl)propylamine congeners was subjected to QSAR study using the linear free energy related (LFER) model of Hansch. Appropriate indicator variables encoding different group contributions and different physicochemical variables such as hydrophobicity (pi), electronic (Hammett sigma), and steric (molar refractivity, STERIMOL values) parameters of phenyl ring substituents of the compounds were used as predictor variables. The Hansch analysis explores the importance of the lipophilicity and electron-donating substituents for the binding affinity. However, this method could not give more insight into the structure-activity relationships because of the diverse molecular features in the data set. 3D-QSAR analyses of the same data set using Molecular Shape Analysis (MSA), Receptor Surface Analysis (RSA), and Molecular Field Analysis (MFA) techniques were also performed. The best model with acceptable statistical quality was derived from the MSA, which showed the importance of the relative negative charge (RNCG): substituents with a high RNCG value have more binding affinity than the unsubstituted piperidine and phenyl (R1 position) congeners. The relative negative charge surface area (RNCS) is detrimental (e.g. R2 = 3,4-Cl2) for the activity. An increase in the length of the molecule in the Z dimension (Lz) is conducive (e.g. R3 = sulfonylmorpholino), while an increase in the area of the molecular shadow in the XZ plane (Sxz) is detrimental (e.g. R1 = N-c-hexylmethyl-5-oxopyrrolidin-3-yl) for the binding affinity. The presence of a chiral center makes the molecule less active (e.g. R1 = N-methyl-5-oxopyrrolidin-3-yl). An increase in the van der Waals area, the molecular volume, and the difference between the volume of the individual molecule and the shape reference compound are conducive (e.g. R3 = (CH3)2NSO2-) for the binding affinity. Substituents with higher JursFPSA_2 values (fractional charged partial surface area) like the N-methylsulfonylpiperidin-4-yl (R1 position) group have better binding affinity than the substituents such as 4-chlorophenylamino (R1 position). Unsubstituted piperidines (R1 position) with less JursFNSA_1 values have lower binding affinity than the 4-chlorophenyl substituted compounds. The MFA derived equation shows interaction energies at different grid points, while the RSA model shows the importance of hydrophobicity and charge at different regions of the molecules. The models were validated through the leave-one-out, leave-15%-out, and leave-25%-out cross-validation techniques. The developed models were also subjected to a randomization test (99% confidence level). Although the MSA derived models had excellent statistical qualities both for the training as well as test sets, RSA and MFA results for the test sets are not comparable statistically with the MSA derived models.

Quantitative structure-activity relationships for small non-peptide antagonists of CXCR2: indirect 3D approach using the frontal polygon method

The chemokine receptor, CXCR2, plays an important role in recruiting granulocytes to sites of inflammation and has been proposed as an important therapeutic target. A number of CXCR2 antagonists have been synthesized and evaluated; however, quantitative structure-activity relationship (QSAR) models have not been developed for these molecules. Most CXCR2 antagonists can be grouped into four related categories: N,N'-diphenylureas, nicotinamide N-oxides, quinoxalines, and triazolethiols. Based on these categories, we developed a QSAR model for 59 nonpeptide antagonists of CXCR2 using a partial 3D comparison of the antagonists with local fingerprints obtained from rigid and flexible fragments of the molecules. Each compound was represented by calculated structural descriptors that encoded atomic charge, molar refraction, hydrophobicity, and geometric features. We obtained good conventional R(2) coefficients, high leave-one-out cross-validated values for the whole dataset (R(cv)(2)=0.785), as well as for the dataset divided into subsets of triazolethiol derivatives (R(cv)(2)=0.821) and joint subset of N'-diphenylureas, nicotinamide N-oxides, N,N'-diphenylureas, and quinoxaline derivatives and quinoxalines derivatives (R(cv)(2)=0.766), indicating a good predictive ability and robustness of the model. Additionally, charge distribution was found to be a significant contributor in modeling whole dataset. Using our model, structural fragments (submolecules) responsible for the antagonist activity were also identified. These data suggest the QSAR models developed here may be useful in guiding the design of CXCR2 antagonists from molecular fragments.

Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor

In the present study, we have used an approach combining protein structure modeling, molecular dynamics (MD) simulation, automated docking, and 3D QSAR analyses to investigate the detailed interactions of CCR5 with their antagonists. Homology modeling and MD simulation were used to build the 3D model of CCR5 receptor based on the high-resolution X-ray structure of bovine rhodopsin. A series of 64 CCR5 antagonists, 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes, were docked into the putative binding site of the 3D model of CCR5 using the docking method, and the probable interaction model between CCR5 and the antagonists were obtained. The predicted binding affinities of the antagonists to CCR5 correlate well with the antagonist activities, and the interaction model could be used to explain many mutagenesis results. All these indicate that the 3D model of antagonist-CCR5 interaction is reliable. Based on the binding conformations and their alignment inside the binding pocket of CCR5, three-dimensional structure-activity relationship (3D QSAR) analyses were performed on these antagonists using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods. Both CoMFA and CoMSIA provide statistically valid models with good correlation and predictive power. The q(2)(r(cross)(2)) values are 0.568 and 0.587 for CoMFA and CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were not included in the training set. Mapping these models back to the topology of the active site of CCR5 leads to a better understanding of antagonist-CCR5 interaction. These results suggest that the 3D model of CCR5 can be used in structure-based drug design and the 3D QSAR models provide clear guidelines and accurate activity predictions for novel antagonist design.