Next Article in Journal
The Role of Microglia in Diabetic Retinopathy: Inflammation, Microvasculature Defects and Neurodegeneration
Next Article in Special Issue
Computational Characterization of Small Molecules Binding to the Human XPF Active Site and Virtual Screening to Identify Potential New DNA Repair Inhibitors Targeting the ERCC1-XPF Endonuclease
Previous Article in Journal
Candidate Genes and Molecular Markers Correlated to Physiological Traits for Heat Tolerance in Fine Fescue Cultivars
Previous Article in Special Issue
Structure and Molecular Dynamics Simulations of Protein Tyrosine Phosphatase Non-Receptor 12 Provide Insights into the Catalytic Mechanism of the Enzyme
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 1
10.3390/ijms19010111
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies

by
Yue Zhou
1,
Na Zhang
2,*,
Xiaoqian Qi
2,
Shan Tang
2,
Guohui Sun
2,
Lijiao Zhao
2,
Rugang Zhong
2 and
Yongzhen Peng
1,*
1
National Engineering Laboratory for Advanced Municipal Wastewater Treatment and Reuse Technology, Engineering Research Center of Beijing, Beijing University of Technology, Beijing 100124, China
2
Beijing Key Laboratory of Environmental & Viral Oncology, College of Life Science and Bioengineering, Beijing University of Technology, Beijing 100124, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(1), 111; https://doi.org/10.3390/ijms19010111
Submission received: 3 November 2017 / Revised: 20 December 2017 / Accepted: 29 December 2017 / Published: 1 January 2018
(This article belongs to the Special Issue Drug-Protein Interactions and Mechanisms of Action by Structural Modifications Driving Protein Function)
Browse Figures
Versions Notes

Abstract

:
Protein kinase is a novel therapeutic target for human diseases. The off-target and side effects of ATP-competitive inhibitors preclude them from the clinically relevant drugs. The compounds targeting the druggable allosteric sites outside the highly conversed ATP binding pocket have been identified as promising alternatives to overcome current barriers of ATP-competitive inhibitors. By simultaneously interacting with the αD region (new allosteric site) and sub-ATP binding pocket, the attractive compound CAM4066 was named as allosteric inhibitor of CK2α. It has been demonstrated that the rigid linker and non-ionizable substituted fragment resulted in significant decreased inhibitory activities of compounds. The molecular dynamics simulations and energy analysis revealed that the appropriate coupling between the linker and pharmacophore fragments were essential for binding of CAM4066 with CK2α. The lower flexible linker of compound 21 lost the capability of coupling fragments A and B to αD region and positive area, respectively, whereas the methyl benzoate of fragment B induced the re-orientated Pre-CAM4066 with the inappropriate polar interactions. Most importantly, the match between the optimized linker and pharmacophore fragments is the challenging work of fragment-linking based drug design. These results provide rational clues to further structural modification and development of highly potent allosteric inhibitors of CK2.

Graphical Abstract

1. Introduction

Protein kinases play predominant regulatory roles in nearly every aspect of cellular function, and have been considered as one of the most attractive drug targets [1,2]. Some FDA-approved anti-cancer drugs of protein kinase, such as Gefitinib, Palbociclib, and Tofacitinib, are ATP-competitive inhibitors. However, most of these inhibitors targeting highly conserved ATP-binding pocket of the protein kinase family were impeded from the clinical-used drugs due to the risk of the off-target and side effects [3,4,5]. It is urgent to discover allosteric inhibitors targeting novel druggable sites outside the catalytic box. Given the over-expression in a range of cancer cell lines, protein kinase CK2 has been regarded as a representative of cancer therapeutic targets [6,7,8]. To date, many ATP-competitive inhibitors with different scaffolds have been developed, such as benzimidazole derivatives, anthraquinone derivatives, tricyclic quinolone derivatives and natural products [9,10,11]. Despite their anti-proliferation effects on malignant tumors, these inhibitors are still facing a stumbling block of lower specificity and diversity deficiencies to be rejected as clinical drugs [12,13,14]. Therefore, it is a challenging task to develop inhibitors with higher affinity and selectivity by targeting the allosteric sites, which may avoid the drawbacks of most conventional ATP-competitive kinase inhibitors. For instance, DRB and W16 have been identified as allosteric inhibitors of CK2. Both are able to disrupt the tetrameric assembly at the CK2α/β interface [15,16]. Recently, Jiang et al. identified another four potential allosteric pockets of CK2α and predicated allosteric pathways between allosteric sites and the active site through bioinformatics methods, such as site 2 (Leu249, Trp281 and tye307), site 3 (Tye125, Met225 and Leu128) and site 4 (His148, Tyr211 and Thr314) as well as site 5 (Trp33, Lys75, Ile78 and Pro109) [17]. These potent sites inspire researchers to design and synthesize the novel allosteric inhibitors of CK2α.
Among these allosteric sites, site 3 is located in the αD region next to the ATP binding pocket. The αD region of CK2α presents the unique features in contrast to other protein kinases, including much more flexibility, partially open conformation different from the closed status of other kinases and the non-conserved residues throughout the whole CMGC family, which includes cyclin-dependent kinases (CDKs), mitogen-activated protein kinases (MAPKs), glycogen synthase kinases (GSKs), and CDK-like kinases (CLKs), and thus appears to be an attractive site for the novel selective CK2 allosteric inhibitors design [17,18,19]. Recently, Brear and Fusco et al. discovered a new potent inhibitor CAM4066 that targeted a previously unseen pocket “αD pocket” and ATP site, simultaneously using the fragments linking strategy [20]. Notably, the comparison of the structure and inhibitory activity of compounds 21 (IC50 = n/a) and Pre-CAM4066 (IC50 = n/a) versus compound CAM4066 (IC50 = 370 nM) indicated that the inappropriate linker and fragments resulted in a significant reduction of inhibitory activities. Therefore, there is an urgent need to systematically investigate the structural basis for the decreased inhibitory activities of CAM4066 derivatives.
In this study, comparative molecular dynamics (MD) analyses were conducted to examine the detailed binding modes of the three compounds (CAM4066, Pre-CAM4066, and compound 21) to CK2α. In addition, molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) binding free-energy calculations were performed to elucidate the instability of the compound 21 and Pre-CAM4066 systems. Our findings may provide valuable information for further structural modification and development of highly potent allosteric inhibitors of CK2.

