Docking Analysis and Multidimensional Hybrid QSAR Model of 1,4-Benzodiazepine-2,5-Diones as HDM2 Antagonists

authors:

avatar Yujie Dai a , * , avatar Nan Chen a , avatar Qiang Wang a , avatar Heng Zheng b , avatar Xiuli Zhang c , avatar Shiru Jia a , avatar Lilong Dong d , avatar Dacheng Feng e

Key Laboratory of Industrial Fermentation Microbiology (Tianjin University of Science and Technology), Ministry of Education, College of Bioengineering, Tianjin University of Science and Technology, Tianjin 300457, P.R. China.
School of Life Science and Technology, China Pharmaceutical University, Nanjing 210009, P.R. China.
Department of Biochemistry, University of Missouri-Columbia, Columbia, MO 65211, USA.
School of Pharmaceutical Sciences, Hebei Medical University, Shijiazhuang 050017, P.R. China.
College of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, P.R. China.
Iranian Journal of Pharmaceutical Research : IJPR: Vol.11, issue 3; 807-830
article type: Original Article
received: May , 2011
accepted: April , 2011
DOI: https://doi.org/10.22037/ijpr.2012.1185

how to cite: Dai Y, Chen N, Wang Q, Zheng H, Zhang X, et al. Docking Analysis and Multidimensional Hybrid QSAR Model of 1,4-Benzodiazepine-2,5-Diones as HDM2 Antagonists. Iran J Pharm Res. 2012;11(3):e125809. https://doi.org/10.22037/ijpr.2012.1185.

Abstract

The inhibitors of p53-HDM2 interaction are attractive molecules for the treatment of wild-type p53 tumors. In order to search more potent HDM2 inhibitors, docking operation with CDOCKER protocol in Discovery Studio 2.1 (DS2.1) and multidimensional hybrid quantitative structure-activity relationship (QSAR) studies through the physiochemical properties obtained from DS2.1 and E-Dragon 1.0 as descriptors, have been performed on 59 1,4-benzodiazepine- 2,5-diones which have p53-HDM2 interaction inhibitory activities. The docking results indicate that π-π interaction between the imidazole group in HIS96 and the aryl ring at 4-N of 1,4-benzodiazepine-2,5-dione may be one of the key factors for the combination of ligands with HDM2. Two QSAR models were obtained using genetic function approximation (GFA) and genetic partial least squares (G/PLS) based on the descriptors obtained from DS2.1 and E-dragon 1.0, respectively. The best model can explain 85.5% of the variance (R 2adj ) while it could predict 81.7% of the variance (R 2cv ). With this model, the bioactivities of some new compounds were predicted.

Introduction

The p53 is a multifunctional protein that regulates genes which can induce either cell cycle arrest or apoptosis (1-3). In human cells, p53 activity is normally modulated by the human double minute-2 (HDM2) protein, the homologue of murine double minute-2 (MDM2) protein (4). HDM2 binds to and blocks the p53 transactivation domain, inhibiting its transcriptional activity. The dissociation of the binding between these two proteins is therefore an attractive therapeutic target for the treatment of wild-type p53 tumors (5). This action can be implemented through neutralizing antibodies (6), p53 peptides (7) or protein fusions (8, 9), which leads to the control of proliferation. Additionally, antisense HDM2 oligonucleotides increase both p53 levels and activity by reducing HDM2 protein levels (10, 11). Since small molecular antagonists of HDM2 are easily synthesized, some investigations on them have been reported, but many of these molecules have low potency and limited cellular activity (12-14). Most recently, some reports describing a series of 1,4-benzodiazepine-2,5-dione compounds (BDPs) that act as the potent antagonists of the HDM2–p53 interaction (15). A library of BDPs was designed by some investigators using directed diversity method (16-18). These compounds were synthesized and screened with thermofluor screening technology (19) and further evaluated with a fluorescence polarization (FP) assay. In order to further understand the antagonizing mechanism of these compounds with HDM2-p53 interaction, Grasberger (11) et al. studied the crystal structure of BDPs with HDM2. The crystal structure confirmed that the BDPs occupy the p53 peptide binding site and they can mimic the α-helix of p53 peptide and may represent a promising scaffold to develop HDM2 antagonists. The three substituted aryl rings of the BDP overlay very closely with PHE19, TRP23 and LEU26 in the p53 binding pocket of HDM2 (20) and some structure-activity relationship (SAR) studies were deduced and explained (15). Recently, Wang et al. (21) conducted a benzodiazepinedione/peptide-based 3D-QSAR analysis using the comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), however, their molecular poses of ligands in the pocket of HDM2 from docking and crystal structure were only considered a little by taking a receptor-guided consensus dynamics alignment since the poses of the alignments method are different from those of docking results. On the other hand, there is no explicit quantitative equation reported for the QSAR of BDPs as inhibitors of the p53-HDM2 although some quantitative structure-activity relation (QSAR) models have been reported for other HDM2 inhibitors (22, 23). Whereas, BDPs are a kind of most potent HDM2 inhibitors, a high quality equation as QSAR model is necessary for finding new BDPs as potent HDM2 inhibitors. It is also important for new HDM2 inhibitors from other structural compounds.

Lately, with the development of computer hardware and computational theory, many molecule structural properties can be obtained through the computational process. Some software products provide these functions, such as Discovery studio (24), MOE (25), SYBYL (26) and E-Dragon (27), etc. For example, Discovery Studio 2.1 can calculate more than 100 kinds of molecular properties and more than 1600 molecular properties can be obtained from E-Dragon 1.0. These properties can be used expediently as descriptors to construct a QSAR model for drug screening. In 2009, Roys (28) calculated 2D and 3D properties of 116 diverse classes of aromatase inhibitors. With the assistance of genetic function approximation (GFA) and genetic partial least squares (G/PLS), QSAR models for non-steroid compounds were constructed using these 2D and 3D descriptors. The best (externally) predictive model was a GFA model with spline option using combined set (2D and 3D) descriptors and its predictive R2 (Rpred2) reached 0.687. Dai et al. (29) used the computed molecular properties and docking scores (CDOCKER interaction energy) constructed a 2D-3D hybrid QSAR model for steroid aromatase inhibitors with the correlation coefficient of R2 = 0.729 and Leave-one-out Cross-Validation of R2 = 0.667. In these two studies, GFA was successfully employed to select the best combined set of descriptors from large number of physiochemical properties for the construction of high quality QSAR models. In order to construct a quantitative predictive model for the bioactivity of new 1,4-benzodiazepine-2,5-dione compounds as HDM2 antagonists, two QSAR models were built in this paper based on the computational properties of 59 1,4-benzodiazepine-2,5-dione compounds with known bioactivities. In addition, the molecular docking of BDPs to the HDM2 cleft was conducted using CDOCKER docking protocol in Discovery Studio 2.1(DS2.1) and the interaction of these compounds was analyzed with LigPlot (30).

