Journal of Chemistry

Volume 2014 (2014), Article ID 323149, 10 pages

http://dx.doi.org/10.1155/2014/323149

Research Article

3D-QSAR and Docking Studies of a Series of β-Carboline Derivatives as Antitumor Agents of PLK1

Jahan B. Ghasemi and Valentin Davoudian

Chemistry Department, Faculty of Sciences, K. N. Toosi University of Technology, 16167 Tehran, Iran

Received 13 May 2013; Accepted 7 October 2013; Published 16 January 2014

Academic Editor: Deniz Ekinci

Copyright © 2014 Jahan B. Ghasemi and Valentin Davoudian. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An alignment-free, three dimensional quantitative structure-activity relationship (3D-QSAR) analysis has been performed on a series of -carboline derivatives as potent antitumor agents toward HepG2 human tumor cell lines. A highly descriptive and predictive 3D-QSAR model was obtained through the calculation of alignment-independent descriptors (GRIND descriptors) using ALMOND software. For a training set of 30 compounds, PLS analyses result in a three-component model which displays a squared correlation coefficient () of 0.957 and a standard deviation of the error of calculation (SDEC) of 0.116. Validation of this model was performed using leave-one-out, of 0.85, and leave-multiple-out. This model gives a remarkably high (0.66) for a test set of 10 compounds. Docking studies were performed to investigate the mode of interaction between -carboline derivatives and the active site of the most probable anticancer receptor, polo-like kinase protein.

1. Introduction

The -carboline alkaloids, which are originally isolated from the medicinal plant Peganum harmala [1–3], comprise a planar tricyclic system with different degrees of aromaticity and various substituents at positions 1, 2, 3, 7, and 9; the presence, location, and nature of these substituents play a crucial role in biological and pharmaceutical properties of these compounds [4–6]. According to previous investigations, -carbolines exert a broad spectrum of pharmacological effects such as anxiolytic, sedative [7–9], antimicrobial [10, 11], antiviral [12], antithrombotic [13, 14], and antiparasitic effects [15, 16]; also they are associated with alcohol dependence [17], and neurological diseases such as Parkinson’s disease [18, 19]. Moreover, a large series of -carbolines have been reported for their affinity with several receptors such as imidazoline [20, 21], benzodiazepine (BZ) [22–24], 5-hydroxytryptamine (5-HT) [25, 26], and dopamine (DA) [27].

Recently, -carboline alkaloids have drawn attention due to their antitumor activity [5, 28–33] and their functions through multiple mechanisms such as their ability to intercalate into DNA helix, inhibiting topoisomerases I and II [34, 35], monoamine oxidase [36–39], CDK (cyclin-dependent kinases) [40, 41], and MK-2 [42]. -Carbolines have the potential to be used as anticancer drug leads. One of the major antitumor mechanisms is apoptosis in cultured HepG2 cells induced by -carbolines with downregulate the expression of Bcl-2 gene and upregulate the expression of death receptor Fas,without altering the level of Bax and P53 [43].

Previous data suggest that -carbolines can inhibit the activity of polo-like kinases (PLKs). The investigations of Han et al. show that PLK1 is attractive as a therapeutic target for inhibition of cancer cell growth [44, 45]. Among the four members of PLKs which are a family of serine-threonine kinases, PLK1 plays the most crucial role in cell proliferation and is an important regulatory enzyme in the G2/M transition through phosphorylation of substrates [46]. PLK1 contains an N-terminal catalytic domain and a C-terminal regulatory polo box domain (PBD) which is composed of two polo boxes. Since the PBD is contributing to PLK1 localization and function, the PBDs of PLK1 are promising alternative targets for designing novel antitumor drugs [47, 48].

The docking study was done using CDOCKER algorithm and the picture of the virtual receptor site was validated by investigation interactions between receptor and some inhibitors. CDOCKER is an implementation of a CHARMm-based docking tool using a rigid receptor that generates several prime random ligand orientations within the receptor active site followed by MD-based simulated annealing and final refinement by minimization [49].

This paper intended to get a better pharmacophoric pattern of these compounds and thus we obtained an alignment-independent 3D-QSAR model for the potency of 40 from total 47 -carboline derivatives which have significant pIC₅₀ toward HepG2 human tumor cell lines (Table 1).

Table 1: Structures and potencies toward HepG2 tumor cell lines of compounds 1–40.

Comp.
Substituents
pIC₅₀ (exp.)
pIC₅₀ (cal.)

R₁
R₂
R₃
R₇
R₉

1
3,4,5-trimethoxyphenyl
—
CO₂C₂H₅
H
H
3.78
3.97

3
3,4,5-trimethoxyphenyl
—
CO₂H
H
n-C₄H₉
3.68
3.64

4
3,4,5-trimethoxyphenyl
—
CONH(CH₂)₂OH
H
n-C₄H₉
4.06
3.92

11
H
—
CONH(CH₂)₂NH₂
H
n-C₄H₉
4.08
4.24

12
H
—
CONH(CH₂)₆NH₂
H
n-C₄H₉
4.02
3.96

CH₃
—
CONH(CH₂)₂OH
H
n-C₄H₉
3.85
4.15

14
CH₃
—
CONH(CH₂)₂NH₂
H
n-C₄H₉
4.66
4.68

CH₃
—
CONH(CH₂)₂NH₂
H
CH₂C₆H₅
4.46
3.77

17
H
—
CONH(CH₂)₂NH₂
H
CH₂C₆H₅
3.98
3.9

H
—
CONH(CH₂)₆NH₂
H
CH₂C₆H₅
4.23
4.27