2. Results

2.1. Molecular Dynamics Simulation Studies

2.1.1. Overall Features of Dynamic Behaviors

The dynamic stability and behavior of three systems were explored by 50 ns conventional MD simulations in explicit water. Figure 1A shows the time-dependent root mean square deviation (RMSD) profile of each trajectory from the minimized structure. The relatively stable RMSD values during the last 20 ns indicated the approaching acceptable equilibrium state. As depicted in the lower panel, the fluctuations of RMSD data of CK2α in all systems followed the similar trend throughout the simulations. The average RMSD values of compound 21 (2.5 Å) and pre-CAM4066 (3.0 Å) systems were higher than those of CK2α–CAM4066 complex (2.0 Å), which suggested that the CK2α structures in the two systems underwent conformational changes to a certain extent. In contrast, some slightly different behaviors could be observed in the upper panel. The RMSD values of CAM4066 remained around 1.0 Å during the MD simulation, while the corresponding data of compound 21 and Pre-CAM4066 gradually reduced from 2.0 Å to 1.5 Å and increased from 1.0 Å to 1.5 Å, respectively. It was speculated that rigid linker of compound 21 and non-ionizable fragment B of Pre-CAM4066 resulted in the inappropriate interactions between ligands and CK2, which induced that the two compounds may deviate away from its original position.
To identify the detailed flexible regions in the protein structure during MD simulation, the thermodynamic stability of the three systems characterized by B-factor are shown in Figure 1B. Judging from the high B-factor values in the region of the N-terminal, G-loop, β4/β5-loop and C-terminal for three systems, these residues appeared to exhibit large fluctuations, as observed in other CK2α-ligand complexes, which have been demonstrated in previous studies [19,21,22,23]. It was speculated that the higher B-factors of αD/hinge region of compound 21 system resulted from the varied binding of fragment A. For Pre-CAM4066 system, as a consequence of the decreased interactions between methyl benzoate substituent and Lys68, the C-loop next to the positive area was more flexible in contrast to other systems.

2.1.2. Comparative Analysis of Different Binding Modes of Three Compounds