Experimental

Dataset and descriptors

The 1,4-benzodiazepine-2,5-dione chiral ligands with known inhibitory activities upon HDM2-p53 interaction reported in the literature (4, 15, 31, 32) were used as the model dataset for the present study (Table 1). The experimental protocol for the determination of inhibitory activities for all the compounds was the FP assay. The inhibitory values in IC50 (μM) were converted to the logarithmic scale pIC50 (mM) before being used for the subsequent QSAR analyses as the response variable. Various physiochemical properties of the ligands obtained using the protocol of Calculate Molecular Properties in DS2.1 and E-Dragon 1.0 respectively, were selected as candidate descriptors used for QSAR construction. These physiochemical properties include 1D (Element counts, functional group counts) 2D (AlogP, LogD, Num_StereoAtoms, Num_Rings, Molecular_Weight, Molecular_SurfaceArea, Num_H_Acceptors, Num_H_Donors, Num_RotatableBonds, 2D autocorrelations, topological descriptors such as BIC, V_ADJ_equ, V_DIST_equ, CHI_3_P, CHI_V_3_P, CHI_1, CHI_3_C, CIC, IAC_Mean, IAC_Total, IC and SIC, etc.) and 3D (3D-NoRSE descriptors, Dipole, Jurs descriptors, shadow indices and Molecular_Volume, etc.) parameters. All the definition of these descriptors can be seen in the Help of DS2.1 and E-Dragon 1.0 software. For the calculation of 3D descriptors, multiple conformations of each molecule were generated using the Generate Conformations protocol with optimal search as a conformational search method. Each conformer was subjected to an energy minimization procedure using smart minimizer under the consistent force field (CFF) to generate the lowest energy conformation for each structure. The charges were calculated according to the Gasteiger method.

Docking process of 1,4-benzodiazepine-2,5-diones to HDM2 protein