CH₃
—
CH₂OH
H
n-C₄H₉
3.89
3.74

23
CH₃
—
CHO
H
n-C₄H₉
3.84
3.69

29
CH₃
CH₂C₆H₅
H
H
H
4.16
4.23

30
CH₃
(CH₂)₃C₆H₅
H
H
H
4.48
4.56

31
CH₃
CH₂C₆H₅
CO₂C₂H₅
H
H
4.28
4.22

32
CH₃
CH₂C₆H₅
H
OCH₃
H
4.26
4.08

33
H
n-C₄H₉
H
H
H
4.03
4.04

H
CH₂C₆H₅
H
H
H
4.11
4.22

35
H
(CH₂)₃C₆H₅
H
H
H
4.35
4.42

40
CH₃
—
H
OH
C₂H₅
3.87
3.89

CH₃
—
H
OH
n-C₄H₉
4.09
4.42

42
CH₃
—
H
OH
i-C₄H₉
3.94
4.03

43
CH₃
—
H
OH
(CH₂)₃C₆H₅
4.55
4.52

44
CH₃
—
H
OC₂H₅
C₂H₅
4.16
4.21

45
CH₃
—
H
OCH₂C₆F₅
C₂H₅
3.63
3.71

CH₃
—
H
OC₂H₅
n-C₄H₉
4.36
4.38

47
CH₃
—
H
OCH(CH₃)₂
n-C₄H₉
4.52
4.44

48
CH₃
—
H
OC₄H₉
n-C₄H₉
4.81
4.74

49
CH₃
—
H
OC₁₀H₂₁
n-C₄H₉
3.90
3.71

50
CH₃
—
H
OC₄H₉
i-C₄H₉
4.92
4.94

51
CH₃
—
H
OCH₂C₆H₅
i-C₄H₉
4.65
4.66

CH₃
—
H
OCH(CH₃)₂
(CH₂)₃C₆H₅
4.84
4.5

53
CH₃
—
H
OC₈H₁₇
(CH₂)₃C₆H₅
3.98
4.12

CH₃
—
H
OCH₂C₆H₅
(CH₂)₃C₆H₅
4.80
4.44

55
CH₃
—
H
OCH₂C₆F₅
(CH₂)₃C₆H₅
3.83
4.05

56
CH₃
CH₂C₆H₅
H
OC₂H₅
C₂H₅
4.84
4.76

57
CH₃
CH₂C₆H₅
H
OCH₂C₆F₅
C₂H₅
5.80
5.65

58
CH₃
CH₂C₆H₅
H
OC₄H₉
i-C₄H₉
5.74
5.6

CH₃
CH₂C₆H₅
H
OCH₂C₆F₅
i-C₄H₉
5.72
5.33

60
CH₃
CH₂C₆H₅
H
OC₈H₁₇
(CH₂)₃C₆H₅
5.41
5.63

Test set.

2. Materials and Methods

2.1. GRIND Descriptors

A new class of molecular descriptors called GRid-INdependent descriptors has been developed by Pastor et al. [50]. There are many advantages in creating such a predictive and robust 3D-QSAR model with a mathematical description of molecules over the other methods. First, in contrast to other methodologies in which one of the most difficult and time-consuming steps is molecular superimposition, in GRIND there is no need for alignment of compounds. Therefore the descriptors obtained are insensitive to coordinate frame of the space. Moreover, the GRIND descriptors are easily interpretive by referring back to the compounds so the original information can be reconstructed. Finally because of smaller variables/objects ratio in this method in comparison with other 3D-QSAR methods, the models obtained by improper use of variable selection with FFD (fractional factorial design) variable selection procedure are less prone to overfitting.

With the aim of getting a virtual molecular interaction fields (MIFs) pattern of the regions (fingerprint of receptor) so-called virtual receptor sites (VRS) to reveal the most common structural groups in the active compounds, first we should compute the MIFs of nodes. Therefore for the derivation of MIFs we used the four most recommended probes. To represent steric and hydrophobic interaction, hydrogen bond acceptor groups, and hydrogen bond donor groups, we use DRY (hydrophobic probe), O (carbonyl oxygen), and N1 (amide nitrogen), respectively. These probes represent strong noncovalent interactions between ligands and receptor. In addition, to consider molecular shape effects in the receptor-ligand interaction process, and as complementary to point-based interaction information, we used a supplementary probe, called TIP (shape probe), that extracts each ligand’s isosurface at 1 kcal/mol from the field of a normal GRID calculation. This method has been described elsewhere in more details [50, 51], but it is worth mentioning the whole procedure briefly.

In the absence of binding site knowledge it is possible to explore the process of ligand-receptor interactions with the help of the MIFs [52, 53]. The procedure for obtaining the GRIND starts with calculation of MIFs of the ligands using the program GRID, in order to obtain VRS described above. Due to the vast number of nodes, between 10000 to 100000, from which we obtain the MIFs for a drug molecule, we should filter these MIFs to obtain the most relevant regions and remove redundant information. This filtering uses an optimization algorithm, called Federov-like optimization algorithm [54], that selects from each field a fixed number of nodes (as a default value of 100) with optimizing a scoring function. MIF points that have low energy value and are as far as possible are retained. Also the scoring function can be tuned by giving to each criterion a different relative weight. In the next step the geometrical relationship between the VRS regions should be encoded into GRIND variables; thus there is no dependence upon the position and orientation of the molecular structure in the 3D space. The encoding is based on auto- and cross-correlation transform method called maximum auto- and cross-correlation (MACC) or, more specifically, MACC2. The continuous distances between pairs of filtered MIF points are divided into discrete “distance bins.” For each distance bin the product of the interaction is calculated but only the maximum energy product is conserved while the others are discarded. The output of MACC2 analysis can be represented directly in correlograms, where each point in it represents the product of two particular nodes within the distance bin separating the nodes of a certain compound. Wide peaks represent intense interactions which are produced by contiguous nodes around the compound. Conversely, weaker interactions are represented by narrow peaks. Each block in auto- and cross-correlogram indicates interactions between couples of nodes (node pairs) generated by the same or different probes, respectively. This approach does not require any scaling or pretreatment. The default settings were used for all other parameters.