As indicated in the crystallized conformation, with the flexible N-(3-oxobutyl)butyramide scaffold as a bridge, CAM4066 simultaneously extended to the positive area of ATP binding pocket and the allosteric site αD region by forming hydrophobic and polar interactions. MD simulation provided the detailed dynamical behaviors of the CK2α-CAM4066 complex (Figure 2). Fragment A was entrapped into the cavity of αD region, where the biphenyl group formed stable hydrophobic packing with side chains of residues Phe121, Leu124, Tyr 125, Leu128, Ile133, Tyr136, Met137, Ile140, Met221 and Met225. The alternative H-bonds between Pro159 (66.59% occurrence) and Val162 (63.87% occurrence) and the NH of fragment A were considered to be an anchor for the binding of fragment A into αD pocket. The carboxylate substituent of fragment B tended to create electrostatic interactions with residue Lys68 of positive area. In contrast to the crystallized CK2α-CAM4066 complex, besides the subtle positional deviation of the linker, the new H-bond (62.54% occurrence) between the NH group of fragment B and the carbonyl oxygen atom of Leu45 was identified. It is speculated that the stable polar interactions between fragment B and Lys68 induced the reorientation of O3 and O4 atoms of the linker, which resulted in the pairing H-bonds between N2 and O3 of linker and Asn118 replaced by the H-bond between O4 atom of linker and NH2 of Asn118 (32.73% occurrence).
Given the limited shortcomings of X-ray crystallization method, the static conformation of compound 21 could not unveil the underlying interaction mechanism with CK2α. Owing to the development of molecular modeling methods, MD analysis can be used to interpret and complement the experimental phenomena. As shown in Figure 3A,B, although compound 21 occupied the same binding pocket as CAM4066 did, the twisted linker and re-orientated fragments A and B can be found easily. In the complex of CK2α-CAM4066, the flexible linker had the capability to couple fragment A and B binding into the αD region and positive area, respectively. With the replacement of ethyl with the methyl group, the rigid linker of compound 21 was incapable of capturing enough conformation to guarantee two fragments targeting two regions simultaneously as the flexible linker of CAM4066 could. It was also supported by the dihedrals angles analysis shown in the Figure 4. During the last 20 ns, the dihedral angles of N1-C8-C9-C10 and C8-C9-C10-N2 of compound 21 were kept fixed near −75° and 150°, respectively, whereas the corresponding two dihedral angles of CAM4066 presented cooperative rotation behaviors. It could be concluded that the flexible linker of CAM4066 covered more chemical spaces than the compound 21 did. The N2 atom of the twisted linker established a stable H-bond with the carbonyl oxygen atom of His160 (85.84% occurrence). Meanwhile, fragment B was also deprived of the ability to make polar interactions with positive area. Fragment A was still entrapped into the hydrophobic pocket of αD region by the alternative H-bond among N3 atom, Val162 (78.70% occurrence) and Pro159 (41.08% occurrence).
As only the methyl benzoate of fragment B is the subtle difference between Pre-CAM4066 and CA4066, the position and orientation of Pre-CAM4066 was similar to those of CAM4066 except the deviated fragment B. As shown in Figure 3C,D, the linkers of two compounds were nearly overlapped. Fragment A of Pre-CAM4066 still located in αD pocket by forming the stable H-bonds with Pro159 and Val162 (71.98% and 60.35% occurrence, respectively), whereas the methyl benzoate of fragment B could not make any interactions with Lys68 and rotated out of the positive area. As triggered by the movement of fragment B, a stable H-bond was established between O3 atom of the deviated linker and NH2 group of Asn118 (85.70% occurrence). To sum up, both the flexible linker and ionizable substituent of fragment B are two determinants for binding of CAM4066. As indicated from the stable interactions between fragment A and αD region of three systems, a fact worth highlighting is that fragment A should be identified as a new pharmacophore group of αD region.

2.2. Energy Analysis

To elucidate the quantitative effect of flexible linker and ionizable substituted fragment B on the compounds binding affinity, molecular mechanics/Poisson−Boltzmann surface area (MM/PBSA) energy terms for each system were calculated, as listed in Table 1. The total ∆Gbinding of the CAM4066 (IC50 = 0.370 μM, Kd = 0.320) was −54.68 kcal/mol, which was 3.25 and 6.51 kcal/mol lower than those of compound 21 (IC50 = n/a, Kd = 1.64) and Pre-CAM4066 (IC50 = n/a, Kd = n/a), respectively. These data indicated that CAM4066 exhibited the highest affinity binding to CK2. Further analysis of the energy components responsible for the binding free energies showed that ∆Eele, ∆Evdw, and ∆Enon-polar was the main driving forces for the three systems. The ∆Eele values of compound CAM4066 (−110.41 kcal/mol) exhibited more favorable contributions than those of compound 21 (−53.42 kcal/mol) and Pre-CAM4066 (−11.17 kcal/mol), which was in accordance with the loss of electrostatic interactions between fragment B and the positive area of two systems. Taking the polar contribution of ∆Gpolar into consideration, the total electrostatic interaction energy (∆Gele) was positive for three compounds. This finding may be interpreted by the fact that the favorable contribution of the ∆Eele was more than compensated by the ∆Gpolar upon binding. Thus, ∆Gpolar was unfavorable to this class of complexes.

3. Discussion

