Structural basis for decreased affinity of Emodin binding to Val66-mutated human CK2 alpha as determined by molecular dynamics

Molecular modelling methods in food safety: Bisphenols as case study

Bisphenol A (BPA), a synthetic compound widely used as a building block for polycarbonate plastics, has been declared in the European Union (EU) as a substance of very high concern (SVHC). A series of BPA alternatives and derivatives (bisphenols/BPs) with similar physical-chemical properties have been produced and used by companies for substituting it. To evaluate the estrogenic and androgenic binding activity of 26 BPs, a non-statistical in silico approach has been applied. The results of molecular docking analyses applied on six different nuclear receptors (NRs) have revealed that: i) some BPA metabolites could lower the harmful effects of BPA exposure; ii) BPS is a lower interactor for all NRs, but it does not appear safer at all for androgen receptor (AR), for which its binding activity is found similar to a pharmacological anti-androgen; iii) only a BP has been found as a safer compound for all NRs considered. Moreover, molecular dynamic simulation of three BPs on ERα have revealed that the presence of negative hydrophobic interactions could induce a decrease in receptor activity. Overall, the present results demonstrate that in silico methods could be a valid approach to screen estrogenic and androgenic activity of food contact materials (FCMs).

The protein kinase CK2 catalytic domain from Plasmodium falciparum: crystal structure, tyrosine kinase activity and inhibition

Malaria causes every year over half-a-million deaths. The emergence of parasites resistant to available treatments makes the identification of new targets and their inhibitors an urgent task for the development of novel anti-malaria drugs. Protein kinase CK2 is an evolutionary-conserved eukaryotic serine/threonine protein kinase that in Plasmodium falciparum (PfCK2) has been characterized as a promising target for chemotherapeutic intervention against malaria. Here we report a crystallographic structure of the catalytic domain of PfCK2α (D179S inactive single mutant) in complex with ATP at a resolution of 3.0 Å. Compared to the human enzyme, the structure reveals a subtly altered ATP binding pocket comprising five substitutions in the vicinity of the adenine base, that together with potential allosteric sites, could be exploited to design novel inhibitors specifically targeting the Plasmodium enzyme. We provide evidence for the dual autophosphorylation of residues Thr⁶³ and Tyr³⁰ of PfCK2. We also show that CX4945, a human CK2 inhibitor in clinical trials against solid tumor cancers, is effective against PfCK2 with an IC₅₀ of 13.2 nM.

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

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.

Exploring the Pivotal Role of the CK2 Hinge Region Sub-Pocket in Binding with Tricyclic Quinolone Analogues by Computational Analysis

Protein kinase CK2 has been considered as an attractive therapeutic target of cancer therapy. The tricyclic quinoline compound CX-4945 is the first representative of CK2 inhibitors used in human clinical trials. The binding of non-2,6-naphtyridine substituted compounds 27e (IC₅₀ > 500 nM) and 27h (IC₅₀ > 1000 nM) to CK2 is abolished. However, the unbinding mechanisms due to the key pharmacophore group replacement of compounds 27e and 27h are unveiled. In the present work, combined computational analysis was performed to investigate the underlying structural basis of the low-affinity of two systems. As indicated in the results, the loss of hydrogen bonds between the non-2,6-naphtyridine and the hinge region destroyed the proper recognition of the two complexes. Besides, the allosteric mechanisms between the deviated ligands and the changed regions (G-loop, C-loop and β4/β5 loop) are proposed. Furthermore, energetic analysis was evaluated by detailed energy calculation and residue-based energy decomposition. More importantly, the summary of known polar pharmacophore groups elucidates the pivotal roles of hinge region sub-pocket in the binding of CK2 inhibitors. These results provide rational clues to the fragment-based design of more potent CK2 inhibitors.

Underlying mechanisms of cyclic peptide inhibitors interrupting the interaction of CK2α/CK2β: comparative molecular dynamics simulation studies

Protein-protein interactions (PPIs) are fundamental to all biological processes. Recently, the CK2β-derived cyclic peptide Pc has been demonstrated to efficiently antagonize the CK2α/CK2β interaction and strongly affect the phosphorylation of CK2β-dependent CK2 substrate specificity. The binding affinity of Pc to CK2α is destroyed to different extents by two single-point mutations of Tyr188 to Ala (Y188A) and Phe190 to Ala (F190A), which exert negative effects on the inhibitory activity (IC50) of Pc against the CK2α/CK2β interaction from 3.0 μM to 54.0 μM and ≫100 μM, respectively. However, the structural influences of Y188A and F190A mutations on the CK2α-Pc complex remain unclear. In this study, comparative molecular dynamics (MD) simulations, principal component analysis (PCA), domain cross-correlation map (DCCM) analysis and energy calculations were performed on wild type (WT), Y188A mutant, and F190A mutant systems. The results revealed that ordered communications between hydrophobic and polar interactions were essential for CK2α-Pc binding in the WT system. In addition to the loss of the hydrogen bond between Gln36 of CK2α and Gly189 of Pc in the two mutants, the improper recognition mechanisms occurred through different pathways. These pathways included the weakened hydrophobic interactions in the Y188A mutant as well as decreased polar and hydrophobic interactions in the F190A mutant. The energy analysis results qualitatively elucidated the instability of the two mutants and energetic contributions of the key residues. This study not only revealed the structural mechanisms for the decreased binding affinity of Y188A and F190A mutant CK2α-Pc complexes, but also provided valuable clues for the rational design of CK2α/CK2β subunit interaction inhibitors with high affinity and specificity.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Bioinformatics and variability in drug response: a protein structural perspective

Marketed drugs frequently perform worse in clinical practice than in the clinical trials on which their approval is based. Many therapeutic compounds are ineffective for a large subpopulation of patients to whom they are prescribed; worse, a significant fraction of patients experience adverse effects more severe than anticipated. The unacceptable risk-benefit profile for many drugs mandates a paradigm shift towards personalized medicine. However, prior to adoption of patient-specific approaches, it is useful to understand the molecular details underlying variable drug response among diverse patient populations. Over the past decade, progress in structural genomics led to an explosion of available three-dimensional structures of drug target proteins while efforts in pharmacogenetics offered insights into polymorphisms correlated with differential therapeutic outcomes. Together these advances provide the opportunity to examine how altered protein structures arising from genetic differences affect protein-drug interactions and, ultimately, drug response. In this review, we first summarize structural characteristics of protein targets and common mechanisms of drug interactions. Next, we describe the impact of coding mutations on protein structures and drug response. Finally, we highlight tools for analysing protein structures and protein-drug interactions and discuss their application for understanding altered drug responses associated with protein structural variants.

Unbiased functional proteomics strategy for protein kinase inhibitor validation and identification of bona fide protein kinase substrates: application to identification of EEF1D as a substrate for CK2

Protein kinases have emerged as attractive targets for treatment of several diseases prompting large-scale phosphoproteomics studies to elucidate their cellular actions and the design of novel inhibitory compounds. Current limitations include extensive reliance on consensus predictions to derive kinase-substrate relationships from phosphoproteomics data and incomplete experimental validation of inhibitors. To overcome these limitations in the case of protein kinase CK2, we employed functional proteomics and chemical genetics to enable identification of physiological CK2 substrates and validation of CK2 inhibitors including TBB and derivatives. By 2D electrophoresis and mass spectrometry, we identified the translational elongation factor EEF1D as a protein exhibiting CK2 inhibitor-dependent decreases in phosphorylation in (32)P-labeled HeLa cells. Direct phosphorylation of EEF1D by CK2 was shown by performing CK2 assays with EEF1D -FLAG from HeLa cells. Dramatic increases in EEF1D phosphorylation following λ-phosphatase treatment and phospho- EEF1D antibody recognizing EEF1D pS162 indicated phosphorylation at the CK2 site in cells. Furthermore, phosphorylation of EEF1D in the presence of TBB or TBBz is restored using CK2 inhibitor-resistant mutants. Collectively, our results demonstrate that EEF1D is a bona fide physiological CK2 substrate for CK2 phosphorylation. Furthermore, this validation strategy could be adaptable to other protein kinases and readily combined with other phosphoproteomic methods.