Inhibitors of PIM-1 kinase: a computational analysis of the binding free energies of a range of imidazo [1,2-b] pyridazines

Dissociation Pathways of the p53 DNA Binding Domain from DNA and Critical Roles of Key Residues Elucidated by dPaCS-MD/MSM

p53 is a transcriptional factor that regulates cell response to a variety of stresses. About a half of all human tumors contain p53 mutations, and the accumulation of mutations in the DNA binding domain of p53 (p53-DBD) can cause destabilization of p53 and its complex with DNA. To identify the key residues of the p53-DBD/DNA binding and to understand the dissociation mechanisms of the p53-DBD/DNA complex, the dissociation process of p53-DBD from a DNA duplex that contains the consensus sequence (the specific target of p53-DBD) was investigated by a combination of dissociation parallel cascade selection molecular dynamics (dPaCS-MD) and the Markov state model (MSM). This combination (dPaCS-MD/MSM) enabled us to simulate dissociation of the two large molecules based on an all-atom model with a short simulation time (11.2 ± 2.2 ns per trial) and to analyze dissociation pathways, free energy landscape (FEL), and binding free energy. Among 75 trials of dPaCS-MD, p53-DBD dissociated first from the major groove and then detached from the minor groove in 93% of the cases, while 7% of the cases unbinding from the minor groove occurred first. Minor groove binding is mainly stabilized by R248, identified as the most important residue that tightly binds deep inside the minor groove. The standard binding free energy calculated from the FEL was -10.9 ± 0.4 kcal/mol, which agrees with an experimental value of -11.1 kcal/mol. These results indicate that the dPaCS-MD/MSM combination can be a powerful tool to investigate dissociation mechanisms of two large molecules. Analysis of the p53 key residues for DNA binding indicates high correlations with cancer-related mutations, confirming that impairment of the interactions between p53-DBD and DNA can be frequently related to cancer.

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

We recently proposed a computational procedure to simulate the dissociation of protein/ligand complexes using the dissociation Parallel Cascade Selection Molecular Dynamics simulation (dPaCS-MD) method and to analyze the generated trajectories using the Markov state model (MSM). This procedure, called dPaCS-MD/MSM, enables calculation of the dissociation free energy profile and the standard binding free energy. To examine whether this method can reproduce experimentally determined binding free energies for a variety of systems, we used it to investigate the dissociation of three protein/ligand complexes: trypsin/benzamine, FKBP/FK506, and adenosine A_2A receptor/T4E. First, dPaCS-MD generated multiple dissociation pathways within a reasonable computational time for all the complexes, although the complexes differed significantly in the size of the molecules and in intermolecular interactions. Subsequent MSM analyses produced free energy profiles for the dissociations, which provided insights into how each ligand dissociates from the protein. The standard binding free energies obtained by dPaCS-MD/MSM are in good agreement with experimental values for all the complexes. We conclude that dPaCS-MD/MSM can accurately calculate the binding free energies of these complexes.

A systematic review on active sites and functions of PIM-1 protein

The Proviral Integration of Molony murine leukemia virus (PIM)-1 protein contributes to the solid cancers and hematologic malignancies, cell growth, proliferation, differentiation, migration, and other life activities. Many studies have related these functions to its molecular structure, subcellular localization and expression level. However, recognition of specific active sites and their effects on the activity of this constitutively active kinase is still a challenge. Based on the close relationship between its molecular structure and functional activity, this review covers the specific residues involved in the binding of ATP and different substrates in its catalytic domain. This review then elaborates on the relevant changes in protein conformation and cell functions after PIM-1 binds to different substrates. Therefore, this intensive study can improve the understanding of PIM-1-regulated signaling pathways by facilitating the discovery of its potential phosphorylation substrates.

Therapeutic targeting of PIM KINASE signaling in cancer therapy: Structural and clinical prospects

BACKGROUND