Fragment-based drug discovery (FBDD) is a robust tool for identify new anti-cancer leads and clinical-used drugs, such as Vemurafenib and Venetoclax [24]. Rather than screening millions of compounds as conducted by high-throughput screening, FBDD begins with the collection of smaller compounds or fragments with high affinity to favorable spots [25]. The small size of the fragment library is much more efficient in covering the chemical spaces than the amount of compound library does [26]. Hence the fragment collection is the most fundamental step for FBDD.
Fusco et al. [20] obtained two effective fragment libraries for the unreported αD site and positive area of ATP sub-pocket with the experimental screening method. As the αD site cavity consisted of hydrophobic residues Phe121, Leu124, Tyr125, Leu128, Ile133, Tyr136, Met137, Ile140, Met221 and Met225, and was also highly solvent accessible, hydrophobic fragments were entrapped into the cavity of αD region and formed two pair H-bonds with carbonyl oxygen atom of Pro159 and Val162, respectively. Meanwhile, these fragments also make the contribution to the selectivity of compounds by specific interacting with αD site rather than the hinge region of ATP binding pocket. Recently, there is increasing interests in the identification of pharmacophore fragments using efficient and economical virtual splitting and screening technologies. In our study, the reported co-crystallized ATP-competitive inhibitors were decomposed into three types of fragments targeting CK2α ATP binding sub-pocket (hydrophobic, hinge and positive area) by a pharmacophore oriented fragmentation algorithm (unpublished). The pharmacophore fragment library targeting sub-binding pocket would enrich the ATP site fragment database and provide structural elements for the rational design of CK2 inhibitor. Similarly, Zhang and Zhao discovered of a novel series of CK2 inhibitors by using a library of virtual fragments with key functionalities via fragmentation of bioactive molecules [27].
Fragment linking is the useful strategy of FBDD, in which two fragments that bind to a distinct site are joined together by the suitable linker. Thus, the length and conformational flexibility of the linker are two decisive factors that influence the binding of fragments. A serial of linkers were obtained by optimization of the alkyl chain and polar groups, and the (N-(3-oxobutyl)butyramide) was chosen as the best candidate. Besides the H-bonds between the best linker and residue Asn118, this linker also had no negative effects on the conformation and binding modes of fragment A and B. This suggested that the linker of CAM4066 met the requirement proposed by Kim [25]. Any variation exiting in the flexible linker may exert the detrimental effect on the binding modes of compounds. For instance, the rigid linker with methyl group of compound 21 destroyed the original proper interactions the αD region and positive area of CK2α t as CAM4066 did. Furthermore, chemical environment and structural topology of two sub-pockets were two other key factors for linker optimization.

4. Materials and Methods

4.1. System Setup and Parameters Preparation

The atomic co-ordinates of CK2α-CAM4066 and CK2α-compound 21 were retrieved directly from the Protein Data Bank [PDB ID: 5CU4 and 5MO8] and hetero atoms, expected for water molecules within 6.5 Å of the ligand, were removed [18,20]. The compound Pre-CAM4066 was constructed based on the compound CAM4066 by using the SYBYL 8.1 program and the Tripos force field was used to energetically minimize the compound [28,29]. The partial atomic charges of three compounds were obtained via quantum electronic structure calculations including an optimization procedure using the Gaussian 03 program at the HF/6-31G* level, electrostatic potential (ESP) generation using the MerzSingh–Kollman van der Waals parameters, the atom-centered charge fitting through the RESP program implemented in the AMBER 10 package [30,31,32]. Subsequently, each system was neutralized by adding suitable counterions and then solvated in a truncated octahedral box of TIP3P water molecules [33,34]. The chemical structures and IC50 values of three compounds are listed in the Table 2.

4.2. Molecular Dynamics Simulations

Molecular dynamics simulations were initiated on the CAM4066, Pre-CAM4066 and compound 21 and each simulation was performed for 50 ns using the Amber 10 package [35]. The force field parameters for protein and ligands were calculated by the AMBER FF03 force field and the general AMBER force field (GAFF), respectively [36,37]. First, the geometric strain and close intermolecular contacts were relieved in the energy minimizations using the steepest descent and conjugate gradient methods. Second, each energy-minimized structure was gradually warmed from 0 to 300 K with weak constraint to the complex (5.0 kcal/mol) over 15 ps, followed by constant temperature equilibration at 300 K for 35 ps with constant volume dynamics. Third, MD simulations were carried out with the periodic boundary condition in the NPT ensemble, using a non-bonded cutoff of 10 Å to truncate the VDW non-bonded interactions [38]. Temperature (300 K) and constant pressure (1 atm) were maintained by Langevin dynamics temperature coupling with a time constant of 1.0 ps and isotropic position scaling with a relaxation time of 2.0 ps, respectively. The long-range electrostatic interactions were calculated based on the particle-mesh Ewald (PME) algorithm, and the SHAKE algorithm was applied to constrain all bonds involving hydrogen atom [39,40].

4.3. MM/PBSA Calculations

