Next Article in Journal
Lowering the Intraocular Pressure in Rats and Rabbits by Cordyceps cicadae Extract and Its Active Compounds
Next Article in Special Issue
A New Series of Aryloxyacetic Acids Endowed with Multi-Target Activity towards Peroxisome Proliferator-Activated Receptors (PPARs), Fatty Acid Amide Hydrolase (FAAH), and Acetylcholinesterase (AChE)
Previous Article in Journal
Melatonin: Regulation of Prion Protein Phase Separation in Cancer Multidrug Resistance
Previous Article in Special Issue
Molecular Docking of SP40 Peptide towards Cellular Receptors for Enterovirus 71 (EV-A71)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations

Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543, Singapore
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(3), 706; https://doi.org/10.3390/molecules27030706
Submission received: 23 December 2021 / Revised: 16 January 2022 / Accepted: 19 January 2022 / Published: 21 January 2022
(This article belongs to the Special Issue Small Molecules in Drug Discovery and Pharmacology)

Abstract

:
Haspin, an atypical serine/threonine protein kinase, is a potential target for cancer therapy. 5-iodotubercidin (5-iTU), an adenosine derivative, has been identified as a potent Haspin inhibitor in vitro. In this paper, quantum chemical calculations and molecular dynamics (MD) simulations were employed to identify and quantitatively confirm the presence of halogen bonding (XB), specifically halogen∙∙∙π (aromatic) interaction between halogenated tubercidin ligands with Haspin. Consistent with previous theoretical finding, the site specificity of the XB binding over the ortho-carbon is identified in all cases. A systematic increase of the interaction energy down Group 17, based on both quantum chemical and MD results, supports the important role of halogen bonding in this series of inhibitors. The observed trend is consistent with the experimental observation of the trend of activity within the halogenated tubercidin ligands (F < Cl < Br < I). Furthermore, non-covalent interaction (NCI) plots show that cooperative non-covalent interactions, namely, hydrogen and halogen bonds, contribute to the binding of tubercidin ligands toward Haspin. The understanding of the role of halogen bonding interaction in the ligand–protein complexes may shed light on rational design of potent ligands in the future.

1. Introduction

Haspin, an atypical serine/threonine protein kinase, phosphorylates histone H3 at ‘Thr-3′ during mitosis [1,2,3]. The mRNA of Haspin was first detected and named cell-specific gene 2 (GSG2) in 1994, which was subsequently renamed as haploid germ cell-specific nuclear protein kinase (Haspin) [4]. Significantly different from other protein kinases, phosphorylation is not required for Haspin to be activated. During mitosis, a docking site is created for chromosome passenger complex (CPC), which plays a crucial role to prevent chromosome misalignment [5,6,7,8,9,10,11]. In addition to that, it is also involved in centromeric cohesion and mitotic spindle stability, making it a potential target for cancer therapy [12,13,14,15]. Human Haspin consists of 798 amino acids. The N-terminal part is a less conserved regulatory domain, while the C-terminal is a well-conserved catalytic kinase domain [16,17]. Recently, Chaikuad, Knapp and co-workers [18] identified 5-iodotubercidin (5-iTU), an adenosine derivative, as a potent Haspin inhibitor in vitro, showing the self-renewal and differentiation effect of mouse embryonic stem cells (ESCs) [19,20,21]. Figure 1a [18] shows the noncovalent interactions between tubercidin and Haspin key residues, including Gly491, Val498, Phe605, Glu606, Phe607, Gly608 and Asp611. 5-iTU inhibits Haspin with an IC50 ranged between 5 and9 nM [18,22,23,24]. The authors identified that potency and slow dissociation were increased with the increasing halogen size of 5-iTU derivatives (F < Cl < Br < I) (Figure 1b), which suggests the presence of the halogen bonding, specifically the halogen–aromatic π interactions.
Halogen bonding (XB) [25], driven by anisotropic charge distributions along the extension of the C–X bonds, is quantum-chemical in origin, with an equatorial ring of negative charge and a hind region of positive charge, termed as the σ-hole [26,27]. Halogens are shown to be useful to optimize ADMET properties and prolong the lifetime of the drug [28,29]. For protein–ligand interactions, halogen bonds can be formed between a halogenated ligand and accessible side chain groups, such as hydroxyls, carboxylates, sulfurs, nitrogen and π systems [30,31,32].
Based on an updated search on the Cambridge Crystallographic Database (CSD), we showed that π-type XB represents the majority (66%) of XB contacts identified [33]. For the aromatic type of XB acceptors, we observed a strong preference (95%) of XB binding towards the rim (i.e., over-bond or over-atom) of an aromatic ring, not towards the centroid. Quantum chemical studies of Cl2 and Br2 XB complexes of 14 polycyclic aromatic compounds, including aromatic amino acids, confirmed the observed rim specificity [33]. Intriguingly, the site specificity of the XB binding sites is identified in all cases. The authors demonstrated that the simple frontier orbital interaction readily rationalizes the rim and site specificities of XB involving aromatic XB-acceptors. In particular, the molecular orbital theory provides a proper description of the important charge transfer contribution in XB formation, supported by various energy decomposition analyses [33,34]. These theoretical studies demonstrate the prevalent role of π-type (particularly aromatic-π) XB acceptors in halogen bonded systems.
Non-covalent interactions play an important role in drug design. As one of the key noncovalent interactions, halogen bonding could contribute significantly to the ligand binding affinity and biological properties. In the present paper, we employed computational chemistry methods, both molecular dynamics (MD) simulations and quantum chemical calculations, to study the ligand–protein interaction between halogenated tubercidin derivatives and Haspin. Specifically, we investigated the role of the halogen aromatic π bond in the drug–protein binding.