PIM kinases are well-studied drug targets for cancer, belonging to Serine/Threonine kinases family. They are the downstream target of various signaling pathways, and their up/down-regulation affects various physiological processes. PIM family comprises three isoforms, namely, PIM-1, PIM-2, and PIM-3, on alternative initiation of translation and they have different levels of expression in different types of cancers. Its structure shows a unique ATP-binding site in the hinge region which makes it unique among other kinases.

SCOPE OF REVIEW

PIM kinases are widely reported in hematological malignancies along with prostate and breast cancers. Currently, many drugs are used as inhibitors of PIM kinases. In this review, we highlighted the physiological significance of PIM kinases in the context of disease progression and therapeutic targeting. We comprehensively reviewed the PIM kinases in terms of their expression and regulation of different physiological roles. We further predicted functional partners of PIM kinases to elucidate their role in the cellular physiology of different cancer and mapped their interaction network.

MAJOR CONCLUSIONS

A deeper mechanistic insight into the PIM signaling involved in regulating different cellular processes, including transcription, apoptosis, cell cycle regulation, cell proliferation, cell migration and senescence, is provided. Furthermore, structural features of PIM have been dissected to understand the mechanism of inhibition and subsequent implication of designed inhibitors towards therapeutic management of prostate, breast and other cancers.

GENERAL SIGNIFICANCE

Being a potential drug target for cancer therapy, available drugs and PIM inhibitors at different stages of clinical trials are discussed in detail.

MovableType Software for Fast Free Energy-Based Virtual Screening: Protocol Development, Deployment, Validation, and Assessment

For decades, the complicated energy surfaces found in macromolecular protein:ligand structures, which require large amounts of computational time and resources for energy state sampling, have been an inherent obstacle to fast, routine free energy estimation in industrial drug discovery efforts. Beginning in 2013, the Merz research group addressed this cost with the introduction of a novel sampling methodology termed “Movable Type” (MT). Using numerical integration methods, the MT method reduces the computational expense for energy state sampling by independently calculating each atomic partition function from an initial molecular conformation in order to estimate the molecular free energy using ensembles of the atomic partition functions. In this work, we report a software package, the DivCon Discovery Suite with the MovableType module from QuantumBio Inc., that performs this MT free energy estimation protocol in a fast, fully encapsulated manner. We discuss the computational procedures and improvements to the original work, and we detail the corresponding settings for this software package. Finally, we introduce two validation benchmarks to evaluate the overall robustness of the method against a broad range of protein:ligand structural cases. With these publicly available benchmarks, we show that the method can use a variety of input types and parameters and exhibits comparable predictability whether the method is presented with “expensive” X-ray structures or “inexpensively docked” theoretical models. We also explore some next steps for the method. The MovableType software is available at http://www.quantumbioinc.com/

Binding affinity and dissociation pathway predictions for a series of USP7 inhibitors with pyrimidinone scaffold by multiple computational methods

Ubiquitin specific protease 7 (USP7) has attracted increasing attention because of its multifaceted roles in different tumor types.

Understanding PIM-1 kinase inhibitor interactions with free energy simulation

The proviral integration site of the Moloney leukemia virus (PIM) family includes three homologous members. PIM-1 kinase is an important target in effective therapeutic interventions of lymphomas, prostate cancer and leukemia. In the current work, we performed free energy calculations to calculate the binding affinities of several inhibitors targeting this protein. The alchemical method with integration and perturbation-based estimators and the end-point methods were compared. The computational results indicated that the alchemical method can accurately predict the binding affinities, while the end-point methods give relatively unreliable predictions. Decomposing the free energy difference into enthalpic and entropic components with MBAR reweighting enabled us to investigate the detailed thermodynamic parameters with which the entropy-enthalpy compensation in this protein-ligand binding case is identified. We then studied the conformational ensemble, and the important protein-ligand interactions were identified. The current work sheds light on the understanding of the PIM-1-kinase-inhibitor interactions at the atomic level and will be useful in the further development of potential drugs.

Dissociation Process of a MDM2/p53 Complex Investigated by Parallel Cascade Selection Molecular Dynamics and the Markov State Model