The MM-PBSA methods was employed to evaluate the three compounds binding energies and the effects of flexibility of linker and ionizable substituted fragment on the compounds binding from an energetic view [41,42]. For each system, the binding energy (ΔGbinding) was calculated for the configurations taken from a single trajectory based on the following equation:
ΔGbinding = Gcomplex − (Gprotein + Gligand) = ΔEgas + ΔGsol − TΔS
where the gas molecular mechanical energy (ΔEgas) is calculated as a sum of internal energies (i.e., bond, angle, and dihedral), van der Waals (Evdw) and electrostatic energies (Eele) using the SANDER module without applying a cutoff for non-bonded interactions. The solvation free energy (ΔGsol) is composed of electrostatic (ΔGpolar) and non-polar (ΔGnon-polar) contributions. The electrostatic contribution to the solvation free energy (ΔGpolar) is determined by PB model as implemented in SANDER, applying dielectric constants of 1 and 80 to represent the solute and the exterior medium phases, respectively. The non-polar component (ΔGnon-polar) is calculated using a linear function of solvent-accessible surface area (SASA) as follows: ΔGnon-polar = λ·SASA + b, where the corresponding parameters λ and b are set to 0.00542 kcal/(mol Å2)and 0.92 kcal/mol, respectively [43]. Given the large computational overhead and low prediction accuracy, the time consuming conformational entropy change (−TΔS) was not considered [44,45]. The entropy term has been neglected, assuming that it will be very similar for all of the systems.

5. Conclusions

MD simulations and energy calculations were performed to elucidate the structural mechanisms through which the rigid linker and non-ionizable substituted fragment influence binding affinity. It seemed that the optimized linker was not only the bridge of the two pharmacophore fragments, but also the adjustor for the binding of fragments into sub-pockets. Both the linker of compound 21 and fragment B of Pre-CAM4066 could not form the proper interactions with CK2α as those of CAM4066, whereas fragment A of three systems maintained stable interactions with αD region of CK2α. In addition, the energy analysis enabled the qualitative investigation of the effect of flexible linker and ionizable substituted fragment B on the three complexes. This will provide the theoretical basis and experiment guidance for the development of potent allosteric inhibitors of CK2.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (Grant No. 81402852), the Beijing Natural Science Foundation (Grant No. 7152013), the Beijing Municipal Commission of Education Research Projects (Grant No. KM201610005031), the China Scholarship Council awardee (Grant No. 201706545010), Beijing Postdoctoral Science Foundation (Grant No. 2017-ZZ-015) and Chaoyang District Postdoctoral Science Foundation (Grant No. 2017ZZ-01-26). We acknowledge Professor David Case for the kind gift of AMBER 10 software.

Author Contributions