2. Computational Methods

In this paper, Molecular Operating Environment (MOE) software and Amber14:EHT force field were used for molecular dynamics simulations [35,36]. The Amber14:EHT force field, which is an all-atomic combination of Extended Hückel Theory, is parameterized for non-bonded interactions. To benchmark the performance of the Amber14:EHT force field in prediction of halogen bonding in protein–ligand complexes, thyroid hormone receptor with ligands were simulated. The Berendsen thermostat was applied to control the simulation temperature [37]. The protocol of this MD simulation was started with 100 ps heat from 10 K to 300 K, followed by 1000 ps of NVT and 1000 ps NPT, at 300 K and 100 kPa, followed by a 100 ps of equilibration under constant temperature (300 K). The production was 8 ns at 300 K and 100 kPa, with a time step of 0.002 fs.
For the MD simulations of halogenated tubercidin ligands with a Haspin receptor, the native all-atomic crystal structures, namely, PDB 6G34 (5-iTU), 6G35 (5-brTU), 6G36 (5-clTU) and 6G37 (5-fTU), were taken as starting geometries. The geometries were optimized before proceeding to MD simulation. Hydrogen atoms were added, and partial charges were calculated based on the AMBER14:EHT forcefield. The partial charge for each atom is stored in the MOE internal atom data structure. Energy minimization was carried out using a succession of three algorithms, namely, steepest descent, conjugate gradient and truncated Newton. The protocol of simulation is the same as those of the benchmarking cases. However, the production run was longer at 25 ns at 300 K and 100 kPa, with a time step of 0.004 fs. Binding of the free energy of ligand–receptor complex was calculated using GBVI/WSA method [38] in MOE.
The counterpoise-corrected interaction energies of modeled halogenated ligands-receptor complexes were calculated using density functional theory (DFT) calculations based on the ωB97X-D functional [39]. The Haspin receptor was truncated to the key residue Phe605 only and the halogenated tubercidin ligands were modeled by halogenated pyrrolo[2,3-d]pyrimidine. The aug-cc-pVTZ basis set is used for non-iodine atoms and Def2-TZVPD basis set for iodine atom. It is important to note the ωB97X-D functional benchmarks well against high-level CCSD(T) method in binding energies of XB complexes involving aromatic acceptors [33,40]. Solvent effect was incorporated using the polarizable continuum model [41] to account for the non-specific (macroscopic) effect of the protein dielectric medium. A dielectric constant (ε) of 6 was used as the average dielectric constant inside protein is relatively low (about 6−7) [42,43]. Visualization of noncovalent interactions in ligand–receptor complexes was carried out using the NCI plot [44,45]. The NCI isosurfaces were visualized with VMD program [46] using data produced by Multiwfn program [47]. The strength of the noncovalent interaction is indicated by the color of the isosurface: green represents attractive while blue denotes strongly attractive. The NCI analysis has been successfully used to shed light on the presence of halogen bonding in various chemical systems [48,49,50,51].

3. Results and Discussion

3.1. Benchmark of AMBER14EHT Force Field on XB Ligand–Protein Complexes

The performance of Molecular Operating Environment (MOE) was initially benchmarked with the thyroid hormone receptor complexed with brominated and iodinated ligands, PDB 2J4A and 1XZX, respectively [52,53]. The halogen bonds were identified between halogen atoms (Br and I) and Phe272 carbonyl oxygen atoms in the ligand–protein complexes. The distances for Br∙∙∙O and I∙∙∙O XBs in the crystal structures are 3.28 Å and 3.23 Å and the corresponding XB angles for Br∙∙∙OC and I∙∙∙OC are 163.6° and 165.9°, respectively. The MD simulations reproduced the halogen bond interaction presented in the crystal structures, with XB distances of 3.54 Å and 3.37 Å, and angles of 158.9° and 162.5°, respectively, for the Br∙∙∙O and I∙∙∙O halogen bonds (Figure 2). However, the XB distances were underestimated slightly by 8% and 4%, for the Br and I systems, respectively. These differences will be used to correct the XB distances in the following MD simulations of the halogenated tubercidin ligands.

3.2. MD Simulations of Halogenated Tubercidin Ligands with Haspin

Molecular dynamics simulation provides a dynamic model which can provide insight on ligand-receptor binding by providing information of internal motions and conformational change. In addition, the relative binding free energy predicted by MD simulation may provide a quantitative estimate of drug–protein interaction. The key halogen bond interaction is reproduced in the MD simulations based on the AMBER14:EHT force field (Figure 3), RMSD plot for 5-iTU-haspin MD simulation was provided in supporting information (Figure S1). Consistent with previous theoretical findings on rim and site specificity [33], the halogen∙∙∙π (aromatic) interaction between the halogen atom and the phenyl ring of gate keeper Phe605 is site specifically at the rim of the phenyl ring and is ortho-directed. The ortho-directed effect can be rationalized in terms of the charge polarization effect of a substituted alkyl group, such as toluene [33]. The X∙∙∙π distances were measure between the halogen atom and the closest carbon of the phenyl ring of Phe605 are 3.78, 3.68, 3.57 and 3.52 Å, for 5-fTU, 5-clTU, 5-brTU and 5-iTU, respectively (Figure 3). It was observed that the closest carbon atoms are always the ortho carbons on the aromatic ring of Phe605.
On the basis of the benchmark MD simulations (Section 3.1), the XB distances are corrected for the significant underestimation. For the I, Br and Cl derivatives, the corrected X∙∙∙π interaction distance is less than the sum of van der Waal radii (Table 1). This indicates the presence of a halogen bond between the halogen atom and the aromatic ring (Phe605). In addition, trajectory plots of XB distance and angle (Figure 4) show that 5-iTU ligand stays tightly to the Haspin binding pocket throughout the simulation. This provides further evidence of a strong halogen bonding interaction between 5-iTU and the receptor. The estimated binding free energies, based on the GBVI/WSA_dG method are −46.1, −47.9, −48.6 and −48.9 kcal/mol, for F, Cl, Br and I derivatives, respectively. The GBVI/WSA_dG method is a forcefield based scoring function, which has been trained on 99 protein–ligand complexes dataset [54]. The values may be overestimated for ligand–protein complex, which do not fall in the dataset. However, the trend of binding energies of the halogenated tubercidin derivatives follows the trend of halogen bonding suggests the role of a halogen bond, specifically the halogen–aromatic π interaction. It is important to note that in addition to halogen bonds, various hydrogen bonds, e.g., between ligand NH and Glu606 carbonyl O atoms (see Figure 3), contribute to the overall binding energy listed in Table 1. Hence, it is instructive to perform DFT calculations (Section 3.3) to further assess the strength of halogen bonding quantitatively.

3.3. Quantum Chemical Calculations of Halogen Bond Interaction Energies

Careful examination of the X-ray structures of the ligand–receptor complexes reveals that the XB binding site is over the carbon atom of phenyl ring, not the centroid. Figure 5 depicts the space-filling model of the complex between XB donor 5-iTU and the aromatic side chain Phe605 and the adjacent Glu606 residue, which clearly indicates that over-atom binding mode at Cα (Phe605). The shortest halogen(X)–π distances over Cα of the phenyl ring in the crystal structures are 3.68, 3.59 and 3.52 Å for 5-clTU, 5-brTU and 5-iTU, respectively (Figure 6). The X∙∙∙π interaction distance is less than the sum of van der Waal radii for Br and I derivatives. The rim and site specificity of the X∙∙∙π (aromatic) interaction of the ligand–receptor complexes are readily confirmed by the MD simulations (Section 3.2).
To further shed light on the role of halogen bonding in this series of halogenated tubercidin ligands, DFT calculations were performed on truncated modeled systems of the four halogenated tubercidin derivatives (Figure 6) to quantify the interaction energy (DE) between the halogen containing moiety (XB donor) and the aromatic XB acceptor. In these modeled systems, only Phe605 key residue was included in the receptor and the halogenated tubercidin ligands were modeled by halogenated pyrrolo[2,3-d]pyrimidine. The moiety containing the NH2 group of tubercidin ligand was not included so that the hydrogen bond interaction is excluded in the interaction energy calculation. The modeled structures (see Figure 6) were derived from the crystal structures (6G34, 6G35, 6G36 and 6G37) without further optimization. We observed a systematic increase in the interaction energy (in vacuo) down Group 17, F (−0.72) < Cl (−1.32) < Br (−1.85) < I (−2.66 kJ/mol) (Figure 6). To stimulate the protein environment, implicit PCM solvation model was employed to examine the dielectric effect (ε = 6). The predicted medium effect on the interaction energies is small (~0.1 kcal/mol) and the trend of interaction energy remains the same (Figure 6). These trends of interaction energies readily support the importance role of halogen bonding in this series of inhibitors. The observed trend is consistent with the experimental observation of the trend of activity within the halogenated tubercidin ligands [18].

3.4. Analysis of Noncovalent Interactions

Finally, we employed the noncovalent interaction (NCI) plots to visualize the key noncovalent interactions between 5-iTU and the Haspin key residues. The noncovalent interactions index enables the visualization of noncovalent interactions through the transformation of the reduced electron density gradient into the surface, in which the color is representative of the nature (repulsive or attractive) and strength of the NCI [44,45]. The NCl plot of 5-iTU (Figure 7) reveals two types of attractive intermolecular NCIs, namely, the hydrogen and halogen bonds. The dark blue surface between the N–H proton of 5-iTU and the carbonyl oxygen of Glu606 moiety demonstrates the stronger N–H···O hydrogen bond. The presence of halogen bond is readily reflected in the green disc between the iodine atom of 5-iTU and the phenyl ring of Phe605. The existence of XB in the Br and Cl analogues and the non-existence of XB in the F analogue were also confirmed in the NCI analysis (see Supporting Information, Figure S2). In summary, cooperative noncovalent interactions, namely, hydrogen and halogen bonds, occur in the interaction between 5-iTU and Phe605-Glu606.

4. Conclusions

Discovered two decades ago, the atypical kinase Haspin plays an essential role in Histone H3 phosphorylation and as well as in CPC activity regulation, making itself an attractive target for cancer therapy. On the basis of QC calculations and MD simulations, we have confirmed the presence of halogen bonding between halogenated tubercidin ligands (Cl, Br and I) with the Haspin receptor. Site specificity of the halogen∙∙∙π (aromatic) interaction over the ortho-carbon (Phe605) is observed in the MD simulations. The trend of DFT calculated interaction energy (F < Cl < Br < I) supports the important role of halogen bonding in this series of halogenated inhibitors. With the understanding of the role of halogen bonding interaction, we hope this type of non-covalent interaction can be further exploited in the rational design of therapeutic drugs in the future.

Supplementary Materials

The following supporting information can be downloaded, Figure S1. Root mean square deviation (RMSD) plot of complex between 5-iTU and Haspin; Figure S2. NCI isosurface between the tubercidin ligand 5-XTU and Haspin key residues (Phe605-Glu606): (a) 5-brTU, (b) 5-clTU, and (c) 5-fTU.

Author Contributions

Conceptualization, M.W.W.; formal analysis and investigation, Y.Z. and M.W.W.; writing—original draft preparation, review and editing, Y.Z. and M.W.W.; funding acquisition, M.W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education (Singapore), MOE2016-T2-1-087.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

No sample is available from the authors.

References

  1. Tanaka, H.; Yoshimura, Y.; Nishina, Y.; Nozaki, M.; Nojima, H.; Nishimune, Y. Isolation and characterization of cDNA clones specifically expressed in testicular germ cells. FEBS Lett. 1994, 355, 4–10. [Google Scholar] [CrossRef] [Green Version]
  2. Tanaka, H.; Yoshimura, Y.; Nozaki, M.; Yomogida, K.; Tsuchida, J.; Tosaka, Y.; Habu, T.; Nakanishi, T.; Okada, M.; Nojima, H.; et al. Identification and characterization of a haploid germ cell-specific nuclear protein kinase (Haspin) in spermatid nuclei and its effects on somatic cells. J. Biol.Chem. 1999, 274, 17049–17057. [Google Scholar] [CrossRef] [Green Version]
  3. Polioudaki, H.; Markaki, Y.; Kourmouli, N.; Dialynas, G.; Theodoropoulos, P.A.; Singh, P.B.; Georgatos, S.D. Mitotic phosphorylation of histone H3 at threonine 3. FEBS Lett. 2004, 560, 39–44. [Google Scholar] [CrossRef] [Green Version]
  4. Higgins, J.M.G. Haspin-like proteins: A new family of evolutionarily conserved putative eukaryotic protein kinases. Protein Sci. 2001, 10, 1677–1684. [Google Scholar] [CrossRef]
  5. Hindriksen, S.; Bramer, A.J.; Truong, M.A.; Vromans, M.J.M.; Post, J.B.; Verlaan-Klink, I.; Snippert, H.J.; Lens, S.M.A.; Hadders, M.A. Baculoviral delivery of CRISPR/Cas9 facilitates efficient genome editing in human cells. PLoS ONE 2017, 12, e0179514. [Google Scholar] [CrossRef]
  6. Higgins, J.M.G. Structure, function and evolution of haspin and haspin-related proteins, a distinctive group of eukaryotic protein kinases. Cell Mol. Life Sci. 2003, 60, 446–462. [Google Scholar] [CrossRef]
  7. Dai, J.; Sultan, S.; Taylor, S.S.; Higgins, J.M.G. The kinase haspin is required for mitotic histone H3 Thr 3 phosphorylation and normal metaphase chromosome alignment. Genes Dev. 2005, 19, 472–488. [Google Scholar] [CrossRef] [Green Version]
  8. Markaki, Y.; Christogianni, A.; Politou, A.S.; Georgatos, S.D. Phosphorylation of histone H3 at Thr3 is part of a combinatorial pattern that marks and configures mitotic chromatin. J. Cell Sci. 2009, 122, 2809–2819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kelly, A.E.; Ghenoiu, C.; Xue, J.Z.; Zierhut, C.; Kimura, H.; Funabiki, H. Survivin reads phosphorylated histone H3 threonine 3 to activate the mitotic kinase Aurora B. Science 2010, 330, 235–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Wang, F.; Dai, J.; Daum, J.R.; Niedzialkowska, E.; Banerjee, B.; Stukenberg, P.T.; Gorbsky, G.J.; Higgins, J.M.G. Histone H3 Thr-3 phosphorylation by haspin positions Aurora B at centromeres in mitosis. Science 2010, 330, 231–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Du, J.; Kelly, A.E.; Funabiki, H.; Patel, D.J. Structural basis for recognition of H3T3ph and Smac/DIABLO N-terminal peptides by human survivin. Structure 2012, 20, 185–195. [Google Scholar] [CrossRef] [Green Version]
  12. Biggins, S.; Murray, A.W. The budding yeast protein kinase Ipl1/Aurora allows the absence of tension to activate the spindle checkpoint. Genes Dev. 2001, 15, 3118–3129. [Google Scholar] [CrossRef] [Green Version]
  13. Tanaka, H.; Iguchi, N.; Nakamura, Y.; Kohroki, J.; de Carvalho, C.E.; Nishimune, Y. Cloning and characterization of human haspin gene encoding haploid germ cell-specific nuclear protein kinase. Mol. Hum. Reprod. 2001, 7, 211–218. [Google Scholar] [CrossRef] [Green Version]
  14. Hauf, S.; Cole, R.W.; LaTerra, S.; Zimmer, C.; Schnapp, G.; Walter, R.; Heckel, A.; van Meel, J.; Rieder, C.L.; Peters, J.-M. The small molecule Hesperadin reveals a role for Aurora B in correcting kinetochore-microtubule attachment and in maintaining the spindle assembly checkpoint. J. Cell Biol. 2003, 161, 281–294. [Google Scholar] [CrossRef] [PubMed]
  15. Pinsky, B.A.; Kung, C.; Shokat, K.M.; Biggins, S. The Ipl1-Aurora protein kinase activates the spindle checkpoint by creating unattached kinetochores. Nat. Cell Biol. 2006, 8, 78–83. [Google Scholar] [CrossRef] [PubMed]
  16. Villa, F.; Capasso, P.; Tortorici, M.; Forneris, F.; de Marco, A.; Mattevi, A.; Musacchio. Crystal structure of the catalytic domain of Haspin, an atypical kinase implicated in chromatin organization. Proc. Natl. Acad. Sci. USA 2009, 106, 20204–20209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Eswaran, J.; Patnaik, D.; Filippakopoulos, P.; Wang, F.; Stein, R.L.; Murray, J.W.; Higgins, J.M.G.; Knapp, S. Structure and functional characterization of the atypical human kinase haspin. Proc. Natl. Acad. Sci. USA 2009, 106, 20198–20203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Heroven, C.; Georgi, V.; Ganotra, G.K.; Brennan, P.; Wolfreys, F.; Wade, R.C.; Fernandez-Montalvan, A.E.; Chaikuad, A.; Knapp, S. Halogen-aromatic π interactions modulate inhibitor residence times. Angew. Chem. Int. Ed. Engl. 2018, 57, 7220–7224. [Google Scholar] [CrossRef]
  19. Newby, A.C. The role of adenosine kinase in regulating adenosine concentration. Biochem. J. 1985, 226, 343–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Fedorov, O.; Marsden, B.; Pogacic, V.; Rellos, P.; Muller, S.; Bullock, A.N. A systematic interaction map of validated kinase inhibitors with Ser/Thr kinases. Proc. Natl. Acad. Sci. USA 2007, 104, 20523–20528. [Google Scholar] [CrossRef] [Green Version]
  21. Balzano, D.; Santaguida, S.; Musacchio, A.; Villa, F. A general framework for inhibitor resistance in protein kinases. Chem. Biol. 2011, 18, 966–975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. De Antoni, A.; Maffini, S.; Knapp, S.; Musacchio, A.; Santaguida, S. A small-molecule inhibitor of Haspin alters the kinetochore functions of Aurora B. J. Cell Biol. 2012, 199, 269–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, F.; Ulyanova, N.P.; Daum, J.R.; Patnaik, D.; Kateneva, A.V.; Gorbsky, G.J.; Higgins, J.M.G. Haspin inhibitors reveal centromeric functions of Aurora B in chromosome segregation. J. Cell Biol. 2012, 199, 251–268. [Google Scholar] [CrossRef] [Green Version]
  24. Kozgunova, E.; Suzuki, T.; Ito, M.; Higashiyama, T.; Kurihara, D. Haspin has multiple functions in the plant cell division regulatory network. Plant Cell Physiol. 2016, 57, 848–861. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The halogen bond. Chem. Rev. 2016, 116, 2478–2601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Clark, T.; Hennemann, M.; Murray, J.S.; Politzer, P. Halogen bonding: The σ-hole. Proceedings of “Modeling interactions in biomolecules II”, Prague, September 5th–9th, 2005. J. Mol. Model. 2007, 13, 291–296. [Google Scholar] [CrossRef]
  27. Politzer, P.; Murray, J.S.; Clark, T. Halogen bonding: An electrostatically-driven highly directional noncovalent interaction. Phys. Chem. Chem. Phys. 2010, 12, 7748–7757. [Google Scholar] [CrossRef] [PubMed]
  28. Wilcken, R.; Zimmermann, M.O.; Lange, A.; Joerger, A.C.; Boeckler, F.M. Principles and applications of halogen bonding in medicinal chemistry and chemical biology. J. Med. Chem. 2013, 56, 1363–1388. [Google Scholar] [CrossRef]
  29. Xu, Z.; Yang, Z.; Liu, Y.; Lu, Y. Chen, K.; Zhu, W. Halogen Bond: Its Role beyond Drug−Target Binding Affinity for Drug Discovery and Development. J. Chem. Inf. Model. 2014, 54, 69–78. [Google Scholar] [CrossRef]
  30. Carter, M.; Rappé, A.K.; Ho, P.S. Scalable Anisotropic Shape and Electrostatic Models for Biological Bromine Halogen Bond. J. Chem. Theory Comput. 2012, 8, 2461–2473. [Google Scholar] [CrossRef]
  31. Scholfield, M.R.; Vander Zanden, C.M.; Carter, M.; Ho, P.S. Halogen bonding (X-bonding): A biological perspective. Prot. Sci. 2013, 22, 139–152. [Google Scholar] [CrossRef]
  32. Parisini, E.; Metrangolo, P.; Pilati, T.; Resnati, G.; Terraneo, G. Halogen bonding in halocarbon-protein complexes: A structural survey. Chem. Soc. Rev. 2011, 40, 2267–2278. [Google Scholar] [CrossRef] [PubMed]
  33. Ang, S.J.; Mak, A.M.; Sullivan, M.B.; Wong, M.W. Site specificity of halogen bonding involving aromatic acceptors. Phys. Chem. Chem. Phys. 2018, 20, 8685–8694. [Google Scholar] [CrossRef] [PubMed]
  34. Ang, S.J.; Mak, A.M.; Wong, M.W. Nature of halogen bonding involving π-systems, nitroxide radicals and carbenes: A highlight of the importance of charge transfer. Phys. Chem. Chem. Phys. 2018, 20, 26463–26478. [Google Scholar] [CrossRef] [PubMed]
  35. Chemical Computing Group ULC. Molecular Operating Environment (MOE); Chemical Computing Group ULC: Montreal, QC, Canada, 2020. [Google Scholar]
  36. Labute, P.; Williams, C. Application of Hückel Theory to Pharmacophore Discovery. CICSJ Bull. 2015, 33, 33–40. [Google Scholar]
  37. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular-Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  38. Labute, P. The generalized Born/volume integral implicit solvent model: Estimation of the free energy of hydration using London dispersion instead of atomic surface area. J. Comput. Chem. 2008, 29, 1693–1698. [Google Scholar] [CrossRef]
  39. Chai, J.-D.; Head-Gordon, D. Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 128, 6615–6620. [Google Scholar] [CrossRef] [Green Version]
  40. Kozuch, S.; Martin, J.M.L. Halogen Bonds: Benchmarks and Theoretical Analysis. J. Chem. Theory Comput. 2013, 9, 1918–1931. [Google Scholar] [CrossRef]
  41. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105, 2999–3093. [Google Scholar] [CrossRef]
  42. Li, L.; Li, C.; Zhang, Z.; Alexov, E. On the Dielectric “Constant” of Proteins: Smooth Dielectric Function for Macromolecular Modeling and Its Implementation in DelPhi. J. Chem. Theory Comput. 2013, 9, 2126–2136. [Google Scholar] [CrossRef]
  43. Siegbahn, P.E.M.; Himo, F. The quantum chemical cluster approach for modeling enzyme reactions. WIREs Comput. Mol. Sci. 2011, 1, 323–336. [Google Scholar] [CrossRef]
  44. Johnson, E.R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A.J.; Yang, W. Revealing Noncovalent Interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506. [Google Scholar] [CrossRef] [Green Version]
  45. Contreras-García, J.; Johnson, E.R.; Keinan, S.; Chaudret, R.; Piquemal, J.-P.; Beratan, D.N.; Yang, W. NCIPLOT: A program for plotting noncovalent interaction regions. J. Chem. Theory Comput. 2011, 7, 625–632. [Google Scholar] [CrossRef] [PubMed]
  46. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  47. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  48. Kee, C.W.; Wong, M.W. In silico design of halogen-bonding-based organocatalyst for Diels–Alder reaction, Claisen rearrangement, and Cope-type hydroamination. J. Org. Chem. 2016, 81, 7459–7470. [Google Scholar] [CrossRef]
  49. Narth, C.; Maroun, Z.; Boto, R.; Chaudret, C.; Bonnet, M.L.; Piquemal, J.-P.; Contreras-Garcia, J. A Complete NCI Perspective: From New Bonds to Reactivity. In Challenges and Advances in Computational Chemistry and Physics; Springer Nature: Basingstoke, UK, 2016; Chapter 18. [Google Scholar]
  50. Song, S.; Wang, L.; Su, J.; Xu, Z.; Hsu, C.-H.; Lyu, P.; Li, J.; Peng, X.; Takahiro, K.; Telychko, M.; et al. Manifold Dynamic Non-Covalent Interactions for Steering Molecular Assembly and Cyclization. Chem. Sci. 2021, 12, 11659–11667. [Google Scholar] [CrossRef]
  51. Masaryk, L.; Moncol, J.; Herchel, R.; Nemec, I. Halogen Bonding in New Dichloride-Cobalt(II) Complex with Iodo Substituted Chalcone Ligands. Crystals 2020, 10, 354. [Google Scholar] [CrossRef]
  52. Koehler, K.; Gordon, S.; Brandt, P.; Carlsson, B.; Backsbro-Saeidi, A.; Apelqvist, T.; Agback, P.; Grover, G.J.; Nelson, W.; Grynfarb, M.; et al. Thyroid Receptor Ligands. 6. A High Affinity “Direct Antagonist” Selective for the Thyroid Hormone Receptor. J. Med. Chem. 2006, 49, 6635–6637. [Google Scholar] [CrossRef]
  53. Sandler, B.; Webb, P.; Apriletti, J.W.; Huber, B.R.; Togashi, M.; Lima, S.T.C.; Juric, S.; Nilsson, S.; Wagner, T.; Fletterick, R.J.; et al. Thyroxine-Thyroid Hormone Receptor Interactions. J. Biol. Chem. 2004, 279, 55801–55808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Corbeil, C.R.; Williams, C.I.; Labute, P. Variability in docking success rates due to dataset preparation. J. Comput. Aided. Mol. Des. 2012, 26, 775–786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Interaction site of tubercidin inhibitors with Haspin key residues (PDB 6G35). (b) Molecular structures of halogenated tubercidin derivatives.
Figure 1. (a) Interaction site of tubercidin inhibitors with Haspin key residues (PDB 6G35). (b) Molecular structures of halogenated tubercidin derivatives.
Molecules 27 00706 g001
Figure 2. (a) Superimpose of MD predicted binding pose (blue) with crystal structure 2J4A (red). (b) Superimpose of MD predicted binding pose (blue) with crystal structure 1XZX (Red).
Figure 2. (a) Superimpose of MD predicted binding pose (blue) with crystal structure 2J4A (red). (b) Superimpose of MD predicted binding pose (blue) with crystal structure 1XZX (Red).
Molecules 27 00706 g002
Figure 3. Key interaction sites of tubercidin inhibitors with Haspin from MD simulations. Interaction distance (Å) between halogen atom and Cα carbon of Phe605 phenyl ring: (a) 5-fTU, (b) 5-clTU, (c) 5-brTU and (d) 5-iTU. The geometries were chosen by the shortest distance from the last 50 snapshots for each MD trajectory.
Figure 3. Key interaction sites of tubercidin inhibitors with Haspin from MD simulations. Interaction distance (Å) between halogen atom and Cα carbon of Phe605 phenyl ring: (a) 5-fTU, (b) 5-clTU, (c) 5-brTU and (d) 5-iTU. The geometries were chosen by the shortest distance from the last 50 snapshots for each MD trajectory.
Molecules 27 00706 g003
Figure 4. Plot of XB parameters of MD simulation trajectory of 5-iTU with Haspin receptor: (A) XB distance (I∙∙∙Cα) and (B) XB angle (<C-I∙∙∙Cα). The RMSD plot is given in Supporting Information (Figure S1).
Figure 4. Plot of XB parameters of MD simulation trajectory of 5-iTU with Haspin receptor: (A) XB distance (I∙∙∙Cα) and (B) XB angle (<C-I∙∙∙Cα). The RMSD plot is given in Supporting Information (Figure S1).
Molecules 27 00706 g004
Figure 5. Space-filling model of close contact between 5-iTU and truncated receptor residues (Phe605 and Glu606) of PDB 6G34.
Figure 5. Space-filling model of close contact between 5-iTU and truncated receptor residues (Phe605 and Glu606) of PDB 6G34.
Molecules 27 00706 g005
Figure 6. Structures and interaction energies (ΔE, kcal/mol) of modeled halogen bonded complexes. Calculated interaction energies in a dielectric medium of ε = 6 are given in parenthesis.
Figure 6. Structures and interaction energies (ΔE, kcal/mol) of modeled halogen bonded complexes. Calculated interaction energies in a dielectric medium of ε = 6 are given in parenthesis.
Molecules 27 00706 g006
Figure 7. NCI isosurface between the tubercidin ligand 5-iTU and Haspin key residues (Phe605-Glu606) of PDB 6G34. Interaction strength increases from green to blue.
Figure 7. NCI isosurface between the tubercidin ligand 5-iTU and Haspin key residues (Phe605-Glu606) of PDB 6G34. Interaction strength increases from green to blue.
Molecules 27 00706 g007
Table 1. Halogen bond distances and binding energies of halogenated tubercidin ligands with Haspin receptor obtained from MD simulations.
Table 1. Halogen bond distances and binding energies of halogenated tubercidin ligands with Haspin receptor obtained from MD simulations.
LigandsX∙∙∙π Closest Carbon Distance (Å) aSum of VDW Radii (Å) bBinding Energy (kcal/moL) c
5-iTu3.52 (3.38)3.68−48.9
5-brTu3.57 (3.28)3.55−48.6
5-clTu3.68 (3.39)3.45−47.9
5-fTu3.78 (3.48)3.17−46.1
(a) Corrected XB distances, in parenthesis, based on benchmark comparison. (b) Van der Waals radii of C, F, Cl, Br and I are 1.70, 1.47, 1.75, 1.85 and 1.98 Å, respectively. (c) Binding free energy derived from molecular dynamic simulations.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Wong, M.W. Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations. Molecules 2022, 27, 706. https://doi.org/10.3390/molecules27030706

AMA Style

Zhou Y, Wong MW. Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations. Molecules. 2022; 27(3):706. https://doi.org/10.3390/molecules27030706

Chicago/Turabian Style

Zhou, Yujing, and Ming Wah Wong. 2022. "Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations" Molecules 27, no. 3: 706. https://doi.org/10.3390/molecules27030706

APA Style

Zhou, Y., & Wong, M. W. (2022). Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations. Molecules, 27(3), 706. https://doi.org/10.3390/molecules27030706

Article Metrics

Back to TopTop