Recently, we efficiently generated dissociation pathways of a protein-ligand complex without applying force bias with parallel cascade selection molecular dynamics (PaCS-MD) and showed that PaCS-MD in combination with the Markov state model (MSM) yielded a binding free energy comparable to experimental values. In this work, we applied the same procedure to a complex of MDM2 protein and the transactivation domain of p53 protein (TAD-p53), the latter of which is known to be very flexible in the unbound state. Using 30 independent MD simulations in PaCS-MD, we successfully generated 25 dissociation pathways of the complex, which showed complete or partial unfolding of the helical region of TAD-p53 during the dissociation process within an average simulation time of 154.8 ± 46.4 ns. The standard binding free energy obtained in combination with one-dimensional-, three-dimensional (3D)- or C_α-MSM was in good agreement with those determined experimentally. Using 3D-MSM based on the center of mass position of TAD-p53 relative to MDM2, the dissociation rate constant was calculated, which was comparable to those measured experimentally. C_α-MSM based on all C_α coordinates of TAD-p53 reproduced the experimentally measured standard binding free energy, and dissociation and association rate constants. We conclude that the combination of PaCS-MD and MSM offers an efficient computational procedure to calculate binding free energies and kinetic rates.

Structural analysis of PIM1 kinase complexes with ATP-competitive inhibitors

PIM1 is an oncogenic kinase overexpressed in a number of cancers where it correlates with poor prognosis. Several studies demonstrated that inhibition of PIM1 activity is an attractive strategy in fighting overexpressing cancers, while distinct structural features of ATP binding pocket make PIM1 an inviting target for the design of selective inhibitors. To facilitate development of specific PIM1 inhibitors, in this study we report three crystal structures of ATP-competitive inhibitors at the ATP binding pocket of PIM1. Two of the reported structures (CX-4945 and Ro-3306) explain the off-target effect on PIM1 of respectively casein kinase 2 and cyclin-dependent kinase 1 dedicated inhibitors. In turn, the structure with CX-6258 demonstrates a binding mode of a potent, selective inhibitor of PIM1, PIM2, PIM3 and Flt-3 kinases. The consequences of our findings for future inhibitor development are discussed.

A combination strategy to inhibit Pim-1: synergism between noncompetitive and ATP-competitive inhibitors

Pim-1 is a serine/threonine kinase critically involved in the initiation and progression of various types of cancer, especially leukemia, lymphomas and solid tumors such as prostate, pancreas and colon, and is considered a potential drug target against these malignancies. In an effort to discover new potent Pim-1 inhibitors, a previously identified ATP-competitive indolyl-pyrrolone scaffold was expanded to derive structure-activity relationship data. A virtual screening campaign was also performed, which led to the discovery of additional ATP-competitive inhibitors as well as a series of 2-aminothiazole derivatives, which are noncompetitive with respect to both ATP and peptide substrate. This mechanism of action, which resembles allosteric inhibition, has not previously been characterized for Pim-1. Notably, further evaluation of the 2-aminothiazoles indicated a synergistic inhibitory effect in enzymatic assays when tested in combination with ATP-competitive inhibitors. A synergistic effect in the inhibition of cell proliferation by ATP-competitive and ATP-noncompetitive compounds was also observed in prostate cancer cell lines (PC3), where all Pim-1 inhibitors tested in showed synergism with the known anticancer agent, paclitaxel. These results further establish Pim-1 as a target in cancer therapy, and highlight the potential of these agents for use as adjuvant agents in the treatment of cancer diseases in which Pim-1 is associated with chemotherapeutic resistance.

Effect of single amino acid substitution observed in cancer on Pim-1 kinase thermodynamic stability and structure