Yue Zhou and Na Zhang conceived and designed the experiments; Yue Zhou and Na Zhang performed the experiments and analyzed the data; Xiaoqian Qi, Shan Tang, Guohui Sun, Lijiao Zhao, Rugang Zhong and Yongzhen Peng contributed analysis tools and helped in the Results and Discussion Section; and Yue Zhou and Na Zhang wrote the paper. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roskoski, R. Classification of small molecule protein kinase inhibitors based upon the structures of their drug-enzyme complexes. Pharmacol. Res. 2016, 103, 26–48. [Google Scholar] [CrossRef] [PubMed]
  2. Roskoski, R. A historical overview of protein kinases and their targeted small molecule inhibitors. Pharmacol. Res. 2015, 100, 1–23. [Google Scholar] [CrossRef] [PubMed]
  3. Klaeger, S.; Gohlke, B.; Perrin, J.; Gupta, V.; Heinzlmeir, S.; Helm, D.; Qiao, H.; Bergamini, G.; Handa, H.; Savitski, M.M.; et al. Chemical Proteomics Reveals Ferrochelatase as a Common Off-target of Kinase Inhibitors. ACS Chem. Biol. 2016, 11, 1245–1254. [Google Scholar] [CrossRef] [PubMed]
  4. Schnieders, M.J.; Kaoud, T.S.; Yan, C.; Dalby, K.N.; Ren, P. Computational insights for the discovery of non-ATP competitive inhibitors of MAP kinases. Curr. Pharm. Des. 2012, 18, 1173–1185. [Google Scholar] [CrossRef] [PubMed]
  5. Wu, P.; Nielsen, T.E.; Clausen, M.H. FDA-approved small-molecule kinase inhibitors. Trends Pharmacol. Sci. 2015, 36, 422–439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ortega, C.E.; Seidner, Y.; Dominguez, I. Mining CK2 in cancer. PLoS ONE 2014, 9, e115609. [Google Scholar] [CrossRef] [PubMed]
  7. Filhol, O.; Giacosa, S.; Wallez, Y.; Cochet, C. Protein kinase CK2 in breast cancer: The CK2beta regulatory subunit takes center stage in epithelial plasticity. Cell. Mol. Life Sci. 2015, 72, 3305–3322. [Google Scholar] [CrossRef] [PubMed]
  8. Guerra, B.; Issinger, O.G. Protein kinase CK2 in human diseases. Curr. Med. Chem. 2008, 15, 1870–1886. [Google Scholar] [CrossRef] [PubMed]
  9. Cozza, G.; Bortolato, A.; Moro, S. How druggable is protein kinase CK2? Med. Res. Rev. 2010, 30, 419–462. [Google Scholar] [CrossRef] [PubMed]
  10. Cozza, G.; Pinna, L.A.; Moro, S. Kinase CK2 Inhibition: An Update. Curr. Med. Chem. 2013, 20, 671–693. [Google Scholar] [CrossRef] [PubMed]
  11. Battistutta, R. Protein kinase CK2 in health and disease: Structural bases of protein kinase CK2 inhibition. Cell. Mol. Life Sci. 2009, 66, 1868–1889. [Google Scholar] [CrossRef] [PubMed]
  12. Sarno, S.; Papinutto, E.; Franchin, C.; Bain, J.; Elliott, M.; Meggio, F.; Kazimierczuk, Z.; Orzeszko, A.; Zanotti, G.; Battistutta, R.; et al. ATP Site-Directed Inhibitors of Protein Kinase CK2: An Update. Curr. Top. Med. Chem. 2011, 11, 1340–1351. [Google Scholar] [CrossRef] [PubMed]
  13. Li, Y.; Luan, Y.; Qi, X.M.; Li, M.; Gong, L.K.; Xue, X.A.; Wu, X.F.; Wu, Y.F.; Chen, M.; Xing, G.Z.; et al. Emodin Triggers DNA Double-Strand Breaks by Stabilizing Topoisomerase II-DNA Cleavage Complexes and by Inhibiting ATP Hydrolysis of Topoisomerase II. Toxicol. Sci. 2010, 118, 435–443. [Google Scholar] [CrossRef] [PubMed]
  14. Pagano, M.A.; Bain, J.; Kazimierczuk, Z.; Sarno, S.; Ruzzene, M.; Di Maira, G.; Elliott, M.; Orzeszko, A.; Cozza, G.; Meggio, F.; et al. The selectivity of inhibitors of protein kinase CK2: An update. Biochem. J. 2008, 415, 353–365. [Google Scholar] [CrossRef] [PubMed]
  15. Raaf, J.; Brunstein, E.; Issinger, O.G.; Niefind, K. The CK2 alpha/CK2 beta interface of human protein kinase CK2 harbors a binding pocket for small molecules. Chem. Biol. 2008, 15, 111–117. [Google Scholar] [CrossRef] [PubMed]
  16. Laudet, B.; Moucadel, V.; Prudent, R.; Filhol, O.; Wong, Y.S.; Royer, D.; Cochet, C. Identification of chemical inhibitors of protein-kinase CK2 subunit interaction. Mol. Cell. Biochem. 2008, 316, 63–69. [Google Scholar] [CrossRef] [PubMed]
  17. Jiang, H.M.; Dong, J.K.; Song, K.; Wang, T.D.; Huang, W.K.; Zhang, J.M.; Yang, X.Y.; Shen, Y.; Zhang, J. A novel allosteric site in casein kinase 2alpha discovered using combining bioinformatics and biochemistry methods. Acta Pharmacol. Sin. 2017, 38, 1691–1698. [Google Scholar] [CrossRef] [PubMed]
  18. Brear, P.; De Fusco, C.; Hadje Georgiou, K.; Francis-Newton, N.J.; Stubbs, C.J.; Sore, H.F.; Venkitaraman, A.R.; Abell, C.; Spring, D.R.; Hyvonen, M. Specific inhibition of CK2alpha from an anchor outside the active site. Chem. Sci. 2016, 7, 6839–6845. [Google Scholar] [CrossRef] [PubMed]
  19. Papinutto, E.; Ranchio, A.; Lolli, G.; Pinna, L.A.; Battistutta, R. Structural and functional analysis of the flexible regions of the catalytic alpha-subunit of protein kinase CK2. J. Struct. Biol. 2012, 177, 382–391. [Google Scholar] [CrossRef] [PubMed]
  20. De Fusco, C.; Brear, P.; Iegre, J.; Georgiou, K.H.; Sore, H.F.; Hyvonen, M.; Spring, D.R. A fragment-based approach leading to the discovery of a novel binding site and the selective CK2 inhibitor CAM4066. Bioorg. Med. Chem. 2017, 25, 3471–3482. [Google Scholar] [CrossRef] [PubMed]
  21. Zhou, Y.; Li, X.; Zhang, N.; Zhong, R. Structural basis for low-affinity binding of non-R2 carboxylate-substituted tricyclic quinoline analogs to CK2alpha: Comparative molecular dynamics simulation studies. Chem. Biol. Drug Des. 2015, 85, 189–200. [Google Scholar] [CrossRef] [PubMed]
  22. Zhou, Y.; Zhang, N.; Zhong, R.G. Exploring the crucial structural elements required for tricyclic quinoline analogs as protein kinase CK2 inhibitors by a combined computational analysis. Med. Chem. Res. 2013, 22, 4410–4422. [Google Scholar] [CrossRef]
  23. Zhang, N.; Zhong, R.G. Structural basis for decreased affinity of Emodin binding to Val66-mutated human CK2 alpha as determined by molecular dynamics. J. Mol. Model. 2010, 16, 771–780. [Google Scholar] [CrossRef] [PubMed]
  24. Keseru, G.M.; Hann, M.M. What is the future for fragment-based drug discovery? Future Med. Chem. 2017, 9, 1457–1460. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, J.; Moon, Y.; Hong, S. Identification of lead small molecule inhibitors of glycogen synthase kinase-3 beta using a fragment-linking strategy. Bioorg. Med. Chem. Lett. 2016, 26, 5669–5673. [Google Scholar] [CrossRef] [PubMed]
  26. Mondal, M.; Radeva, N.; Fanlo-Virgos, H.; Otto, S.; Klebe, G.; Hirsch, A.K. Fragment Linking and Optimization of Inhibitors of the Aspartic Protease Endothiapepsin: Fragment-Based Drug Design Facilitated by Dynamic Combinatorial Chemistry. Angew. Chem. 2016, 55, 9422–9426. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, N.; Zhao, H. Enriching screening libraries with bioactive fragment space. Bioorg. Med. Chem. Lett. 2016, 26, 3594–3597. [Google Scholar] [CrossRef] [PubMed]
  28. CERTARA. Available online: https://www.certara.com/ (accessed on 19 May 2017).
  29. Clark, M.; Cramer, R.D.; Vanopdenbosch, N. Validation of the General-Purpose Tripos 5.2 Force-Field. J. Comput. Chem. 1989, 10, 982–1012. [Google Scholar] [CrossRef]
  30. Fox, T.; Kollman, P.A. Application of the RESP methodology in the parametrization of organic solvents. J. Phys. Chem. B 1998, 102, 8070–8079. [Google Scholar] [CrossRef]
  31. Besler, B.H.; Merz, K.M.; Kollman, P.A. Atomic Charges Derived from Semiempirical Methods. J. Comput. Chem. 1990, 11, 431–439. [Google Scholar] [CrossRef]
  32. Gaussian.com. Available online: http://gaussian.com/ (accessed on 19 May 2017).
  33. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  34. Zhuang, S.; Wang, H.; Ding, K.; Wang, J.; Pan, L.; Lu, Y.; Liu, Q.; Zhang, C. Interactions of benzotriazole UV stabilizers with human serum albumin: Atomic insights revealed by biosensors, spectroscopies and molecular dynamics simulations. Chemosphere 2016, 144, 1050–1059. [Google Scholar] [CrossRef] [PubMed]
  35. Case, D.A.; Cheatham, T.E., III; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  37. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef] [PubMed]
  38. Berendsen, H.J.; Postma, J.V.; van Gunsteren, W.F.; DiNola, A.R.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  39. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H.J.C. Numerical integration of the cartesian equations of motion of a system with constraints:molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 321–341. [Google Scholar] [CrossRef]
  40. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef]
  41. Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models. Acc. Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef] [PubMed]
  42. Irina Massova, P.A.K. Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding. Perspect. Drug Discov. Des. 2000, 18, 113–135. [Google Scholar] [CrossRef]
  43. Sitkoff, D.; Sharp, K.A.; Honig, B. Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models. J. Phys. Chem. 1994, 98, 1978–1988. [Google Scholar] [CrossRef]
  44. Wang, W.; Kollman, P.A. Computational study of protein specificity: The molecular basis of HIV-1 protease drug resistance. Proc. Natl. Acad. Sci. USA 2001, 98, 14937–14942. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, J.; Morin, P.; Wang, W.; Kollman, P.A. Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA. J. Am. Chem. Soc. 2001, 123, 5221–5230. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 00111 g001 550
