Received Date: August 30, 2017 Accepted Date: September 25, 2017 Published Date: Setember 27, 2017
Citation: Hammad S, Meghnem R, Meziane D, Bouaziz-Terrachet S (2017) In Silico Design of Protein Kinase 2 Inhibitors Using Molecular Docking and Pharmacophore Modelling. Int J Drug Dev & Res 9: 33-40
Flavonoids derivatives are a new class of natural substances that effectively inhibit CK2 protein. This protein represents an attractive pharmacological target for the treatment of cancer. With the aim to better understand the structural and chemical features responsible for the recognition mechanism of these compounds, and to explore their binding modes of interaction, a series of thirty-seven hydroxy flavones derivatives were docked into CK2 active site using FlexX software. The superposition of the docked conformation of the crystal ligand to that in the CK2 crystal structure are close since the RMSD between the two conformations was only 0.75 Å. The docking results show that important conserved residues, in particular Val116, Asp120, Arg47, Asn118, Ile66, Met163 and Val45, could be essential for the binding of the ligands to CK2 protein. The docking energies are in good agreement with the corresponding experimental activities PIC50 values of studied compounds since a significant correlation coefficient (R=0.866) was obtained. Then, the obtained docking conformations of some active compounds were used to generate pharmacophore models with the help of Ligand Scout. Therefore, the information provided in this study can be useful for designing new CK2 based on flavonoids compounds.
Casein kinase 2; CK2; Flavonoids; Breast cancer; Molecular docking; Pharmacophore modelling; Structure based drug design
Protein kinase II, also known as casein kinase II (CK2), is a constitutively active, ubiquitously expressed serine/threonine protein kinase. This protein exists as a heterotetrameric complex and is composed of an alpha, an alpha-prime, and two beta subunits. The alpha subunits contain the catalytic activity while the beta subunits undergo autophosphorylation. It has demonstrated that CK2 is involved in a wide variety of cellular processes such as gene expression, protein synthesis, proliferation, apoptosis and differentiation/ transformation [1-6]. Thus, perturbations in expression or activity of CK2 are associated with human disease [7-10]. Recently, it has demonstrated that this protein is involved in several human cancers, including cancer of the breast, lung, colon, and prostate, as well as hematologic malignancies [11,12]. In the light of these observations, CK2 represents an attractive pharmacological target for the treatment of cancer and the development of specific inhibitors of CK2. Therefore, the development and design of more promising CK2 inhibitors would be highly desirable.
A huge number of CK2 inhibitors are available, most of them are able to directly target the adenosine triphosphate (ATP) binding site, such as 5,6-Dichloro-1-β-D-ribofuranosyl-benzimidazole (DRB), emodin, tetrabromobenzotriazole (TBB), apigenin (AGI), and [5- oxo-5,6- dihydroindolo(1,2-a) quinazolin-7-yl] acetic acid (IQA) [13,14]. AGI is a flavonoid compound, belonging to a group of natural substances, which effectively inhibits CK2. This inhibitor, like most protein kinase inhibitors, is competitive with respect to the ATP. This indicates that AGI is directed towards the conserved ATP binding site . However, AGI is more potent, but as other general flavonoids, its specificity is quite broad [14,16]. Thus, there is a need to design new compounds which are able to inhibit CK2 with selectivity higher than those of AGI.
The aim of this investigation is to identify new CK2 inhibitors. In this context, we are interested in the ligand-protein binding, which is important not only as an essential molecular recognition process but also for the discovery of new CK2 inhibitors, between a series of flavonoids derivatives and CK2 protein. Molecular docking is an in-silico approach that is used to predict the best binding mode of compound within the active site of the specified enzyme. This technique has been applied previously in drug discovery [17-19] and many molecular modeling softwares are available among them the FlexX program .
In the present study, we performed a molecular docking study to search for best binding mode of a series of thirty-seven hydroxyflavones derivatives in the CK2 active site by using FlexX program. The obtained conformations of some active compounds are used to generate a pharmacophore model which can be applied to drug design and virtual screening [21,22].
Ligands preparation: The dataset used in this study, characterized by adequate biological and structural diversity, belongs to flavonoids family (Figure 1). Thirty-seven hydroxy flavones derivatives, which exhibited experimental activities ranged from 0.012 to 30 μM, were selected as ligands . These activities, are given in terms of IC50, the drug concentration that inhibits 50% of the enzymatic activity of CK2, which were converted to pIC50. The two-dimensional structure of these compounds (Table 1) were drawn with ISIS 2.4 package . The obtained geometries were energy optimized using the steepest descent approximation and were converged to a gradient of 0.02 using the UCSF Chimera software . Table 1 provides the negative logarithms of the observed activities, denoted PIC50.
Table 1: Structures and PIC50 of thirty seven flavonoids derivatives.
Receptor preparation: There are many CK2 complex crystal structures available in Protein Data Bank . Based on the resolution, the crystal structure of CK2 bound with AGI (Apigenin) (code 4DGM determined at 1.65 Å of resolution) was selected as a receptor . The protein was prepared using Receptor Preparation Wizard of LeadIT software from BioSolveIT, GmbH Germany (version 2.2.0)  by adding polar hydrogen.
Binding site determination: In LeadIT software, the selection of the binding site can be performed using two methods; the reference ligand and the selection of sphere to get the best score. In this work, the active site of the CK2 protein was determined based on the cocrystallized ligand (AGI) in the crystal structure. The active site radius was taken to be 6.5 Å from the center of co-crystallized ligand at the site of docking; so that the amino acid residues namely Val116, Glu114, Ile66, Met163, Val95, Val53, Val45 and Ile174 were included.
Docking protocol: The docking of flavonoids derivatives in the CK2 active site was performed using the FlexX docking algorithm  implemented in LeadIT program. The SDF file of all compounds was loaded in FlexX as a docking library. Then, the compounds were docked in the CK2 active site using default docking parameters and 50 highest scoring docked solutions of each compound were selected for further analysis. The interactions present in the selected dock poses were represented in 2D format using protein-ligand interaction map generating widget, PoseView  of LeadIT program. In order to validate obtained docking results, the correlation between the docking score values in Kcal/mol and the inhibitory activity of the selected compounds was computed.
Binding free energy calculations: FlexX is one of the best suited software for the analysis of macromolecule-ligand Interactions . The docking algorithm in FlexX uses flexible ligand and rigid protein . Ligand flexibility limited to torsion angles and ring conformations. This algorithm is based on the model of molecular interactions defined by Böhm [32,33]. This model can be divided into three areas: conformational flexibility, protein-ligand interactions, and the scoring function [20,34]. The ranking of the generated solutions is performed using a Böhm scoring function which estimates the free binding energy (ΔGbinding) of the protein-ligand complex [33-35]. ΔGbinding can be calculated using the equation 1 (Eq.1). ΔG0 corresponds to the loss of translational and rotational entropy of the ligand upon binding. ΔGpolar and ΔGapolar correspond to polar and a polar interaction, respectively. ΔGsolv and ΔGflexi are energies of solvation and flexibility of the ligand structure.
In the present study, the structure-based pharmacophore model was generated using the LigandScout software package  to identify the interaction of critical residues between the protein and ligand. LigandScout software uses an algorithm that studies and interprets ligand-receptor interactions and then automatically creates an advanced pharmacophore model. As input for structure-based pharmacophore generation, the obtained docked complexes of some active compounds with CK2 protein were chosen. Following this, a new structure based pharmacophore was generated by selecting the most relevant structural features of these compounds and removing the overlapped chemical features.
Comparison between the docked and X-ray crystal structures
The docked conformation of AGI was close to the crystal structure (Figure 2), since the RMSD between the two conformations was only 0.75 Å, which is quite satisfactory. This suggests that the FlexX docking method is reasonable and reliable. Therefore, the studied compounds can be docked correctly to the CK2 active site.
Docking interaction analysis
Molecular docking was carried out to interpret the best binding pose of the ligands and search the similar chemical entities that can bind strongly with the active site residues of CK2 protein. Based on the obtained docking results, ligands bind at the same ATP-binding site of CK2 protein as determined experimentally (Figure 3). Furthermore, the ligands poses are stabilized by a combination of hydrogen bonds and hydrophobic interactions with active site residues. Analysis of results revealed that water molecule Hoh581 and Val116 were the common residues in most interactions between carbonyl group (CO) of studied compounds and CK2 protein. Asp175, Lys68, Val45, Asn118, Arg47, Val116, Glu114 and Asp120 were the main amino acid residues involved in hydrogen bonding with the compounds which exhibited low activities (more potent compounds) toward CK2 (Figure 4). Nitro group (NO2) of compounds F1 and F4 was very active; it favors the hydrogen bond interactions. This group is pointed toward the positive charge region that is located near the amine function (NH3) of the Lys68 side chain. It interacts with protein CK2 forming hydrogen bonds with backbone nitrogen of Asp175 and amino group (NH3) of Lys68. This result is in good agreement with previous studies in which the authors demonstrate that the ligands that carry the negatively charged moiety adopt the conformation in which the negative group is pointed toward the positively charged region of the ATP binding site located in the buried area of the cavity (near the well conserved salt bridge Lys68- Glu81) . This region is responsible for polar interactions with the negatively charged moiety of the ligand, contributing to the different orientation of the ligand itself .
The hydroxyl oxygen at position 3’ (Figure 1) of compounds F2, F5, F13, F16 and F20 is hydrogen bonded to the backbone nitrogen of Asp120 and Arg47. The compounds F17, F16, F18 and F19 interact with protein forming hydrogen bond between hydroxyl group (-OH) and residue Asp120 (Figure 4). The backbone oxygen and backbone nitrogen of Arg47 showed strong hydrogen bonds with the methoxy group (-OMe) in position 2’ and hydroxyl oxygen in position 3’ of compounds F10, F17, F16, F18, F19, respectively. The residues Asn118, Gly46, Met163, Val45, Phe113, Val95, Val53, Ile174 and Ile66 made hydrophobic interactions with these compounds (Figure 4). The aromatic rings A of compounds F1, F4 and F16 are embedded into hydrophobic cavity consisting of Asn118, Gly46, Met163 and Val45. This indicates that all these compounds fit well in the binding pocket of CK2 protein which may explain their highest affinities toward CK2. It is worth noting that, for many of the studied compounds, hydrophobic interactions and hydrogen bonds are the most important driving force that stabilizes their binding to CK2 and, hence, are mainly responsible for the potency of the inhibitor . However, we judged the binding mode proposed by FlexX for these compounds as quite plausible.
The compounds which exhibited the PIC50 ranging from 0 to 1 (0<PIC50<1), interact with protein forming hydrogen bonds between their hydroxyl oxygen at 3’ position (2’ for compound F23) and residues Asp120 (Figure 5). The visualization of docked conformation of compound F22 shows that Asp120 had two hydrogen bonds, one with hydroxyl oxygen (–OH) in position 3’ and one other with methoxy group (–OCH3) at position 2’ (Figure 5). This allows multiple hydrophobic interactions to be made between the polycyclic core of the molecule and residues Asp120, Val53, Met163, Ile66, Val45, Ile174, Thr119 and Asn118 that normally form the environment of the adenosine part of ATP in the enzyme pocket. From Figure 3, we can see that the aromatic rings A of compounds F22 and F24 are surrounded by Ile66, Val45, Ile174, Val53 and Met163.
H-bonds formed by interaction of compounds with high activities (PIC50>1) and CK2 residues Val116 and Hoh581. The hydroxyl oxygen (-OH) at position R8 of compound F34 is involved in hydrogen bonding with backbone nitrogen of Val116 and backbone oxygen of Glu114. This group at position 3 of compounds F35 and F36 is hydrogen bonded with backbones oxygen and nitrogen of residue Val116 (same interactions are observed for compounds F31, F32 and F33), as shown in Figure 6. It is worth noting that, for some CK2 inhibitors which exhibited high activities, the preferred interaction is with the backbone carbonyls of Glu114 and Val116. Other previous docking studies of CK2 inhibitors also discovered according results . Moreover, Sarno et al.  clarified that backbone of residues Glu114 and Val116 was the key amino acid residue inside the binding site of CK2 responsible for the interaction with pentabromoindazole (K64). Residues Met163 is the common residue linked to the ring A of compounds F34, F35 and F36 via hydrophobic interactions.
Docking and experimental energies of selected compounds: The predicted docking energies score, obtained by FlexX program for studied compounds and the corresponding experimental PIC50 values are listed in Table 2. More negative values pointed out more favourable binding.
The PIC50 values of ligands can be grouped into three classes, active, inactive or moderately active, based on whether the PIC50 value is less than -1, more than 1 or ranging from -1 to 1, respectively. The performance of docking models is mostly measured in terms of their ability to discriminate between active and inactive compounds by calculating the docking energies. Thus, the actives and inactive compounds must be classified in active and inactive classes, respectively. Analysis of data presented in the Table 2 shows that compounds F1, F2, F4, F5 and F6 are well classified in the active class and compounds F2- F35 and F37 are well classified in the inactive class. Thus, we can say that the active molecules are 72% correctly classified, while the inactive molecules are 80% correctly classified. The correlation between the calculated energies score values of the studied compounds and the experimental PIC50 values is studied and plotted in Figure 7, which displays a good linear relationship. It is obvious from this figure that all the docking energies are in good agreement with the corresponding experimental activities since a significant correlation coefficient R=0.866 (R2=0.750) was obtained.
Table 2: Docking energies score (kcal mol-1) and PIC50 values of the studied compounds.
Structure based pharmacophore modelling
Obtained docking results were imported in LigandScout. Molecule and side chains of each obtained complex were energy minimized using MMFF94 to refine the FlexX output. Then, important molecular interactions to binding site residues of CK2 protein such as hydrogen bonds, hydrophobic interactions, and so on were analyzed in order to identify important features for ligand binding. The obtained pharmacophore models show interactions between important amino acid residues in the binding site and ligand and then elucidate the favorable binding position of the docked compounds. Pharmacophore features in the generated models consist of hydrogen bond donor HBD (green arrow), hydrogen bond acceptor HBA (red arrow) and hydrophobic interaction Hy (yellow sphere); as shown in Figure 8 and listed in Table 3.
|mol||R2’||R3’||R4’||R6||R7||R8||Ring A||Ring C|
Table 3: Pharmacophoric features of each well classified compound in the active class.
The final pharmacophore features at each position in the common structure of the five well classified compounds, in the active class (compounds F1, F2, F4, F5 and F6) (Figure 8), were analyzed and clustered according to their interaction pattern with the CK2 protein. The most relevant structural features of these compounds were selected and included in the final pharmacophore model.
The analysis of the results presented in Table 3 shows that HBA and HBD in 3’ position is present in almost all compounds (four out of five). Three compounds out of five make an HBA between carbonyl group (CO) and CK2 active site residues. Two hydrophobic interactions between two rings A and C and CK2 protein should be taken into consideration. Hense, the final pharmacophore model consists of two HBA at positions 3’ and 4, one HBD at position 3’ and two hydrophobic interactions centers, ring A and ring C (Figure 9).
Figure 9: Pharmacophore hypothesis for flavonoids as CK2 inhibitors. (a) Pharmacophore model with five features: two hydrogen bond acceptors (HBA), one hydrogen bond donors (HBD) and two hydrophobic group (H); (b) predicted pharmacophore features distance range. All features in this pharmacophore model are color coded with green: HBA; red: HBD; blue spheres: Hy. Distance between certain features is represented in Å.
Thirty-seven flavonoids derivatives were docked on CK2 active site using the FlexX software. The docked conformation of the apigenin, extracted from the crystallographic complex 4DGM, was close to the crystal structure. Visualization of docking results shows that flavonoids derivatives interact with CK2 protein predominantly through hydrogen bonding. The majority of studied compounds were observed to form hydrogen bonds between either the carbonyl group (CO) of each compound and Val116 of the CK2 and with water molecule Hoh581 or between the hydroxyl oxygen and Asp120 and Arg47. The residues Asn118, Asp120, Ile66, Met163 and Val45 are linked to the ligands via hydrophobic interactions with the phenyl rings A and C. In order to add robustness to the employed docking model, the correlation between the docking energies of studied compounds and experimental PIC50 values was studied and plotted. The results show that the docked energies are in good agreement with the corresponding experimental activities since a significant correlation coefficient (R=0.866) was obtained.
Obtained docking conformations of the well classified active compounds were used to generate pharmacophore models with the help of LigandScout. the final pharmacophore model consists of two HBA at positions 3’ and 4, one HBD at position 3’ and two hydrophobic interaction at ring A and ring C. Therefore, the information provided in this study can be useful for designing new CK2 based on flavonoids compounds.
The authors would like to thank Inte:Ligand Software-Entwicklungs und Consulting GmbH and BioSolveIT team for providing an academic free license for LigandScout 4.1.7 and LeadIT 2.2.0 programs.