Pim-1 kinase, a serine/threonine protein kinase encoded by the pim proto-oncogene, is involved in several signalling pathways such as the regulation of cell cycle progression and apoptosis. Many cancer types show high expression levels of Pim kinases and particularly Pim-1 has been linked to the initiation and progression of the malignant phenotype. In several cancer tissues somatic Pim-1 mutants have been identified. These natural variants are nonsynonymous single nucleotide polymorphisms, variations of a single nucleotide occurring in the coding region and leading to amino acid substitutions. In this study we investigated the effect of amino acid substitution on the structural stability and on the activity of Pim-1 kinase. We expressed and purified some of the mutants of Pim-1 kinase that are expressed in cancer tissues and reported in the single nucleotide polymorphisms database. The point mutations in the variants significantly affect the conformation of the native state of Pim-1. All the mutants, expressed as soluble recombinant proteins, show a decreased thermal and thermodynamic stability and a lower activation energy values for kinase activity. The decreased stability accompanied by an increased flexibility suggests that Pim-1 variants may be involved in a wider network of protein interactions. All mutants bound ATP and ATP mimetic inhibitors with comparable IC50 values suggesting that the studied Pim-1 kinase mutants can be efficiently targeted with inhibitors developed for the wild type protein.

Flexibility of the P-loop of Pim-1 kinase: observation of a novel conformation induced by interaction with an inhibitor

The serine/threonine kinase Pim-1 is emerging as a promising target for cancer therapeutics. Much attention has recently been focused on identifying potential Pim-1 inhibitor candidates for the treatment of haematopoietic malignancies. The outcome of a rational drug-design project has recently been reported [Nakano et al. (2012), J. Med. Chem. 55, 5151-5156]. The report described the process of optimization of the structure-activity relationship and detailed from a medicinal chemistry perspective the development of a low-potency and nonselective compound initially identified from in silico screening into a potent, selective and metabolically stable Pim-1 inhibitor. Here, the structures of the initial in silico hits are reported and the noteworthy features of the Pim-1 complex structures are described. A particular focus was placed on the rearrangement of the glycine-rich P-loop region that was observed for one of the initial compounds, (Z)-7-(azepan-1-ylmethyl)-2-[(1H-indol-3-yl)methylidene]-6-hydroxy-1-benzofuran-3(2H)-one (compound 1), and was also found in all further derivatives. This novel P-loop conformation, which appears to be stabilized by an additional interaction with the β3 strand located above the binding site, is not usually observed in Pim-1 structures.

Protein structure along the order-disorder continuum

Thermal fluctuations cause proteins to adopt an ensemble of conformations wherein the relative stability of the different ensemble members is determined by the topography of the underlying energy landscape. “Folded” proteins have relatively homogeneous ensembles, while “unfolded” proteins have heterogeneous ensembles. Hence, the labels “folded” and “unfolded” represent attempts to provide a qualitative characterization of the extent of structural heterogeneity within the underlying ensemble. In this work, we introduce an information-theoretic order parameter to quantify this conformational heterogeneity. We demonstrate that this order parameter can be estimated in a straightforward manner from an ensemble and is applicable to both unfolded and folded proteins. In addition, a simple formula for approximating the order parameter directly from crystallographic B factors is presented. By applying these metrics to a large sample of proteins, we show that proteins span the full range of the order–disorder axis.

Specific CLK inhibitors from a novel chemotype for regulation of alternative splicing

There is a growing recognition of the importance of protein kinases in the control of alternative splicing. To define the underlying regulatory mechanisms, highly selective inhibitors are needed. Here, we report the discovery and characterization of the dichloroindolyl enaminonitrile KH-CB19, a potent and highly specific inhibitor of the CDC2-like kinase isoforms 1 and 4 (CLK1/CLK4). Cocrystal structures of KH-CB19 with CLK1 and CLK3 revealed a non-ATP mimetic binding mode, conformational changes in helix αC and the phosphate binding loop and halogen bonding to the kinase hinge region. KH-CB19 effectively suppressed phosphorylation of SR (serine/arginine) proteins in cells, consistent with its expected mechanism of action. Chemical inhibition of CLK1/CLK4 generated a unique pattern of splicing factor dephosphorylation and had at low nM concentration a profound effect on splicing of the two tissue factor isoforms flTF (full-length TF) and asHTF (alternatively spliced human TF).