2.2. Dataset

The dataset was adopted from the work of Cao et al. [55] in which they reported the cytotoxic potential of a number of synthesized -carbolines against a panel of human tumor cell lines (Table 1). The 3D structures of dataset were constructed in SYBYL7.3 molecular modeling package (Tripos Inc., St. Louis, USA), and the energy minimizations were performed using Tripos force field with a distance-dependent dielectric and the Powell conjugate gradient algorithm convergence criterion of 0.01 kcal/mol Å. Gasteiger-Hückel method is used to calculate partial atomic charge of all compounds. All other computational procedures were performed within the software ALMOND (version 3.3.0) running on Red Hat Enterprise Linux 4.7 workstation.

The dataset was partitioned into training and test set using most descriptive compound (MDC) method in which the compounds are weighted according to their population density [56]. The training set is composed of 30 compounds which were used to adjust the parameters of the models. The corresponding pIC₅₀ values which are listed in Table 1 range from 3.63 for the most weak compound 45 to a value of 5.8 for the most potent compound 57 and cover a spectrum of approximately 3-log units. The 3D-QSAR model derived was successfully validated by using a test set of 10 similar compounds, with comparable pIC₅₀.

2.3. ALMOND Model

The -matrix was generated, with compounds as its rows and the GRIND variables grouping into correlogram’s blocks as its columns. The information contained in this matrix can be directly used for chemometric analysis such as principal component analysis (PCA) and partial least square (PLS) analysis. A principal component analysis (PCA) was carried out to make a closer inspection on the distributions of the variables and compounds. The optimum number of PLS components (latent variables, LV) was chosen by monitoring changes in the model’s predicting index ( and leave-one-out) evaluated by applying the cross-validation procedure available in ALMOND. The optimum number was given at three PLS components. In order to select the most informative variables in the -matrix, FFD is applied that yields 145 -variables out of total 437 original -variables. The same data matrix was subjected to PLS analysis. Statistical parameters of the model obtained are summarized in Table 2. PLS analysis resulted in a 3-component model; a correlation coefficient of 0.96 and a standard deviation of the error of calculation (SDEC) of 0.12 were found. Also, the validation of this model using cross-validation method and an external test set shows = 0.85 of the validated model and the for the external test set was 0.66. The plot of experimental versus calculated values of pIC₅₀ in Figure 1 proves the good quality of the PLS model obtained.

Table 2: ALMOND-PLS results of the model developed for the series of -carboline alkaloids.

LV

SDEC
SDEP

1
0.63
0.46
0.45
0.34
0.41

2
0.94
0.81
0.77
0.14
0.24

3
0.96
0.85
0.82
0.12
0.22

4
0.97
0.80
0.81
0.09
0.23

Figure 1: Plot of predicted versus experimental pIC₅₀ values for ALMOND model.

2.4. Molecular Docking

Molecular docking studies were carried out by using CDOCKER (CHARMm-based DOCKER) [57]. The crystal structure of PLK1 was retrieved from the RCSB Protein Data Bank (PDB entry code: 2OWB). For the preparation step of ligands all structures sketched and modified in SYBYL 7.3 molecular modeling package were transferred into Discovery Studio 2.5 (Accelrys Inc, San Diego, CA, USA) and typed with CHARMm force field and Momany-Rone partial charges calculation method [58]. Energy minimization was performed by SMART minimization algorithm, followed by conjugate gradient minimization [49]. The active site of the target protein is created as a spherical region with a radius of 9.18 Å. The next step is protein preparation in which all ligands and water molecules in PLK1(2OWB) were deleted and hydrogen atoms were added to the original crystal structure.

The CDOCKER score (-CDOCKER ENERGY) as a negative value includes receptor-ligand interaction energy and internal ligand strain energy [58]. Default values were assigned to the other parameters. High value of 25.45 for compound 57 as its dock score indicates favorable binding between this active inhibitor and kinase domain of PLK1.

3. Results and Discussion

3.1. The Results of ALMOND Model

PLS coefficients for the model are calculated using the DRY, O, N1, and TIP probes [59]. PLS coefficient plot representing the contribution of each single descriptor to the model with respect to the value of is shown in Figure 2. Since the regions of the MIF corresponding to each important variable can be visualized with the 3D plots included in ALMOND, the interpretation of these variables is straightforward. From PLS coefficient plot it can be seen that the most important variables which have positive impact on inhibitory activity are DRY-DRY: 24 and 25, TIP-TIP: 52 and 56, N1-TIP: 23 and 24, and N1-N1: 28. On the contrary, the analysis of all the distances at higher PLS coefficient showed that the variables TIP-TIP: 42, 43, and 44, N1-TIP: 10 and 7, and N1-N1: 8 correlate negatively with the activity. According to this plot the probe O has a marginal contribution.

Figure 2: PLS coefficient plots for 3D-QSAR model. Direct and reverse correlation with the activity are indicated with positive and negative PLS coefficients, respectively. Bars with the most intensive height in the PLS plots have the most profound impact on the model obtained.

3.2. ALMOND and Docking Interpretation

Variables 11–24 and 11–25 which have the strongest positive impact on antitumor potency are within the block of DRY-DRY node pairs. The variable 24 represents a distance of 9.6 Å between hydrophobic nodes due to ring A and phenyl ring in 2-benzyl substituent in the two compounds 57 and 58 (Figure 3(b)) (the bold letter faces show the number of compounds). The variable 25 is the same for compound 57 but in a distance of 10 Å in compound 58 this variable is due to hydrophobic properties of ring B and 2-benzyl in this series of -carboline derivatives.

(a)

(b)

Figure 3: Association of structural fragments with GRIND variables. Nodes were selected using ALMOND filtering for DRY and TIP probes. Distances in red are favorable whereas distances in blue are not favorable for antitumor potency.

According to Figure 2 all DRY-DRY interactions are positively related to antitumor potency and differ in their relative intensity. These variables referring to favorable regions for hydrophobic groups mainly arise from interactions of triple ring system of all -carbolines and the aromatic moiety far from template. However the long aliphatic saturated chains are hydrophobic regions in many compounds (e.g., the compounds 33, 47, 48, 50, and 51). These variables can be seen clearly in the active compounds such as the compounds 56, 57, 58, and 60. Those variables are also present in compound 45, due to its pentafluorobenzyl moiety, although it is one of the most inactive compounds.

The results of molecular docking -carboline derivatives of PLK1 are overlaid in Figure 4; it is clearly observed that compound 57, which is one of the most active compounds, interacts with the active site of PLK1 through an extensive hydrophobic region, including mainly the following residues: PHE58, LEU59, GLY60, ALA80, VAL114, LEU130, LEU132, PHE183, and GLY193; thus hydrophobic groups would benefit the potency.

Figure 4: Docked binding mode of compound 57 into active site of the PLK1 (PDB code 2WOB). The protein key residues that form the main interactions with the different structural units of the inhibitor are labeled. Hydrogen bonds are represented as red dots. The stacking interaction between PHE183 and compound 57 is shown as an orange line.

Compound 57 also exhibits a strong binding affinity through stacking with the PHE183 residue of PLK1. The stacking interaction between PHE183 of PLK1 and aromatic rings either in the backbone or one of substitutions of compounds like pentafluorobenzoxyl moiety has crucial role in enhancing the binding affinity since it is absent from most of the compounds with lower pIC₅₀.

Within the block of TIP-TIP node pairs, the largest positive impact on antitumor potency is attributed to variables 52 and 56 (node-pair distances of 20.8 and 22.4 Å, resp.). The variables 44–52 and 44–56 are due to the terminal benzyl ring and the aromatic or ether group on the other end of the molecule in compounds 56, 57, 58, and 60. Compound 50 is active but since it has no benzyl group this variable is zero according to the correlogram. Actually these variables are the most discriminating variables that can differentiate active from inactive compounds (Figure 5). On the other hand, GRIND variables 42, 43, and 44 that correspond to TIP-TIP node-pair distances of 16.8, 17.2, and 17.6 Å, respectively, exhibit strong inverse correlation to the biological activity (). So these variables are present in compounds 3 and 45 (Figure 3(a)).

Figure 5: Visual inspection of variables to identify the most discriminative ones. Red and blue boxes include active and inactive compounds, respectively. The numbers written in red represent the variables which have a high positive impact on antitumor activity and the variables with a significant negative impact are written in blue. If a variable is present (✓) or absent (—) in a compound as expected, the relative box is highlighted in green; otherwise it is highlighted in yellow.

Within the N1-N1 block in PLS coefficient plot, it is worth mentioning that the negatively correlated bars are located on the left side (representing smaller node-node distances) and the positively correlated variables are positioned on larger node-node distances, that is, on the right side. It is due to the fact that structural elements which exert negative impact on potency are located closer to each other in molecules than those having positive impact. According to the N1-N1 autocorrelogram, regions favorable for hydrogen bond acceptor groups separated by 3.2 Å (i.e., variable 33–8) seem to be characteristic of inactive compounds. They are associated to compounds 4, 23, 3, 17, 11, 55, and 45 in respect to their energy products value. These variables correspond to the interactions of the polar carbonyl and carboxyl groups as hydrogen bond acceptors separated by 3.2 Å, being of similar energy product for the two compounds 3 and 23.

Figure 6 shows that the carboxyl group at position 3 of compound 3 binds via a hydrogen bond to amino group of ARG136 which indicates that the presence of H-bond donor groups in this region would emphasize the lower affinity of this compound to the active site of enzyme.

Figure 6: Molecular docking between compound 3 and PLK1, in which key residues of protein are shown as sticks. The dotted red line indicates hydrogen bonding interaction.

Also there is a similarity in node-node energy products between compounds 11 and 17 due to primary amine nitrogen, and also between compounds 45 and 55 due to pentafluoro benzoxyl ring. The highest value of energy interaction product for the variable 33–8 in compound 4 which is clearly evident in correlogram means that the N1 probe interacts more strongly with nitrogen atom in –C(O)NHCH₂CH₂OH than the groups described above.

According to the docking results there are two hydrogen bonds in the –C(O)NHCH₂CH₂OH moiety. One is between oxygen of CONH and alcoholic hydrogen within the molecule itself, and the other is between alcoholic oxygen at the end of the moiety and LYS61 amino acid.

Overall, the hydrogen bond acceptor group in the side chain of ring A at position 7 forms a hydrogen bond with LYS82 or ASP194, such as oxygen of ether group in compounds 56, 32, and 43 or fluorine of pentafluorobenzoxyl in compound 57, which plays an important role in the PLK1 inhibitory activity of -carboline. By contrast, presence of hydrogen donor in this position forms unfavorable hydrogen bond with LYS61 like in compounds 40, 42, and 43. The next hydrogen bond formation represents negative influence on inhibitory potency observed for compounds 49 and 53 due to pyridinium nitrogen (ring C) so this nitrogen atom in its protonated form would be more desirable for potent inhibitory activity.

The interpretation of all relevant peaks is summarized in Table 3.

Table 3: The variables with the highest impact (positive or negative) in the final GRIND model and the structural elements of the most active (57 and 58) and inactive (3 and 45) compounds associated with these variables.

Node pair
GRIND
Interpretation

Variable no.

(Å)
Impact (coefficient)

TIP-TIP
42, 43

(16.8–17.2)
Inverse

(−0.23, −0.24)
Shape of ring A and one of the methoxy substituents in 3,4,5-trimethoxyphenyl at position 1 for 3 and shape of pentafluorobenzoxyl at position 7 and hydrogens of ring C for 45

TIP-TIP
52

(20.8)
Direct

(0.21)
Shape of pentafluorobenzoxyl at position 7 and N²-arylated substituent for 57 and shape of 7-n-butoxy and 2-benzyl for 58

DRY-DRY
24

(9.6)
Direct

(0.2)
Hydrophobic properties of ring A and phenyl ring in 2-benzyl substituent for 57 and 58

N1-N1
8

(3.2)
Inverse

(−0.18)
Interaction of carboxyl at position 3 with the probe N1 for 3 and interaction of oxygen atom in pentafluorobenzoxyl group at position 7 with probe N1 for 45

TIP-TIP
44

(17.6)
Inverse

(−0.18)
Shape of ring A and one of the methoxy substituents in 3,4,5-trimethoxyphenyl at position 1 for 3 and shape of pentafluorobenzoxyl at position 7 and hydrogens of ring C for 45

N1-TIP
10

(4)
Inverse

(−0.18)
Shape of 3-carboxyl and the interaction of ring C with probe N1 for 3 and interaction of hydrogens of ring A with probe N1 and the shape of ethyl substituent in ring B (TIP) for 45

DRY-DRY
25

(10)
Direct

(0.17)
Hydrophobic properties of ring A and phenyl ring in 2-benzyl substituent for 57 and hydrophobic properties of ring B and 2-benzyl for 58

N1-TIP
24

(9.6)
Direct

(0.15)
Shape of 2-benzyl ring (TIP) and nitrogen of ring C (N1) for 57 and shape of i-butyl substituent in ring B and interaction of nitrogen in ring C with probe N1 for 58

N1-TIP
7

(2.8)
Inverse

(−0.15)
Shape of 3-carboxyl and the interaction of ring C with probe N1 for 3. This variable is absent from compound 45

TIP-TIP
56

(22.4)
Direct

(0.13)
Shape of pentafluorobenzoxyl at position 7 and N²-arylated substituent for 57. This variable is absent from compound 58

N1-TIP
23

(9.2)
Direct

(0.13)
Shape of N²-benzyl and the interaction of oxygen in pentafluorobenzoxyl group with probe N1 for 57 and shape of 7-n-butoxy and the interaction of N in ring C with the probe N1 for 58

N1-N1
28

(11.2)
Direct

(0.1)
Interaction of nitrogen in ring C and O in pentafluorobenzoxyl group with probe N1 for 57 and interaction of nitrogen in ring C and O in 7-n-butoxy group with probe N1 for 58

At a deeper insight, in order to determine the most relevant descriptors which can properly distinguish between active and inactive compounds, 6 compounds of training set (1, 3, 45, 57, 58, and 60) and 2 compounds of test set (21 and 59) were selected; then existence of all variables which have high positive or negative PLS coefficients described above is carefully inspected for these 8 compounds, including 4 active and 4 inactive -carboline derivatives (Figure 5). Variables with positive PLS coefficient should be present in active compounds and should be absent from inactive compounds and vice versa. Green boxes in Figure 5 represent presence or absence of variables as expected, and yellow boxes represent unexpected results. As it is shown, variables 44–52, 33–8, and 34–10 are the most discriminative variables. Since more than 65% of boxes are highlighted in green it can be concluded that the obtained model is very predictive and efficient. Also the yellow cells are a consequence of attempting a univariate interpretation of a multivariate model.

3.3. Summary of Structural Insights from 3D-QSAR Models and Docking Studies

The final results of analysis of 3D-QSAR models and docking studies are summarized in Figure 7. As discussed above, the presence, spatial shape and chemical characteristics of R2 and R7, and their mutual distance have high impact on -carbolines inhibitory activity. Pyridinium nitrogen of ring C can form hydrogen bond with ASP194 and LYS82; however there is negative association between presence of this moiety and the activity. Also there should be an electron donating group like benzyl to form N1-N1:28. This substituent should be bulky enough to form the TIP-TIP variable in a distance of 20.8–22.4 Å from R7. Otherwise it may represent the TIP-TIP variable in a distance of 16.8–17.6 Å, with negative impact. At position 3, carbonyl or carboxyl groups as hydrogen bond acceptor are disfavored, because they form the N1-N1: 8 variable with negative impact.

(a)

(b)

Figure 7: Schematic representation of the results revealed by the present study. (a) As variables of 3D-QSAR model. Yellow, green, and blue circles represent hydrophobic (DRY), shape (TIP), and hydrogen bond acceptor (N1) regions, respectively. Red and blue lines represent the variables with positive and negative impact, respectively. (b) As structural requirements of -carboline derivatives as PLK1 inhibitors.

4. Conclusion

Within this paper we aimed at development and validation of a ligand-based 3D-QSAR model in order to get a deeper insight into molecular structure-antitumor potency relationship in a series of -carboline derivatives. The method used in this paper is based on alignment independent descriptors whose particular value is the ability of backtracking of all interactions to the original filtered GRID field.

The hydrophobicity (associated with the DRY probe), shape effects (associated with TIP probe), and hydrogen bond acceptor-donor interactions (associated with N1 probe) are the main factors that determine antitumor potency toward HepG2 cell lines, within studied set. In addition, molecular docking is carried out to map the binding pocket of the PLK1 and its characteristics. And the interactions described with docking studies are similar to those described through analysis of 3D-QSAR model. Results obtained in the present study can be used as guidance for design of novel drugs.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

J. Cheng and K. R. Mitchelson, “Improved separation of six harmane alkaloids by high-performance capillary electrophoresis,” Journal of Chromatography A, vol. 761, no. 1-2, pp. 297–305, 1997.

M. Kartal, M. L. Altun, and S. Kurucu, “HPLC method for the analysis of harmol, harmalol, harmine and harmaline in the seeds of Peganum harmala L.,” Journal of Pharmaceutical and Biomedical Analysis, vol. 31, no. 2, pp. 263–269, 2003.

Q. Chen, R. Chao, H. Chen et al., “Antitumor and neurotoxic effects of novel harmine derivatives and structure-activity relationship analysis,” International Journal of Cancer, vol. 114, no. 5, pp. 675–682, 2005.

R. Cao, H. Chen, W. Peng et al., “Design, synthesis and in vitro and in vivo antitumor activities of novel β-carboline derivatives,” The European Journal of Medicinal Chemistry, vol. 40, no. 10, pp. 991–1001, 2005.

R. Cao, W. Yi, Q. Wu et al., “Synthesis and cytotoxic activities of 1-benzylidine substituted β-carboline derivatives,” Bioorganic and Medicinal Chemistry Letters, vol. 18, no. 24, pp. 6558–6561, 2008.

R. Cao, W. Peng, Z. Wang, and A. Xu, “β-carboline alkaloids: biochemical and pharmacological functions,” Current Medicinal Chemistry, vol. 14, no. 4, pp. 479–500, 2007.

C. Braestrup, M. Nielsen, and C. E. Olsen, “Urinary and brain β-carboline-3-carboxylates as potent inhibitors of brain benzodiazepine receptors,” Proceedings of the National Academy of Sciences of the United States of America, vol. 77, no. 4, pp. 2288–2292, 1980.

W. Schlecker, A. Huth, E. Ottow, and J. Mulzer, “Regioselective metalation of 9-methoxymethyl-β-carboline-3-carboxamides with amidomagnesium chlorides,” Synthesis, vol. 131, no. 10, pp. 1225–1227, 1995.

A. Batch and R. H. Dodd, “Ortho-directed metalation of 3-carboxy-β-carbolines: use of the SmI₂-cleavable 9-N-(N′,N′-dimethylsulfamoyl) blocking group for the preparation of 9-N-deprotected 4-amino derivatives via azide introduction or a palladium-catalyzed cross-coupling reaction,” Journal of Organic Chemistry, vol. 63, no. 3, pp. 872–877, 1998.

P. Schupp, T. Poehner, R. A. Edrada et al., “Eudistomins W and X, two new β-carbolines from the micronesian tunicate Eudistoma sp.,” Journal of Natural Products, vol. 66, no. 2, pp. 272–275, 2003.

Y. Iinuma, S. Kozawa, H. Ishiyama et al., “Gesashidine A, a β-carboline alkaloid with an imidazole ring from a thorectidae sponge,” Journal of Natural Products, vol. 68, no. 7, pp. 1109–1110, 2005.

A. S. Nazari Formagio, P. R. Santos, K. Zanoli et al., “Synthesis and antiviral activity of β-carboline derivatives bearing a substituted carbohydrazide at C-3 against poliovirus and herpes simplex virus (HSV-1),” The European Journal of Medicinal Chemistry, vol. 44, no. 11, pp. 4695–4701, 2009.

N. Lin, M. Zhao, C. Wang, and S. Peng, “Synthesis and antithrombotic activity of carbolinecarboxyl RGD sequence,” Bioorganic and Medicinal Chemistry Letters, vol. 12, no. 4, pp. 585–587, 2002.

M. Zhao, L. Bi, W. Bi et al., “Synthesis of new class dipeptide analogues with improved permeability and antithrombotic activity,” Bioorganic and Medicinal Chemistry, vol. 14, no. 14, pp. 4761–4774, 2006.

C. di Giorgio, F. Delmas, E. Ollivier, R. Elias, G. Balansard, and P. Timon-David, “In vitro activity of the β-carboline alkaloids harmane, harmine, and harmaline toward parasites of the species Leishmania infantum,” Experimental Parasitology, vol. 106, no. 3-4, pp. 67–74, 2004.

S. Nafisi, M. Bonsaii, P. Maali, M. A. Khalilzadeh, and F. Manouchehri, “β-carboline alkaloids bind DNA,” Journal of Photochemistry and Photobiology B, vol. 100, no. 2, pp. 84–91, 2010.

C. Melchior and M. A. Collins, “The route and significance of endogenous synthesis of alkaloids in animals,” Critical Reviews in Toxicology, vol. 9, no. 4, pp. 313–356, 1982.

M. A. Collins, “Alkaloids, alcohol and Parkinson's disease,” Parkinsonism and Related Disorders, vol. 8, no. 6, pp. 417–422, 2002.

D. A. Gearhart, P. F. Toole, and J. W. Beach, “Identification of brain proteins that interact with 2-methylnorharman: an analog of the parkinsonian-inducing toxin, MPP+,” Neuroscience Research, vol. 44, no. 3, pp. 255–265, 2002.

S. M. Husbands, R. A. Glennon, S. Gorgerat et al., “β-carboline binding to imidazoline receptors,” Drug and Alcohol Dependence, vol. 64, no. 2, pp. 203–208, 2001.

R. A. Glennon, B. Grella, R. J. Tyacke, A. Lau, J. Westaway, and A. L. Hudson, “Binding of β-carbolines at imidazoline I₂ receptors: a structure-affinity investigation,” Bioorganic and Medicinal Chemistry Letters, vol. 14, no. 4, pp. 999–1002, 2004.

M. Cain, R. W. Weber, F. Guzman et al., “β-carbolines: synthesis and neurochemical and pharmacological actions on brain benzodiazepine receptors,” Journal of Medicinal Chemistry, vol. 25, no. 9, pp. 1081–1091, 1982.

S. P. Hollinshead, M. L. Trudell, P. Skolnick, and J. M. Cook, “Structural requirements for agonist actions at the benzodiazepine receptor: studies with analogues of 6-(benzyloxy)-4-(methoxymethyl)-β-carboline-3-carboxylic acid ethyl ester,” Journal of Medicinal Chemistry, vol. 33, no. 3, pp. 1062–1069, 1990.

R. H. Dodd, C. Ouannès, L. Prado de Carvalho et al., “3-amino-β-carboline derivatives and the benzodiazepine receptor. Synthesis of a selective antagonist of the sedative action of diazepam,” Journal of Medicinal Chemistry, vol. 28, no. 6, pp. 824–828, 1985.

B. Grella, M. Teitler, C. Smith, K. Herrick-Davis, and R. A. Glennon, “Binding of β-carbolines at 5-HT2 serotonin receptors,” Bioorganic and Medicinal Chemistry Letters, vol. 13, no. 24, pp. 4421–4425, 2003.

R. A. Glennon, M. Dukat, B. Grella et al., “Binding of β-carbolines and related agents at serotonin (5-HT2 and 5- HT(1A)), dopamine (D2) and benzodiazepine receptors,” Drug and Alcohol Dependence, vol. 60, no. 2, pp. 121–132, 2000.

A. F. M. Abdel-Fattah, K. Matsumoto, H. A. K. Gammaz, and H. Watanabe, “Hypothermic effect of harmala alkaloid in rats: involvement of serotonergic mechanism,” Pharmacology Biochemistry and Behavior, vol. 52, no. 2, pp. 421–426, 1995.

J. Ishida, H. K. Wang, K. F. Bastow, C. Q. Hu, and K. H. Lee, “Antitumor agents 201. Cytotoxicity of harmine and β-carboline analogs,” Bioorganic and Medicinal Chemistry Letters, vol. 9, no. 23, pp. 3319–3324, 1999.

M. Zhao, L. Bi, W. Wang et al., “Synthesis and cytotoxic activities of β-carboline amino acid ester conjugates,” Bioorganic and Medicinal Chemistry, vol. 14, no. 20, pp. 6998–7010, 2006.

R. Cao, W. Peng, H. Chen et al., “Synthesis and in vitro cytotoxic evaluation of 1,3-bisubstituted and 1,3,9-trisubstituted β-carboline derivatives,” The European Journal of Medicinal Chemistry, vol. 40, no. 3, pp. 249–257, 2005.

R. Cao, Q. Chen, X. Hou et al., “Synthesis, acute toxicities, and antitumor effects of novel 9-substituted β-carboline derivatives,” Bioorganic and Medicinal Chemistry, vol. 12, no. 17, pp. 4613–4623, 2004.

Q. Wu, R. Cao, M. Feng et al., “Synthesis and in vitro cytotoxic evaluation of novel 3,4,5-trimethoxyphenyl substituted β-carboline derivatives,” The European Journal of Medicinal Chemistry, vol. 44, no. 2, pp. 533–540, 2009.

H. Guan, H. Chen, W. Peng et al., “Design of β-carboline derivatives as DNA-targeting antitumor agents,” The European Journal of Medicinal Chemistry, vol. 41, no. 10, pp. 1167–1179, 2006.

R. Cao, W. Peng, H. Chen et al., “DNA binding properties of 9-substituted harmine derivatives,” Biochemical and Biophysical Research Communications, vol. 338, no. 3, pp. 1557–1563, 2005.

A. M. Deveau, M. A. Labroli, C. M. Dieckhaus, M. T. Barthen, K. S. Smith, and T. L. MacDonald, “The synthesis of amino-acid functionalized β-carbolines as topoisomerase II inhibitors,” Bioorganic and Medicinal Chemistry Letters, vol. 11, no. 10, pp. 1251–1255, 2001.

T. Herraiz and C. Chaparro, “Human monoamine oxidase enzyme inhibition by coffee and β-carbolines norharman and harman isolated from coffee,” Life Sciences, vol. 78, no. 8, pp. 795–802, 2006.

T. Herraiz and C. Chaparro, “Human monoamine oxidase is inhibited by tobacco smoke: β-carboline alkaloids act as potent and reversible inhibitors,” Biochemical and Biophysical Research Communications, vol. 326, no. 2, pp. 378–386, 2005.

K. Hoon, S. O. Sablin, and R. R. Ramsay, “Inhibition of monoamine oxidase A by β-carboline derivatives,” Archives of Biochemistry and Biophysics, vol. 337, no. 1, pp. 137–142, 1997.

J. Reniers, S. Robert, R. Frederick, B. Masereel, S. Vincent, and J. Wouters, “Synthesis and evaluation of β-carboline derivatives as potential monoamine oxidase inhibitors,” Bioorganic and Medicinal Chemistry, vol. 19, no. 1, pp. 134–144, 2011.

Y. Song, J. Wang, S. F. Teng et al., “β-carbolines as specific inhibitors of cyclin-dependent kinases,” Bioorganic and Medicinal Chemistry Letters, vol. 12, no. 7, pp. 1129–1132, 2002.

Y. Song, D. Kesuma, J. Wang et al., “Specific inhibition of cyclin-dependent kinases and cell proliferation by harmine,” Biochemical and Biophysical Research Communications, vol. 317, no. 1, pp. 128–132, 2004.

J. I. Trujillo, M. J. Meyers, D. R. Anderson et al., “Novel tetrahydro-β-carboline-1-carboxylic acids as inhibitors of mitogen activated protein kinase-activated protein kinase 2 (MK-2),” Bioorganic and Medicinal Chemistry Letters, vol. 17, no. 16, pp. 4657–4663, 2007.

Z. Chen, R. Cao, L. Yu et al., “Synthesis, cytotoxic activities and DNA binding properties of β-carboline derivatives,” The European Journal of Medicinal Chemistry, vol. 45, no. 11, pp. 4740–4745, 2010.

X. Han, J. Zhang, L. Guo et al., “A series of beta-carboline derivatives inhibit the kinase activity of PLKs,” PloS ONE, vol. 7, no. 10, pp. 1–11, 2012.

J. Zhang, Y. Li, L. Guo et al., “Dh166, a beta-carboline derivative, inhibits the kinase activity of PLK1,” Cancer Biology and Therapy, vol. 8, no. 24, pp. 2374–2383, 2009.

B. García-Álvarez, G. de Cárcer, S. Ibañez, E. Bragado-Nilsson, and G. Montoya, “Molecular and structural basis of polo-like kinase 1 substrate recognition: implications in centrosomal localization,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 9, pp. 3107–3112, 2007.

S. Cao, “QSAR, molecular docking studies of thiophene and imidazopyridine derivatives as polo-like kinase 1 inhibitors,” Journal of Molecular Structure, vol. 1020, no. 8, pp. 167–176, 2012.

R. K. Chourasiya, A. R. Rao, and R. K. Agrawal, “QSAR and docking studies of novel β-carboline derivatives as anticancer,” Medicinal Chemistry Research, vol. 22, no. 6, pp. 2991–3001, 2013.

G. Wu, D. H. Robertson III, C. L. Brooks, and M. Vieth, “Detailed analysis of grid-based molecular docking: a case study of CDOCKER-A CHARMm-based MD docking algorithm,” Journal of Computational Chemistry, vol. 24, no. 13, pp. 1549–1562, 2003.

M. Pastor, G. Cruciani, I. McLay, S. Pickett, and S. Clementi, “GRid-INdependent descriptors (GRIND): a novel class of alignment-independent three-dimensional molecular descriptors,” Journal of Medicinal Chemistry, vol. 43, no. 17, pp. 3233–3243, 2000.

F. Fontaine, M. Pastor, and F. Sanz, “Incorporating molecular shape into the alignment-free GRid-INdependent descriptors,” Journal of Medicinal Chemistry, vol. 47, no. 11, pp. 2805–2815, 2004.

Á. Durán, G. C. Martínez, and M. Pastor, “Development and validation of AMANDA, a new algorithm for selecting highly relevant regions in molecular interaction fields,” Journal of Chemical Information and Modeling, vol. 48, no. 9, pp. 1813–1823, 2008.

P. J. Goodford, “A computational procedure for determining energetically favorable binding sites on biologically important macromolecules,” Journal of Medicinal Chemistry, vol. 28, no. 7, pp. 849–857, 1985.

P. Cratteri, M. N. Romanelli, G. Cruciani, C. Bonaccini, and F. Melani, “GRIND-derived pharmacophore model for a series of α-tropanyl derivative ligands of the sigma-2 receptor,” Journal of Computer-Aided Molecular Design, vol. 18, no. 5, pp. 361–374, 2004.

R. Cao, X. Guan, B. Shi et al., “Design, synthesis and 3D-QSAR of β-carboline derivatives as potent antitumor agents,” The European Journal of Medicinal Chemistry, vol. 45, no. 6, pp. 2503–2515, 2010.

B. D. Hudson, R. M. Hyde, E. Rahr, J. Wood, and J. Osman, “Parameter based methods for compound selection from chemical databases,” Quantitative Structure-Activity Relationships, vol. 15, no. 4, pp. 285–289, 1996.

F. A. Momany and R. J. Rone, “Validation of the general purpose QUANTA 3.2/CHARMm force field,” Journal of Computational Chemistry, vol. 13, no. 7, pp. 888–900, 1992.

Discovery Studio, Accelrys Software, San Diego, Calif, USA, 2009.

L. Afzeliu, C. M. Masimirembwa, A. Karlén, T. B. Anderson, and I. Zamora, “Discriminant and quantitative PLS analysis of competitive CYP2C9 inhibitors versus non-inhibitors using alignment independent GRIND descriptors,” Journal of Computer-Aided Molecular Design, vol. 16, no. 7, pp. 443–458, 2002.