Figure 1. (A) The time dependence of RMSD of inhibitors (upper) and CK2α (lower) for CK2 in complex with CAM4066, Pre-CAM4066 and compound 21; and (B) calculated per-residue B-factor of three complexes systems.
Figure 1. (A) The time dependence of RMSD of inhibitors (upper) and CK2α (lower) for CK2 in complex with CAM4066, Pre-CAM4066 and compound 21; and (B) calculated per-residue B-factor of three complexes systems.
Ijms 19 00111 g001
Ijms 19 00111 g002 550
Figure 2. (A) Superimposition of co-crystallized pose (gray stick) and the average structure (white stick) of CAM4066; and (B) interactions between CAM4066 and the residues in the active site.
Figure 2. (A) Superimposition of co-crystallized pose (gray stick) and the average structure (white stick) of CAM4066; and (B) interactions between CAM4066 and the residues in the active site.
Ijms 19 00111 g002
Ijms 19 00111 g003 550
Figure 3. (A) Superimposition of co-crystallized pose (gray stick) and the average structure (white stick) ofcompound 21; (B) stable binding mode of compound 21 with the residues in the active site; (C) superimposition of co-crystallized pose (gray stick) and the average structure (white stick) of compound Pre-CAM4066; and (D) stable binding mode of compound Pre-CAM4066 with the residues in the active site.
Figure 3. (A) Superimposition of co-crystallized pose (gray stick) and the average structure (white stick) ofcompound 21; (B) stable binding mode of compound 21 with the residues in the active site; (C) superimposition of co-crystallized pose (gray stick) and the average structure (white stick) of compound Pre-CAM4066; and (D) stable binding mode of compound Pre-CAM4066 with the residues in the active site.
Ijms 19 00111 g003
Ijms 19 00111 g004 550
Figure 4. The dihedral N1-C8-C9-C10 and C8-C9-C10-N2 as the time evolution of CAM4066 and compound 21.
Figure 4. The dihedral N1-C8-C9-C10 and C8-C9-C10-N2 as the time evolution of CAM4066 and compound 21.
Ijms 19 00111 g004
Table
Table 1. Energy terms of MM/PBSA results for three CK2α-inhibitor complexes systems.
Table 1. Energy terms of MM/PBSA results for three CK2α-inhibitor complexes systems.
Energy Term (kcal/mol)CAM4066Compound 21Pre-CAM4066
ΔEele−110.41 ± 2.83−53.42 ± 3.87−11.17 ± 4.26
ΔEvdw−75.00 ± 3.16−65.18 ± 3.64−68.35 ± 4.09
ΔEgas a−185.41 ± 5.23−118.60 ± 3.64−79.52 ± 4.52
ΔGnonpolar−8.56 ± 0.13−8.29 ± 0.19−8.89 ± 0.18
ΔGpolar139.29 ± 4.1875.46 ± 4.4240.24 ± 2.94
ΔGsol b130.72 ± 3.1567.16 ± 3.3631.35 ± 2.96
ΔGele c28.88 ± 5.3222.04 ± 3.1529.07 ± 3.98
ΔGbinding d−54.68 ± 4.87−51.43 ± 4.42−48.17 ± 3.33
ΔΔGbinding03.256.51
a ΔEgas = ΔEele + ΔEvdw; b ΔGsol = ΔGpolar + ΔGnonpolar; c ΔGele = ΔEele + ΔGpolar; d ΔGbinding = ΔEele + ΔEvdw + ΔGsol.
Table
Table 2. Chemical structures, Kd and IC50 values of CAM4066, Pre-CAM4066 and compound 21.
Table 2. Chemical structures, Kd and IC50 values of CAM4066, Pre-CAM4066 and compound 21.
CompoundFragment ALinkerFragment BIC50 (μM)KdPDB
CAM4066 Ijms 19 00111 i001 Ijms 19 00111 i002 Ijms 19 00111 i0030.3700.3205CU4
Pre-CAM4066 Ijms 19 00111 i004 Ijms 19 00111 i005 Ijms 19 00111 i006n/an/an/d
21 Ijms 19 00111 i007 Ijms 19 00111 i008 Ijms 19 00111 i009n/a1.645MO8
n/a = not active; n/d = not determined.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Zhang, N.; Qi, X.; Tang, S.; Sun, G.; Zhao, L.; Zhong, R.; Peng, Y. Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies. Int. J. Mol. Sci. 2018, 19, 111. https://doi.org/10.3390/ijms19010111

AMA Style

Zhou Y, Zhang N, Qi X, Tang S, Sun G, Zhao L, Zhong R, Peng Y. Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies. International Journal of Molecular Sciences. 2018; 19(1):111. https://doi.org/10.3390/ijms19010111

Chicago/Turabian Style

Zhou, Yue, Na Zhang, Xiaoqian Qi, Shan Tang, Guohui Sun, Lijiao Zhao, Rugang Zhong, and Yongzhen Peng. 2018. "Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies" International Journal of Molecular Sciences 19, no. 1: 111. https://doi.org/10.3390/ijms19010111

APA Style

Zhou, Y., Zhang, N., Qi, X., Tang, S., Sun, G., Zhao, L., Zhong, R., & Peng, Y. (2018). Insights into the Impact of Linker Flexibility and Fragment Ionization on the Design of CK2 Allosteric Inhibitors: Comparative Molecular Dynamics Simulation Studies. International Journal of Molecular Sciences, 19(1), 111. https://doi.org/10.3390/ijms19010111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop