Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Small Conformational Changes Underlie Evolution of Resistance to NNRTI in HIV Reverse Transcriptase

Despite achieving considerable success in reducing the number of fatalities due to acquired immunodeficiency syndrome, emergence of resistance against the reverse transcriptase (RT) inhibitor drugs remains one of the biggest challenges of the human immunodeficiency virus antiretroviral therapy (ART). Non-nucleoside reverse transcriptase inhibitors (NNRTIs) form a large class of drugs and a crucial component of ART. In NNRTIs, even a single resistance mutation is known to make the drugs completely ineffective. Additionally, several inhibitor-bound RTs with single resistance mutations do not exhibit any significant variations in their three-dimensional structures compared with the inhibitor-bound RT but completely nullify their inhibitory functions. This makes understanding the structural mechanism of these resistance mutations crucial for drug development. Here, we study several single resistance mutations in the allosteric inhibitor (nevirapine)-bound RT to analyze the mechanism of small structural changes leading to these large functional effects. In this study, we have shown that in absence of significant conformational variations in the inhibitor-bound wild-type RT and RT with single resistance mutations, the protein contact network analysis of their static structures, along with molecular dynamics simulations, can be a useful approach to understand the functional effect of small local conformational variations. The simple network analysis exposes the localized contact changes that lead to global rearrangement in the communication pattern within RT. Furthermore, these conformational changes have implications on the overall dynamics of RT. Using various measures, we show that a single resistance mutation can change the network structure and dynamics of RT to behave more like unbound RT, even in the presence of the inhibitor. This combined coarse-grained contact network and molecular dynamics approach promises to be a useful tool to analyze structure-function studies of proteins that show large functional changes with negligible variations in their overall conformation.

Tios: The Internet of Simulations. Turning Molecular Dynamics into a Data Streaming Web Application

The configuration of most current academic high-performance computing (HPC) resources tends to enforce ways of working with, and thinking about, molecular dynamics (MD) simulations that are not always optimal. For example, when the aim of the simulation(s) is to produce a representative sample of a Boltzmann weighted ensemble, the ideal scenario would be to be able to do just that-i.e. to tap into a running simulation of indefinite length, collect data from it in real time, and only terminate the simulation once the quality of a sample was assured. Current approaches, based on batch jobs of proscribed maximum length, and a postprocessing style of data analysis, inhibit this. In the spirit of the Internet of Things, we have developed Tios, a Python application that turns MD simulations into remotely discoverable and accessible streaming web applications to which researchers can connect and download data as they please. Tios is freely available, works with standard MD codes, and requires no modifications to them. In this paper we outline how Tios works and present a number of test cases that demonstrate its capabilities.

Diacritic Binding of an Indenoindole Inhibitor by CK2α Paralogs Explored by a Reliable Path to Atomic Resolution CK2α' Structures

CK2α and CK2α' are the two isoforms of the catalytic subunit of human protein kinase CK2, an important target for cancer therapy. They have similar, albeit not identical functional and structural properties, and were occasionally reported to be inhibited with distinct efficacies by certain ATP-competitive ligands. Here, we present THN27, an indeno[1,2-b]indole derivative, as a further inhibitor with basal isoform selectivity. The selectivity disappears when measured using CK2α/CK2α' complexes with CK2β, the regulatory CK2 subunit. Co-crystal structures of THN27 with CK2α and CK2α' reveal that subtle differences in the conformational variability of the interdomain hinge region are correlated with the observed effect. In the case of CK2α', a crystallographically problematic protein so far, this comparative structural analysis required the development of an experimental strategy that finally enables atomic resolution structure determinations with ab initio phasing of potentially any ATP-competitive CK2 inhibitor and possibly many non-ATP-competitive ligands as well bound to CK2α'.

The mechanism of PI3Kα activation at the atomic level

PI3K lipid kinases phosphorylate PIP₂ to PIP₃ in the PI3K/Akt/mTOR pathway to regulate cellular processes. They are frequently mutated in cancer. Here we determine the PI3Kα activation mechanism at the atomic level. Unlike protein kinases where the substrate abuts the ATP, crystal structures indicate that in PI3Kα, the distance between the γ phosphate of the ATP and the PIP₂ lipid substrate is over 6 Å, much too far for the phosphoryl transfer, raising the question of how catalysis is executed. PI3Kα has two subunits, the catalytic p110α and the regulatory p85α. Our simulations show that release of the autoinhibition exerted by the nSH2 domain of the p85α triggers significant conformational change in p110α, leading to the exposure of the kinase domain for membrane interaction. Structural rearrangement in the C-lobe of the kinase domain reduces the distance between the ATP γ-phosphate and the substrate, offering an explanation as to how phosphoryl transfer is executed. An alternative mechanism may involve ATP relocation. This mechanism not only explains how oncogenic mutations promote PI3Kα activation by facilitating nSH2 release, or nSH2-release-induced, allosteric motions; it also offers an innovative, PI3K isoform-specific drug discovery principle. Rather than competing with nanomolar range ATP in the ATP-binding pocket and contending with ATP pocket conservation and massive binding targets, this mechanism suggests blocking the PI3Kα sequence-specific cavity between the ATP-binding pocket and the substrate binding site. Targeting isoform-specific residues in the cavity may prevent PIP₂ phosphorylation.

2-Aminothiazole Derivatives as Selective Allosteric Modulators of the Protein Kinase CK2. 1. Identification of an Allosteric Binding Site

CK2 is a ubiquitous Ser/Thr protein kinase involved in the control of various signaling pathways and is known to be constitutively active. In the present study, we identified aryl 2-aminothiazoles as a novel class of CK2 inhibitors, which displayed a non-ATP-competitive mode of action and stabilized an inactive conformation of CK2 in solution. Enzyme kinetics studies, STD NMR, circular dichroism spectroscopy, and native mass spectrometry experiments demonstrated that the compounds bind in an allosteric pocket outside the ATP-binding site. Our data, combined with molecular docking studies, strongly suggested that this new binding site was located at the interface between the αC helix and the flexible glycine-rich loop. A first hit optimization led to compound 7, exhibiting an IC₅₀ of 3.4 μM against purified CK2α in combination with a favorable selectivity profile. Thus, we identified a novel class of CK2 inhibitors targeting an allosteric pocket, offering great potential for further optimization into anticancer drugs.

2-Aminothiazole Derivatives as Selective Allosteric Modulators of the Protein Kinase CK2. 2. Structure-Based Optimization and Investigation of Effects Specific to the Allosteric Mode of Action

Protein CK2 has gained much interest as an anticancer drug target in the past decade. We had previously described the identification of a new allosteric site on the catalytic α-subunit, along with first small molecule ligands based on the 4-(4-phenylthiazol-2-ylamino)benzoic acid scaffold. In the present work, structure optimizations guided by a binding model led to the identification of the lead compound 2-hydroxy-4-((4-(naphthalen-2-yl)thiazol-2-yl)amino)benzoic acid (27), showing a submicromolar potency against purified CK2α (IC₅₀ = 0.6 μM). Furthermore, 27 induced apoptosis and cell death in 786-O renal cell carcinoma cells (EC₅₀ = 5 μM) and inhibited STAT3 activation even more potently than the ATP-competitive drug candidate CX-4945 (EC₅₀ of 1.6 μM vs 5.3 μM). Notably, the potencies of our allosteric ligands to inhibit CK2 varied depending on the individual substrate. Altogether, the novel allosteric pocket was proved a druggable site, offering an excellent perspective to develop efficient and selective allosteric CK2 inhibitors.

Cell-based screen identifies a new potent and highly selective CK2 inhibitor for modulation of circadian rhythms and cancer cell growth

Compounds targeting the circadian clock have been identified as potential treatments for clock-related diseases, including cancer. Our cell-based phenotypic screen revealed uncharacterized clock-modulating compounds. Through affinity-based target deconvolution, we identified GO289, which strongly lengthened circadian period, as a potent and selective inhibitor of CK2. Phosphoproteomics identified multiple phosphorylation sites inhibited by GO289 on clock proteins, including PER2 S693. Furthermore, GO289 exhibited cell type-dependent inhibition of cancer cell growth that correlated with cellular clock function. The x-ray crystal structure of the CK2α-GO289 complex revealed critical interactions between GO289 and CK2-specific residues and no direct interaction of GO289 with the hinge region that is highly conserved among kinases. The discovery of GO289 provides a direct link between the circadian clock and cancer regulation and reveals unique design principles underlying kinase selectivity.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.