The crystal structure of HDM2 in complex with a benzodiazepine ligand ((S)-2-(4-chlorophenyl)-2-((S)-3-(4-chlorophenyl)-7-iodo-2,5-dioxo-2,3-dihydro-1H-benzo[e][1,4]diazepin-4(5H)-yl) acetic acid, ligand 25 in Table 1, Figure 1) has been obtained from the RCSB protein data bank (http://www.pdb.org, ID:1T4E).

The docking studies were conducted using CDOCKER of Receptor-ligand Interactions protocol in DS 2.1(33). The ligands and the HDM2 protein were all pretreated initially. For ligand preparation, the 3D structures of all the 1,4-benzodiazepine-2,5-diones were generated with ChemBioOffice 2008 (34) and optimized with AM1 method. Thereafter, the ligands were treated with Prepare Ligands protocol in DS 2.1. All the duplicate structures were removed and the options for ionization change, tautomer generation, isomer generation and 3D generator have been set true. For protein preparation, the hydrogen atoms were added with the pH in the range of 6.5-8.5 with DS2.1. HDM2 protein was defined as a total receptor and the site sphere was built with diameter of 10 Å based on the ligand 25. Then, the preexisting ligand 25 was removed. From the receptor-ligand interaction section of DS 2.1, CDOCKER was chosen and all the default operating parameters were used unless pre-declared. During the docking process, a freshly prepared ligand (compound from the dataset in Table 1) prepared by us was placed. CHARMm was selected as the force field. The molecular docking was performed with a simulated annealing method to minimize the CDOCKER energy (ECD) for obtaining an optimum pose. This method was independently described by Scott Kirkpatrick et al. in 1983 (35) and by Vlado Cerny in 1985 (36). The method is a technique suitable for the optimization problems of large scale, especially ones where a desired global extremum is hidden among many, poorer, local extrema. The name comes from annealing in metallurgy, a procedure involving heating and controlled cooling of a material to increase the crystal orders and reduce their defects. The heat makes the atoms to become unstuck from their initial positions with a local minimum of internal energy and move randomly through states of higher energy. The slow cooling procedure gives them ample time for redistribution of the atoms to find configurations with lower internal energy than the initial one. Similar to this physical process, each step of the simulated annealing algorithm changes the current solution through a random “nearby” solution with a probability that relies on both the difference between the corresponding function values and a global parameter T (named as temperature), which is gradually decreased during the process. The changes of current solution depend on T and are almost random when T is large, but increasingly “downhill” (for a minimization value) as T comes to zero. The chance for “uphill” moves potentially and prevents the method from becoming stuck at a local optimal value. In this study, the heating steps were set as 2000 with 700, heating target temperature. The cooling steps were set as 5000 with 300 as the cooling target temperature. Ten molecular docking poses saved for each ligand were ranked according to -CDOCKER energy (-ECD). The pose with the highest -ECD was chosen as the most suitable pose for the subsequent pose analysis. After the end of molecular docking, the interactions of the docked receptor (HDM2) with ligand were analyzed with LigPlot.

QSAR model development

For a comparatively small quantity of samples with large available variables, the number of variables (descriptors) should be reduced far smaller than that of samples so that the obtained equation has the statistical meaning. But it is difficult to choose which property is more suitable as the descriptor to build QSAR models. Recently, this problem can be solved with genetic function approximation (GFA) technique. The principles about GFA can be seen elsewhere (37, 38). It involves the multivariate adaptive regression algorithm combining with the genetic algorithm (GA) to evolve the population of equations (each containing only a subset of variables) that best fit the training set data. With this method, a series of potential solutions to a problem (the population of organisms) are derived and tested repeatedly until an approximate optimal solution is found. In this study, 59 1,4-benzodiazepine-2,5-dione chiral compounds in Table 1 were selected as the training set. There are more than 120 2D and 3D physiochemical properties obtained from the protocol of Calculate Molecular Properties in DS 2.1 and 1621 physiochemical properties obtained from E-Dragon 1.0. But only a subset of these properties is statistically significant according to the correlation with the compounds activities. Some redundant parameters from the protocol of Calculate Molecular Properties in DS 2.1 were omitted using GFA protocol in DS 2.1 in order to save time and space and to decrease the complexity of the QSAR model (The reduction of 1D-3D physiochemical descriptors from E-Dragon 1.0 to build the QSAR model was conducted using the Material Studio 4.0 (39)).

Statistical quality evaluation and model validation

For a successful QSAR model, it should be robust enough to be capable of making accurate and reliable predictions of the biological activities, thus, the predictive capacity of the developed QSAR models from the training set should be validated. There are several methods to confirm the quality of QSAR model. In our experiments, Friedman lack-of-fit (LOF) was used for the selection of GFA-derived equations, while for G/PLS equations, the correlation coefficient R2 and the adjusted R2(R 2adj ), were taken as objective functions to select an equation. Leave-one-out cross-validation R2 (R 2cv ) was employed to validate the predictivity of generated QSAR model equations. The LOF [40] is designed to resist over-fitting, which is a problem often encountered in constructing statistical models. Since the number of descriptors available in HDM2 inhibitor QSAR analysis normally exceeds the number of observations (training set compounds), the ability to prevent over-fitting of GFA is critical to the successful construction of a statistically significant QSAR model. The smoothing factor was set to 0.5, which controls the model size, and GFA was used to optimize the QSAR models having different numbers of descriptor terms. For a given smoothing factor, the optimization of a QSAR model was considered to be realized when descriptor usage became constant and independent of an increasing number of crossover operations. All the descriptors in the QSAR trial descriptor pool were used as linear terms during the GFA to generate QSAR models.

Results and Discussion

Analysis of the binding site

To understand the molecular details of HDM2 binding through the BDPs, the binding site view was made using the DS visualizer 2.5 (Figure 2).

HDM2 binding model of Ligand 25 (25 is displayed in stick. H, white; C, gray; O, red; N, blue; Cl, green; iodo, and purple; PDB ID: 1T4E). (A) The specific residues surrounding 25 within 4Å (HIS96 is displayed in stick. π-π interaction is displayed in orange line and H-bond is in green dash line); (B) Comparison of Ligand 25 poses from X-ray crystal structure (red) and from CDOCKER docking (gray) and HIS96 is also shown in stick as location reference (The image was made with Chimera 1.3(41)).

HDM2 binding model of Ligand 25 (25 is displayed in stick. H, white; C, gray; O, red; N, blue; Cl, green; iodo, and purple; PDB ID: 1T4E). (A) The specific residues surrounding 25 within 4Å (HIS96 is displayed in stick. π-π interaction is displayed in orange line and H-bond is in green dash line); (B) Comparison of Ligand 25 poses from X-ray crystal structure (red) and from CDOCKER docking (gray) and HIS96 is also shown in stick as location reference (The image was made with Chimera 1.3(41)).

Figure 2

The specific cleft to which the ligand binds (within 4 Å), contains both polar (GLY16, SER17, GLY58, ILE61, TYR67, GLN72, HIS73, HIS96, ILE99, TYR100) and nonpolar (LEU54, PHE55, MET62, VAL75, VAL93, PHE86, PHE97) amino acids and the inhibitor (contains phenyls A, B and C) occupies the same pockets as the peptide side-chains PHE 19, TRP 23, and LEU 26 of p53 (Figure 2A) ( 20). The HDM2 interaction with the ligand 25 was analyzed with LigPlot and it was found that the nonspecific hydrophobic contacts are largely responsible for their interaction. The hydrophobic contacts of HIS 96, LEU 54, GLY 16 with the ring of phenyl C, ILE 99, Leu 57, GLY 58 with phenyl B, and ILE 61, MET 62 with phenyl A are shown in Figure 3A. There are also two hydrogen bonds between the carboxyl group of Ligand 25 and hydroxyl group of SER17, iodine of Ligand 25 and hydroxyl group of GLN 72 respectively. It can be concluded that the binding cleft of HDM2 is predominantly hydrophobic and largely nonspecific Van Der Waals contacts are responsible for the interaction between the ligand and HDM2 hydrophobic pocket (11).

The 1,4-benzodiazepine-2,5-dione compounds (BDPs) used to construct QSAR models are shown in Table 1. All the ligands listed in Table 1 are different from each other in R1 to R4, which lead to the different conformations adopted in the cleft of HDM2, therefore, the interactions with HDM2 are different and the bioactivities are varied.

Table 1

Structural features of the 1,4-benzodiazepine-2,5-diones (4, 15, 27, 28) having HDM2 inhibitory activity

Analysis of the optimum molecular docking poses of different ligands

The CDOCKER algorithm was applied to dock the ligands listed in Table 1. This docking technique needs a site sphere surrounding the ligand (the radius was set to 10 Å in this study) and this technique was tested with the ligand 25. Our docking pose of freshly prepared model of ligand 25 with CDOCKER also corroborates with that in crystal structure (RMSD = 0.425) indicating the reliability of this docking procedure (Figure 2B).

Each ligand in Table 1 was docked as described previously and the pose with the highest -ECD was considered as the optimum pose for each ligand. For the ligands in the top bioactivity range, such as ligands 25, 51 and 55 (IC50 < 0.5 μM), there exists π-π interaction between the imidazole group in HIS96 and the aryl ring in R3 (the occasion for ligand 25 is seen in Figure 2A). It can be considered as one of the key factors for the combination of the ligands with HDM2. Although the pocket of HDM2 is mainly hydrophobic, the additional hydrogen bonds can be formed. By introducing an amine functional group in the ortho position of aryl ring in R3, an additional hydrogen bond with VAL93 is formed for ligand 51 or 55, which is consistent with the previous prediction in literature (Figure 3B) (33).The extra hydrogen bond formation increased the inhibitory potency of these ligands.

The interaction of the ligands with HDM2 protein (The image was made with LigPlot (30). A: ligand 25, B: ligand 51).

The interaction of the ligands with HDM2 protein (The image was made with LigPlot (30). A: ligand 25, B: ligand 51).

Figure 3

The R1 changes in fatty acid, alkyl ether, alkyl (acyl) morpholine, alkyl (acyl) piperazine etc. (ligands 48-59). They are introduced in the ligand merely for improving the water-solubility of the ligands (33). From the docking position, we can see that these groups have nearly little contact with the acceptor protein (Figure 4A shows the morpholine group in ligand 51 stretches out of the HDM2 cleft), and in turn, have little effect on enhancing the ligand affinity. On the other hand, they improve the solvent effect, which has the negative effect on the ligands’ affinity to the target.

Poses of different ligands in the active site of HDM2. (A) Surface representation of Hdm-ligand 51 complex the protein (the morpholine group of ligand 51 is out of the HDM2 cleft). (B) Docking pose difference between ligands 51 (violet) and 52 (grey). Instead of the aryl ring in R3 of 51, the aryl ring in R2 formed the π-π interaction with the imidazole group of HIS96 in ligand 52. (C) Comparison of docking poses of ligand 25 (red, I in purple) with 35 (blue, Cl in green). (D) docking poses comparison of ligand 12 (in line) with ligand 15 (in stick).

Poses of different ligands in the active site of HDM2. (A) Surface representation of Hdm-ligand 51 complex the protein (the morpholine group of ligand 51 is out of the HDM2 cleft). (B) Docking pose difference between ligands 51 (violet) and 52 (grey). Instead of the aryl ring in R3 of 51, the aryl ring in R2 formed the π-π interaction with the imidazole group of HIS96 in ligand 52. (C) Comparison of docking poses of ligand 25 (red, I in purple) with 35 (blue, Cl in green). (D) docking poses comparison of ligand 12 (in line) with ligand 15 (in stick).

Figure 4

The chirality change at R2 and α-C of R3 also alters the bioactivity seriously, for example, ligands 49 and 50, 51 and 52, 53 and 54, are different from each other in R2 or α-C of R3 correspondingly (Figure 4B shows the docking pose difference between the ligands 51 and 52). Instead of the aryl ring in R3 of 51 (in violet), the aryl ring in R2 formed the π-π interaction with the imidazole group of HIS96 in ligand 52 (grey). This conformation makes the extra hydrogen bond existed in ligands 51 which can’t be formed for ligand 52 that may lower the affinity of the ligand with HDM2.

Figure 4C shows the comparison of docking poses of I substituted (ligand 25) and Cl substituted (ligand 35) phenyl A ligands. From the comparison of the docking results, it was found that iodine group in phenyl A is suitable for the cleft space. When the iodine group of phenyl A was replaced by other halogen, cyanogen, amide, alkyl or alkyne groups, the ligand sizes could not match the HDM2 cleft as exactly as iodine, therefore, the interaction between the ligand and HDM2 was decreased and as a result, the HDM2 inhibitory activities were lowered.

The FP activity of ligands can be changed by substituted groups in aryl groups of R2 and R3. The absence of a substituent in phenyl group resulted in a dramatic loss of potency as exemplified by ligand 28. Proper substituent size in the para-position of the phenyl group will make the ligand match the HDM2 cleft more suitably and the results show that the ethyl group (ligand 10), Cl, CF3 and OCF3 (ligands 12-14) are optimal for these substituents. All three aryl groups are in the cleft of HDM2 for these ligands. Substitution at either the ortho- or meta-positions of phenyl group in R2 results in a sharp loss in activity due to the spatial clag generated for HDM2-liangd interaction (ligands 14 and 15). Figure 4D shows the docking poses comparison of ligand 12 with ligand 15. The phenyl group A in 4-10 is driven out of the cleft. To get a better understanding of the effect of 1,4-benzodiazepine-2,5-diones compounds 3D structure on their bioactivities, the docked positions with the minimized CDOCKER interaction energies (ECD) of the 59 compounds are shown in Figure 5.

Docked positions with the minimized CDOCKER interaction energies (ECD) of the 59 1,4-benzodiazepine-2,5-diones compounds

Docked positions with the minimized CDOCKER interaction energies (ECD) of the 59 1,4-benzodiazepine-2,5-diones compounds

Figure 5

Although the structural information can be obtained by viewing the ligand poses in the cleft of HDM2 and there are many SAR analyses everywhere (4, 15), these SARs cannot directly give the bioactivity values. The QSAR can give an explicit predicted value for a new compound and QSAR studies are vital to enable the prioritization of analogues resulting from iterative virtual screening and thus, the design of focused small molecule libraries around active ligands. Recently, there are some reports on QSAR study on MDM2 inhibitors (22, 42) and most recently Wang et al. (12) built some QSAR models using CoMFA and CoMSIA for BDPs and peptides as HDM2 inhibitors. However, there are no explicit quantitative equation model reported on 1,4-benzodiazepine-2,5-dione compounds though they are a kind of most potent entities as HMD2 inhibitors. A quantitative equation QSAR model of HMD2 inhibitors is more convenient to be used for the discovery of new antitumor drugs based on p53-HDM2 interaction. Therefore in this article, we used the protocol of Calculate Molecular Properties in DS 2.1 and E-Dragon 1.0 to calculate various physiochemical properties as the candidate descriptors for the construction of QSAR models of the substituted 1,4-benzodiazepine-2,5-dione compounds. Since the number of properties is larger than the kinds of obtained ligands, it should be reduced to a much small number than that of ligands in order to obtain an equation statistically. Due to the large number of properties, GFA method was employed to eliminate the properties with less relation to the bioactivity of ligand and keep the subset consisting of the most correlation to the potency of ligands. In this study, the physiochemical properties of BDPs from DS 2.1 were directly treated with the GFA protocol of this software, but the molecular descriptors from E-Dragon were treated with GFA in Material Studio 4.0. The linear term was used for the development of models with Friedman LOF smoothness parameter of 0.5 and the population size of 1000. The obtained QSAR model using linear term was then further treated with G/PLS and the model on the descriptors from DS 2.1 is as follows:

pIC50 = 14.568 + 0.388 LogD - 0.166 Num_ RotatableBonds - 0.670 Num_StereoAtoms + 0.00278 V_DIST_equ - 1.446 CHI_1 - 0.0471 Dipole_X + 0.230 Shadow_Xlength - 0.0328 Shadow_XZ (Equation 1)

The sample number N = 59, LOF = 0.198, R2 = 0.750, R 2adj= 0.672 = R 2cv ,0.712 = R 2adj, F = 19.54. The standardized regression coefficient for each variable is 0.624, - 0.450, - 0.477, 4.01, - 3.30, - 0.546, 0.492 and -0.394 respectively.

In our study, R2, R 2cv , R 2adj and F were used to evaluate the regression model. Equation 1 can explain 71.2% of the variance (R 2 cv ) while it could predict 67.2% of the variance (R 2cv ). F > F(a = 0.05) = 2.13 shows that the model is in the confidence interval of 95%. It can be seen from Equation 1 that LogD, V_DIST_equ and Shadow_Xlength have positive contribution to the bioactivity of the ligands, however, Num_ RotatableBonds, Num_StereoAtoms, Dipole_X, CHI_1, and Shadow_XZ have negative effect on the bioactivities of the ligands. The relative importance of the descriptors is in the following order according to their standardized regression coefficients:

V_DIST_equ > CHI_1 >> LogD > Dipole_X > Shadow_Xlength > Num_StereoAtoms > Num_RotatableBonds > Shadow_XZ

From this order, we can see that V_DIST_equ and CHI_1 play the key role in determining the bioactivity of ligands, however, since CHI_1 and Shadow_XZ have roughly the same change tendency as V_DIST_equ, their effect on the bioactivities of ligands is mainly counteracted by V_DIST_equ. Although ligands 26, 27, 56 and 57 have comparatively high values of CHI_1 and Shadow_XZ, they possess significant inhibitory activity due to the high V_DIST_equ values. Ligands 23, 26-29 with R1 substituents have the high LogD and the ligands 52, 54 and 55 with the higher Shadow_Xlength also have higher pIC50 values. Num_StereoAtoms reflects that the fewer chiral atoms a ligand has, the higher the pIC50 value it possesses (for example, ligand 1). The observed and predicted pIC50 results and the values of physiochemical properties of the 59 ligands are listed in Table 2.

Table 2

Observed and predicted HDM2 inhibitory activities, physiochemical properties of different ligands from DS 2.1 used for the construction of QSAR models

Ligand No.LogDNum_RotatableBondsNum_StereoAtomsV_DIST_equCHI_1Dipole_XShadow_XlengthShadow_XZpIC50 (Obsa)pIC50(preda)Residual
15.968814724.7217.242.22616.104112.1693.0712.9950.076
26.346824975.4117.668-0.23217.246106.9273.1553.0990.056
36.624935948.6618.8798.91117.044105.9343.0092.8780.131
47.801103635719.35217.42617.251102.2092.9833.388-0.404
56.869825267.1218.2064.82716.575100.8923.1493.1400.009
66.819925586.7318.7060.52516.594112.0332.9552.960-0.005
74.35534029.7616.074-9.63614.33588.6272.6582.2060.452
82.993422898.3714.469-6.56913.22784.0851.421.444-0.024
93.367423165.4914.863-9.40313.42987.2551.8761.8370.039
103.823523466.2415.401-12.11613.10286.9752.1251.9670.158
114.602533750.5115.774-11.52513.01888.2931.7451.759-0.014
123.544423165.4914.863-4.63315.80686.2592.6022.2600.342
134.929634399.1716.548-4.98815.71485.6612.8792.8000.079
143.935533784.0616.091-3.58214.31884.7670.9031.175-0.272
153.618533903.0716.074-10.05513.47789.0931.3471.377-0.030
163.866433165.4914.863-7.03514.80188.0551.1941.538-0.344
174723034.5714.329-5.41217.16398.3721.8542.298-0.444
184.251733290.3414.684-2.95418.145106.531.9211.7650.156
193.909523473.3515.346-2.05613.72683.8341.9211.8720.049
204.659423442.815.257-11.19714.26290.4032.7992.7120.087
215.371534050.0716.168-3.29417.25899.5232.812.5380.272
225.191534342.9616.468-3.43118.061100.2422.6443.015-0.371
235.116534342.9616.468-9.90415.33890.3643.1742.9890.185
243.117423442.815.257-3.54516.2787.5751.7832.306-0.523
254.793423442.815.257-7.8115.66392.0373.3772.8720.505
264.333423442.815.257-5.8317.14991.6693.2082.9540.254
275.071534329.6216.468-8.79616.68397.3573.2082.9620.246
285.593634709.416.941-7.53217.3998.7833.1193.425-0.306
294.316423402.7815.257-7.62314.34490.5432.5692.3130.256
304.222423377.3415.274-2.81814.73285.3081.6992.216-0.517
313.531423204.9814.86311.82417.04294.7861.831.5940.236
324.869523713.5215.795-1.20616.84693.7992.8242.6120.212
335.121533982.1816.168-9.20316.69796.6582.9912.4950.496
345.327534252.116.468-8.90816.41695.6852.6442.845-0.201
354.919423442.815.257-8.05716.86794.1022.8153.142-0.327
364.419423442.815.257-5.74216.66189.3943.0132.9460.067
374.53423419.0915.274-7.34515.00789.472.2012.592-0.391
385.109423662.1515.684-8.06314.6791.6433.1552.7840.371
393.322523713.5215.795-11.79216.85694.8822.4812.4770.004
403.643624021.1616.295-12.17816.61296.0582.622.4900.130
413.127524319.0316.65120.78115.82196.5090.9031.020-0.117
424.171423713.5215.795-7.89216.85393.8531.832.822-0.992
435.008423713.5215.795-4.95516.94695.9863.062.9600.100
442.627423442.815.257-7.77515.55590.1932.8332.0650.768
453.847423442.815.257-7.71715.45590.2312.1312.512-0.381
463.901423897.8716.34616.63217.37898.0991.261.2590.001
474.649423944.9816.329-9.05415.08694.9992.1022.490-0.388
486.346824975.4117.66812.05115.95292.3663.0682.7010.367
494.1162522718.329-0.34515.15194.62.5692.2340.335
504.1162522718.3290.5315.82190.1382.2042.493-0.289
514.382625512.3118.74-0.68818.153106.9813.4352.8370.598
524.382625512.3118.74-0.86915.882107.3511.8832.312-0.429
535.581925015.0517.812-0.68317.066100.7922.6222.718-0.096
545.581925015.0517.812-1.10117.187113.3431.9032.354-0.451
554.803925294.318.222-0.45619.139102.3283.4053.0130.392
565.404726151.4519.7780.65618.773117.2983.1043.0820.022
573.694726364.8619.634-0.71516.25199.8473.2633.276-0.013
582.553525791.519.151-0.66617.56107.3722.2682.323-0.055
593.913525791.519.151-0.32115.98893.5122.812.929-0.119

The plot of the observed pIC50 vs. the predicted data is shown in Figure 6.

Plot of observed vs. predicted HDM2 inhibitory activities of different ligands in Table 1 with Equation 1

Plot of observed vs. predicted HDM2 inhibitory activities of different ligands in Table 1 with Equation 1

Figure 6

It can be seen that the predicted data by this model is basically in accordance with the experimental results. As a whole, it is only considered as a moderate QSAR model. In order to further improve the model quality, obtaining more descriptors is necessary. Thus, we collected 1620 kinds of molecular descriptors of BDPs using E-Dragon online tool. The QSAR model was obtained using GFA in MS 4.0. The QSAR model obtained is as follows:

pIC50 = 7.858 IDE - 4.855 MATS7v - 1.198 DP09 - 0.448 Mor14m + 1.481 Mor30p - 25.917 G2e + 1.678 E2e + 0.319 Tp + 61.575 R5u - 61.64 BELp6 + 0.272 SeaC2C3aa - 12.276

(Equation 2)

The sample number N = 59, LOF = 0.161, R2 = 0.882, 0.817 = R 2cv ,0.855 = R 2 adj , F = 31.98

The observed and predicted pIC50 results and the values of physiochemical properties of the 59 ligands for Equation 2 are listed in Table 3.

Table 3

Observed and predicted HDM2 inhibitory activities, physiochemical properties of different ligands from E-Dragon 1.0 used for the construction of QSAR models

Ligand No.IDEMATS7vDP09Mor14mMor30pG2eE2eTpR5uBELp6SeaC2C3aapIC50(Obsa)pIC50(preda)Residual
13.551-0.1079.4551.209-0.0320.1650.41518.7650.0311.478123.0712.5110.560
23.640.02710.3420.0290.0280.1420.4521.3910.0251.478123.1553.236-0.081
33.7070.01711.604-0.689-0.1170.1490.40422.480.031.478123.0092.8030.206
43.566-0.05910.281-0.0780.1540.150.28220.910.0321.515122.9832.9360.047
53.593-0.0910.5930.0330.1540.1670.37722.840.0331.555123.1493.0190.130
63.394-0.0278.432-0.047-0.1330.1510.51715.9050.0271.248102.9552.8060.150
73.394-0.0279.424-0.994-0.1020.1510.64717.590.0271.248102.6582.842-0.184
83.113-0.1567.7990.2-0.0120.1630.46415.1460.0291.22381.421.1450.275
93.191-0.1627.9920.486-0.1540.160.50615.6390.0261.232101.8761.8260.050
103.287-0.1668.950.0640.0150.1660.6118.7110.0321.389102.1252.271-0.146
113.35-0.179.2160.823-0.1420.1470.4919.9240.0271.412101.7452.119-0.374
123.191-0.1628.596-0.0680.1180.1730.66216.7910.0291.223102.6022.2870.315
133.504-0.0438.75-1.0980.0310.1680.43215.7640.0291.357102.8792.890-0.011
143.175-0.0628.281-0.2030.0150.1560.45716.0980.0241.30180.9030.4930.410
153.2810.0718.697-0.369-0.1570.1510.49617.120.031.276101.3471.594-0.247
163.191-0.1628.6110.4690.0130.1580.36917.5440.0271.23681.1941.260-0.066
173.27-0.1888.8660.7350.0470.1670.3719.0870.0371.29881.8542.118-0.264
183.326-0.1978.8120.8420.130.1470.36317.7670.0251.32581.9211.9190.002
193.294-0.1819.724-0.902-0.2310.1880.56719.3890.0291.224101.9211.967-0.046
203.263-0.1729.120.4050.1090.1690.32319.0860.0261.226122.7992.6550.144
213.412-0.1789.6760.6480.0740.1750.33221.8440.0241.402122.812.5360.274
223.457-0.1818.0380.767-0.0410.1610.35516.0940.0241.409122.6443.168-0.524
233.457-0.0429.973-0.727-0.0510.1510.46717.5330.0271.237123.1743.0020.172
243.263-0.0419.1130.2540.2580.1720.35218.1190.0261.223121.7831.997-0.214
253.263-0.1728.089-0.4360.070.1630.32615.730.0231.223123.3773.1340.243
263.263-0.2499.22-0.433-0.0880.1630.41318.6720.0241.223123.2083.0630.145
273.446-0.0429.62-0.840.0820.160.3518.4760.0261.248123.2083.326-0.118
283.549-0.0599.6340.031-0.0710.160.5217.2390.0311.357123.1193.0980.021
293.223-0.1628.958-0.5780.2070.1530.33317.80.0251.223122.5693.050-0.481
303.197-0.167.767-0.2560.0960.1830.35114.7960.0261.223101.6991.767-0.068
313.234-0.1858.014-0.2990.0680.1730.4614.4060.0281.253101.832.113-0.283
323.315-0.1728.436-0.2190.10.1580.30116.1520.0281.401122.8242.4840.340
333.352-0.1598.837-0.0780.1920.1560.31717.4050.0291.412122.9912.7750.216
343.381-0.1489.993-0.420.1220.1570.36521.1640.0261.416122.6442.657-0.013
353.263-0.1869.034-0.8380.1940.1530.30216.2220.0231.236122.8152.7270.088
363.263-0.1928.902-1.0760.0130.1630.49815.550.0291.231123.0133.0100.003
373.24-0.1958.036-0.6920.0350.1680.37214.1390.0311.252102.2012.398-0.197
383.266-0.1728.2970.3080.160.1580.35615.7620.0321.223103.1552.9080.247
393.315-0.28.874-0.6790.0980.150.33215.6850.031.367122.4812.746-0.265
403.3870.0228.081-0.8060.0270.1480.50915.3190.0221.41122.622.6050.015
413.4370.0029.481-0.6890.1970.1960.3117.9950.0281.423120.9031.181-0.278
423.315-0.29.572-0.5570.2180.1720.39216.8280.0241.316121.831.880-0.050
433.315-0.1729.301-1.0910.2310.1660.34217.8120.031.337123.062.9490.111
443.263-0.1218.1210.3920.2560.1740.3415.4880.0331.223112.8332.7580.075
453.263-0.1029.1580.138-0.0570.1530.39918.0890.0251.223112.1312.0520.079
463.277-0.1068.706-0.81-0.0620.160.41616.8390.0221.406111.261.2540.006
473.319-0.1878.964-0.322-0.140.1770.53817.7230.0261.367122.1022.144-0.042
483.551-0.10710.267-0.3590.1720.1550.50620.4930.0261.478123.0683.197-0.129
493.537-0.0839.3211.0530.0730.1420.41318.8660.0251.497122.5692.806-0.237
503.537-0.08310.8340.970.2220.1620.40823.1220.0291.497122.2042.327-0.123
513.543-0.03310.826-0.0340.2860.1480.40523.2390.0291.499123.4353.0680.367
523.543-0.03311.070.9380.4490.1660.37823.4060.0261.499121.8831.939-0.056
533.581-0.06110.6640.9380.2570.1630.41623.2930.0271.511122.6222.667-0.045
543.581-0.06110.6910.4930.0240.1430.39422.1180.021.511121.9032.164-0.261
553.585-0.0189.081-0.2110.0840.1490.45718.3930.021.512123.4053.0770.328
563.5810.0019.405-0.3760.1450.1390.50119.1130.0181.571123.1042.7980.306
573.711-0.0198.9680.7220.1360.150.44418.3730.0191.578123.2633.336-0.073
583.533-0.03510.559-0.2530.2230.1560.3721.5960.0251.494122.2682.319-0.051
593.540.0029.656-0.410.2180.1410.48919.5930.0251.507122.813.208-0.398
I3.577-0.1339.4380.2070.2530.1540.4619.5230.0221.617143.451
II3.565-0.069.2740.3280.3230.1470.41418.4970.0221.617143.025

Equation 2 can explain 85.5% of the variance (R2cv ) while it could predict 81.7% of the variance (R2cv ). F > F(a = 0.05) = 2.13 shows that the model is in the confidence interval of 95%. This model shows IDE, Mor30p, E2e, Tp, R5u and SeaC2C3aa. Count16 give positive contribution to the bioactivity of the ligands but MATS7v, Mor14m, G2e and BELp6 have the negative effect on the bioactivities of BDPs. The standardized regression coefficient for each variable is 1.772, -0.559, - 1.621, - 0.416, 0.300, -0.453, 0.217, 1.192, 0.336, -1.106 and 0.520 respectively. The relative importance of the descriptors according to their standardized regression coefficients is in the following order:

IDE > DP09 > Tp > BELp6 > MATS7v > SeaC2C3aa.Count16 > G2e > Mor14m > R5u > Mor30p > E2e

The plot of the observed pIC50 vs. the predicted data with Equation 2 is shown in Figure 7 and it illustrates that Equation 2 is a better model compared with Equation 1.

The plot of the observed pIC50 vs. the predicted data

The plot of the observed pIC50 vs. the predicted data

Figure 7

Prediction of some new HDM2 inhibitors

Based on the mentioned analysis, we designed some new compounds (Figure 8, I-III) and evaluated them with Equation 2. The predicted bioactivities are listed in Table 3. These compounds have higher inhibitory potency and their activities can be verified through chemosynthesis and FP assay lately. The docked pose of III using CDOCKER is shown in Figure 9.

From the structure of these coumpounds, we can see that the cyclopentane group gives the stability of the ligand conformation accormadating the cleft space. The extra aryl group at R3 improved the hydrophobic interaction between HDM2 and the ligands. It is to be clarified that this model is based on the molecule-based (FP assay) potency. Considering the cellular activity of these compounds, the water-solubility should be improved through introducing some polar groups. From this point, compound III may be more applicable as potent HDM2 inhibitors.

Conclusions

In order to explore appropriate binding mode of different 1,4-benzodiazepine-2,5-dione compounds as HDM2 inhibitors, CDOCKER protocol in DS 2.1 was employed to conduct the docking process using a dataset of 59 1,4-benzodiazepine-2,5-dione compounds from literatures (4, 15, 31, 32) and the binding situation was analyzed. Our docking pose of freshly prepared model of ligand 25 corroborates to the crystal structure, which indicates the reliability of this docking procedure. The docking results indicate that hydrophobic interaction between the imidazole group in HIS96 and the aryl ring in R3 may be one of the key factors for the combination of the ligands with HDM2. The binding cleft of HDM2 is predominantly hydrophobic and largely nonspecific Van Der Waals contacts are responsible for the interaction between the ligand and HDM2 hydrophobic pocket. Although the pocket of HDM2 is mainly hydrophobic, the hydrogen bonds can be formed with VAL93 residue in HDM2 for some ligands. The extra hydrogen bond formation can increase the inhibitory potency of the ligands to some extent. The chirality species of substituents in R1 to R4 influence the ligands’ potency seriously. Various 1D-3D physiochemical properties such as structural and Jurs parameters, and also topological descriptors from DS2.1 and E-Dragon 1.0 were selected as candidate descriptors to build two QSAR models separately. The QSAR model Equation 1 (with eight descriptors) built using more than 100 descriptors from DS2.1 was obtained. It can explain 71.2% of the variance (R2adj ) while it could predict 67.2% of the variance (R2adj ). F > F(a = 0.05) = 2.13 shows that the model is in the confidence interval of 95%. In order to improve the QSAR model quality, more than 1600 descriptors obtained from E-Dragon 1.0 were used to build the QSAR model with 11 descriptors (Equation 2). It can explain 85.5% of the variance (R2adj ) while it could predict 81.7% of the variance (R2cv). For both models, F > F(a = 0.05) = 2.13 shows that they are in the confidence interval of 95%. With these models, the bioactivities of some new compounds were predicted. The results show that the formation of cyclopentane ring at α-C or introduction of extra aryl group at original phenyl group of R3 can improve the ligands’ predicted pIC50 values, which may improve the HDM2 inhibitory potency.

Acknowledgements

References

  • 1.

    Vogelstein B, Kinzler KW. p53 Function and dysfunction. Cell. 1992;70:523-526. [PubMed ID: 1505019].

  • 2.

    Harris CC. Structure and function of the p53 tumor suppressor gene: clues for rational cancer therapeutic strategies. J. Natl. Cancer Inst. 1996;88:1442-1455. [PubMed ID: 8841019].

  • 3.

    Kubbutat MH, Vousden KH. Keeping an old friend under control: regulation of p53 stability. Mol. Med. Today. 1998;4:250-256. [PubMed ID: 9679243].

  • 4.

    Leonard K, Marugan JJ, Raboisson P, Calvo R, Gushue JM, Koblish HK, Lattanze J, Zhao S, Cummings MD, Player MR, Maroney AC, Lu T. Novel 1,4-benzodiazepine-2,5-diones as Hdm2 antagonists with improved cellular activity. Bioorg. Med Chem. Lett. 2006;16:3463-3468. [PubMed ID: 16647257].

  • 5.

    Chene P. Inhibiting the p53-MDM2 interaction: an important target for cancer therapy. Nat. Rev. Cancer. 2003;3:102-109. [PubMed ID: 12563309].

  • 6.

    Blaydes JP, Wynford-Thomas D. The proliferation of normal human fibroblasts is dependent upon negative regulation of p53 function by mdm2. Oncogene. 1998;16:3317-3322. [PubMed ID: 9681831].

  • 7.

    Chène P, Fuchs J, Bohn J, GarcIa-EcheverrIa C, Furet P, Fabbro D. A small synthetic peptide, which inhibits the p53- hdm2 interaction, stimulates the p53 pathway in tumour cell lines. J. Mol. Biol. 2000;299:245-253. [PubMed ID: 10860736].

  • 8.

    Bottger A, Bottger V, Sparks A, Liu WL, Howard SF, Lane DP. Design of a synthetic Mdm2-binding mini protein that activates the p53 response in-vivo. Curr. Biol. 1997;7:860-869. [PubMed ID: 9382809].

  • 9.

    Wasylyk C, Salvi R, Argentini M, Dureuil C, Delumeau I, Abecassis J, Debussche L, Wasylyk B. cells overexpressing MDM2 by an inhibitor of MDM2 interaction with p53. Oncogene. 1999;18:1921-1934. [PubMed ID: 10208414].

  • 10.

    Wang H, Nan L, Yu D, Agrawal S, Zhang R. oligonucleotides as a novel therapeutic approach to human breast cancer: In-vitro and in-vivo activities and mechanisms. Clin. Cancer Res. 2001;7:3613-3624. [PubMed ID: 11705884].

  • 11.

    Grasberger BL, Lu T, Schubert C, Parks DJ, Carver TE, Koblish HK, Cummings MD, LaFrance LV, Milkiewicz KL, Calvo RR, Maguire D, Lattanze J, Franks CF, Zhao S, Ramachandren K, Bylebyl GR, Zhang M, Manthey CL, Petrella EC, Pantoliano MW, Deckman IC, Spurlino JC, Maroney AC, Tomczuk BE, Molloy CJ, Bone RF. Discovery and cocrystal structure of benzodiazepinedione HDM2 antagonists that activate p53 in cells. J. Med. Chem. 2005;48:909-912. [PubMed ID: 15715460].

  • 12.

    Stoll R, Renner C, Hansen S, Palme S, Klein C, Belling A, Zeslawski W, Kamionka M, Rehm T, Mühlhahn P, Schumacher R, Hesse F, Kaluza B, Voelter W, Engh RA, Holak TA. Chalcone derivatives antagonize interactions between the human oncoprotein. MDM2 and p53. Biochem. 2000;40:336-344.

  • 13.

    Zhao J, Wang M, Chen J, Luo A, Wang X, Wu M, Yin D, Liu Z. The initial evaluation of non-peptidic small-molecule HDM2 inhibitors based on p53–HDM2 complex structure. Cancer Let. 2002;183:69-77. [PubMed ID: 12049816].

  • 14.

    Galatin PS, Abraham DJ. p53-dependent transcription in mdm2-overexpressing cells. J. Med. Chem. 2004;47:4163-4165. [PubMed ID: 15293988].

  • 15.

    Parks DJ, LaFrance LV, Calvo RR, Milkiewicz KL, Gupta V, Lattanze J, Ramachandren K, Carver TE, Petrella EC, Cummings MD, Maguire D, Grasberger BL, Lu T. 1,4-Benzodiazepine-2,5-diones as small molecule antagonists of the HDM2-p53 interaction: Discovery and SAR. Bioorg. Med. Chem. Lett. 2005;15:765-770. [PubMed ID: 15664854].

  • 16.

    Agrafiotis DK. Stochastic algorithms for maximizing molecular diversity. J. Chem. Inf. Comput. Sci. 1997;37:841-851.

  • 17.

    Rassokhin DN, Agrafiotis DK. Kolmogorov-Smirnov statistic and its application in library design. J. Mol. Graphics Model. 2000;18:368-382.

  • 18.

    Agrafiotis DK, Lobanov VS, Salemme FR. Combinatorial informatics in the post-genomics era. Nat. Rev. Drug Discov. 2002;1:337-346. [PubMed ID: 12120409].

  • 19.

    Pantoliano MW, Petrella EC, Kwasnoski JD, Lobanov VS, Myslik J, Graf E, Carver T, Asel E, Springer BA, Lane P, Salemme FR. High-density miniaturized thermal shift assays as a general strategy for drug discovery. J. Biomol. 2001;6:429-440.

  • 20.

    Cummings MD, Schubert C, Parks DJ, Calvo RR, LaFrance LV, Lattanze J, Milkiewicz KL, Lu T. Substituted 1,4-benzodiazepine-2,5-diones as alpha-helix mimetic antagonists of the HDM2-p53 protein-protein interaction. Chem. Biol. Drug Des. 2006;67:201-205. [PubMed ID: 16611213].

  • 21.

    Wang F, Li Y, Ma Z, Wang X, Wang Y. Structural determinants of benzodiazepinedione/peptide-based p53-HDM2 inhibitors using 3D-QSAR, docking and molecular dynamics. J. Mol. Model. 2011. DOI: 10.1007/s00894-011-1041-4.

  • 22.

    Galatin PS, Abraham DJ. The tumor suppressor p53: from structures to drug discovery. Proteins. 2001;45:169-175. [PubMed ID: 11599019].

  • 23.

    Xu H, Agrafiotis DK. Retrospect and prospect of virtual screening in drug lead discovery. Curr. Topics Med. Chem. 2002;2:1305-1320.

  • 24.

    Discovery Studio 2.1. San Diego: Accelrys Inc; 2007.

  • 25.

    MOE Molecular Operating Environment. Montreal: Chemical Computing Group Inc; 2008.

  • 26.

    SYBYL 8.0. Missouri: Tripos International; 2008.

  • 27.

    Tetko IV, Gasteiger J, Todeschini R, Mauri A, Livingstone D, Ertl P, Palyulin VA, Radchenko EV, Zefirov NS, Makarenko AS, Tanchuk VY, Prokopenko VV. Virtual computational chemistry laboratory - design and description. J. Comput. Aid Mol. Des. 2005;19:453-463.

  • 28.

    Roy P, Roy K. Docking and 3D-QSAR studies of diverse classes of human aromatase (CYP19) inhibitors. J. Mol. Model. 2010;16:1597-1616. [PubMed ID: 20195667].

  • 29.

    Dai Y, Wang Q, Zhang X, Jia S, Zheng H, Feng D, Yu P. Molecular docking and QSAR study on steroidal compounds as aromatase inhibitors. Eur. J. Med. Chem. 2010;45:5612-5620. [PubMed ID: 20926163].

  • 30.

    Wallace AC, Roman AL, Thornton Janet M. LIGPLOT: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995;8:127-134. [PubMed ID: 7630882].

  • 31.

    Parks DJ, LaFrance LV, Calvo RR, Milkiewicz KL, Marugan JJ, Raboisson P, Schubert C, Koblish HK, Zhao S, Franks CF, Lattanze J, Carver TE, Cummings MD, Maguire D, Grasberger BL, Maroney AC, Lu T. Enhanced pharmacokinetic properties of 1,4-benzodiazepine-2,5-dione antagonists of the HDM2-p53 protein–protein interaction through structure-based drug design. Bioorg. Med. Chem. Lett. 2006;16:3310-3314. [PubMed ID: 16600594].

  • 32.

    Marugan JJ, Leonard K, Raboisson P, Gushue JM, Calvo R, Koblish HK, Lattanze J, Zhao S, Cummings MD, Player MR, Schubert C, Maroney AC, Lu T. Enantiomerically pure 1,4-benzodiazepine-2,5-diones as Hdm2 antagonists. Bioorg. Med. Chem. Lett. 2006;16:3115-3120. [PubMed ID: 16630722].

  • 33.

    Wu G, Robertson DH, Brooks CL, Vieth M. Study of CDOCKER-a CHARMm-based MD docking algorithm. J. Comput. Chem. 2003;24:1549-1562. [PubMed ID: 12925999].

  • 34.

    ChemBioOffice Ultra Version 11.0. Cambridge: CambridgeSoft Inc; 2008.

  • 35.

    Kirkpatrick S, Gelatt CD, Vecchi MP. Optimization by simulated annealing. Science. 1983;220:671-680. [PubMed ID: 17813860].

  • 36.

    Kirkpatrick S. Optimization by Simulated Annealing: Quantitative Studies. J. Stat. Phy. 1984;34:975-986.

  • 37.

    So SS, Karplus M. Optimization in quantitative structure-activity. relationship: an application of genetic neural networks. J. Med. Chem. 1996;39:1521-1530. [PubMed ID: 8691483].

  • 38.

    Caballero J, Fernández L, Garriga M, Abreu JI, Collina S, Fernández M. Proteometric study of ghrelin receptor function variations upon mutations using amino acid sequence autocorrelation vectors and genetic algorithm-based least square support vector machines. J. Mol. Graphics Model. 2007;26:166-178.

  • 39.

    Materials Studio Modeling 4.0. San Diego: Accelrys Software Inc; 2007.

  • 40.

    Friedman JH. Multivariate adaptive regression splines. Ann. Statistics. 1991;19:1-67.

  • 41.

    Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera-a visualization system for exploratory research and analysis. J. Comput. Chem. 2004;25:1605-1612. [PubMed ID: 15264254].

  • 42.

    Dezi C, Carotti A, Magnani M, Baroni M, Padova A, Cruciani G, Macchiarulo A, Pellicciari R. Molecular interaction fields and 3d-qsar studies of p53−mdm2 inhibitors suggest additional features of ligand-target interaction. J. Chem. Inf. Model. 2010;50:1451-1465. [PubMed ID: 20726601].