NMR characterization of kinase p38 dynamics in free and ligand-bound forms

Dimerization Promotes PKR Activation by Modulating Energetics of αC Helix Conversion between Active and Inactive Conformations

Protein kinase R (PKR) functions in the eukaryotic innate immune system as a first-line defense against viral infections. PKR binds viral dsRNA, leading to autophosphorylation and activation. In its active state, PKR can phosphorylate its primary substrate, eIF2α, which blocks the initiation of translation in the infected cell. It has been established that PKR activation occurs when the kinase domain dimerizes in a back-to-back configuration. However, the mechanism by which dimerization leads to enzymatic activation is not fully understood. Here, we investigate the structural mechanistic basis and energy landscape for PKR activation, with a focus on the αC helix─a kinase activation and signal integration hub─using all-atom equilibrium and enhanced sampling molecular dynamics simulations. By employing window-exchange umbrella sampling, we compute free-energy profiles of activation, which show that back-to-back dimerization stabilizes a catalytically competent conformation of PKR. Key hydrophobic residues in the homodimer interface contribute to stabilization of the αC helix in an active conformation and the position of its critical glutamate residue. Using linear mutual information analysis, we analyze allosteric communication connecting the protomers' N-lobes and the αC helix dimer interface with the αC helix.

Exploring the conformational landscapes of protein kinases: perspectives from FRET and DEER

Conformational changes of catalytically-important structural elements are a key feature of the regulation mechanisms of protein kinases and are important for dictating inhibitor binding modes and affinities. The lack of widely applicable methods for tracking kinase conformational changes in solution has hindered our understanding of kinase regulation and our ability to design conformationally selective inhibitors. Here we provide an overview of two recently developed methods that detect conformational changes of the regulatory activation loop and αC-helix of kinases and that yield complementary information about allosteric mechanisms. An intramolecular Förster resonance energy transfer-based approach provides a scalable platform for detecting and classifying structural changes in high-throughput, as well as quantifying ligand binding cooperativity, shedding light on the energetics governing allostery. The pulsed electron paramagnetic resonance technique double electron-electron resonance provides lower throughput but higher resolution information on structural changes that allows for unambiguous assignment of conformational states and quantification of population shifts. Together, these methods are shedding new light on kinase regulation and drug interactions and providing new routes for the identification of novel kinase inhibitors and allosteric modulators.

Styrylquinazoline derivatives as ABL inhibitors selective for different DFG orientations

Among tyrosine kinase inhibitors, quinazoline-based compounds represent a large and well-known group of multi-target agents. Our previous studies have shown interesting kinases inhibition activity for a series of 4-aminostyrylquinazolines based on the CP-31398 scaffold. Here, we synthesised a new series of styrylquinazolines with a thioaryl moiety in the C4 position and evaluated in detail their biological activity. Our results showed high inhibition potential against non-receptor tyrosine kinases for several compounds. Molecular docking studies showed differential binding to the DFG conformational states of ABL kinase for two derivatives. The compounds showed sub-micromolar activity against leukaemia. Finally, in-depth cellular studies revealed the full landscape of the mechanism of action of the most active compounds. We conclude that S4-substituted styrylquinazolines can be considered as a promising scaffold for the development of multi-kinase inhibitors targeting a desired binding mode to kinases as effective anticancer drugs.

Molecular Dynamics and Machine Learning Give Insights on the Flexibility-Activity Relationships in Tyrosine Kinome

Tyrosine kinases are a subfamily of kinases with critical roles in cellular machinery. Dysregulation of their active or inactive forms is associated with diseases like cancer. This study aimed to holistically understand their flexibility-activity relationships, focusing on pockets and fluctuations. We studied 43 different tyrosine kinases by collecting 120 μs of molecular dynamics simulations, pocket and residue fluctuation analysis, and a complementary machine learning approach. We found that the inactive forms often have increased flexibility, particularly at the DFG motif level. Noteworthy, thanks to these long simulations combined with a decision tree, we identified a semiquantitative fluctuation threshold of the DGF+3 residue over which the kinase has a higher probability to be in the inactive form.

Tackling Hysteresis in Conformational Sampling: How to Be Forgetful with MEMENTO

The structure of proteins has long been recognized to hold the key to understanding and engineering their function, and rapid advances in structural biology and protein structure prediction are now supplying researchers with an ever-increasing wealth of structural information. Most of the time, however, structures can only be determined in free energy minima, one at a time. While conformational flexibility may thus be inferred from static end-state structures, their interconversion mechanisms─a central ambition of structural biology─are often beyond the scope of direct experimentation. Given the dynamical nature of the processes in question, many studies have attempted to explore conformational transitions using molecular dynamics (MD). However, ensuring proper convergence and reversibility in the predicted transitions is extremely challenging. In particular, a commonly used technique to map out a path from a starting to a target conformation called steered MD (SMD) can suffer from starting-state dependence (hysteresis) when combined with techniques such as umbrella sampling (US) to compute the free energy profile of a transition. Here, we study this problem in detail on conformational changes of increasing complexity. We also present a new, history-independent approach that we term "MEMENTO" (Morphing End states by Modelling Ensembles with iNdependent TOpologies) to generate paths that alleviate hysteresis in the construction of conformational free energy profiles. MEMENTO utilizes template-based structure modelling to restore physically reasonable protein conformations based on coordinate interpolation (morphing) as an ensemble of plausible intermediates, from which a smooth path is picked. We compare SMD and MEMENTO on well-characterized test cases (the toy peptide deca-alanine and the enzyme adenylate kinase) before discussing its use in more complicated systems (the kinase P38α and the bacterial leucine transporter LeuT). Our work shows that for all but the simplest systems SMD paths should not in general be used to seed umbrella sampling or related techniques, unless the paths are validated by consistent results from biased runs in opposite directions. MEMENTO, on the other hand, performs well as a flexible tool to generate intermediate structures for umbrella sampling. We also demonstrate that extended end-state sampling combined with MEMENTO can aid the discovery of collective variables on a case-by-case basis.

Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset of Resistance Mutations in Kinases

Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all-atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti-cancer drug Imatinib (Gleevec) against wild-type and the N368S mutant of Abl kinase. We show how this point mutation triggers far-reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.

Multiscale computational study of ligand binding pathways: Case of p38 MAP kinase and its inhibitors

Protein kinases are one of the most important drug targets in the past 10 years. Understanding the inhibitor association processes will profoundly impact new binder designs with preferred binding kinetics. However, after more than a decade of effort, a complete atomistic-level study of kinase inhibitor binding pathways is still lacking. As all kinases share a similar scaffold, we used p38 kinase as a model system to investigate the conformational dynamics and free energy transition of inhibitor binding toward kinases. Two major kinase conformations, Asp-Phe-Gly (DFG)-in and DFG-out, and three types of inhibitors, type I, II, and III, were thoroughly investigated in this work. We performed Brownian dynamics simulations and up to 340 μs Gaussian-accelerated molecular dynamics simulations to capture the inhibitor binding paths and a series of conformational transitions of the p38 kinase from its apo to inhibitor-bound form. Eighteen successful binding trajectories, including all types of inhibitors, are reported herein. Our simulations suggest a mechanism of inhibitor recruitment, a faster ligand association step to a pre-existing DFG-in/DFG-out p38 protein, followed by a slower molecular rearrangement step to adjust the protein-ligand conformation followed by a shift in the energy landscape to reach the final bound state. The ligand association processes also reflect the energetic favor of type I and type II/III inhibitor binding through ATP and allosteric channels, respectively. These different binding routes are directly responsible for the fast (type I binders) and slow (type II/III binders) kinetics of different types of p38 inhibitors. Our findings also echo the recent study of p38 inhibitor dissociation, implying that ligand unbinding could undergo a reverse path of binding, and both processes share similar metastates. This study deepens the understanding of molecular and energetic features of kinase inhibitor-binding processes and will inspire future drug development from a kinetic point of view.

Automated affinity selection for rapid discovery of peptide binders

In-solution affinity selection (AS) of large synthetic peptide libraries affords identification of binders to protein targets through access to an expanded chemical space. Standard affinity selection methods, however, can be time-consuming, low-throughput, or provide hits that display low selectivity to the target. Here we report an automated bio-layer interferometry (BLI)-assisted affinity selection platform. When coupled with tandem mass spectrometry (MS), this method enables both rapid de novo discovery and affinity maturation of known peptide binders with high selectivity. The BLI-assisted AS-MS technology also features real-time monitoring of the peptide binding during the library selection process, a feature unattainable by current selection approaches. We show the utility of the BLI AS-MS platform toward rapid identification of novel nanomolar (dissociation constant, K_D < 50 nM) non-canonical binders to the leukemia-associated oncogenic protein menin. To our knowledge, this is the first application of BLI to the affinity selection of synthetic peptide libraries. We believe our approach can significantly accelerate the use of synthetic peptidomimetic libraries in drug discovery.

This work reports an automated affinity selection-mass spectrometry (AS-MS) approach amenable to both de novo peptide binder discovery and affinity maturation of known binders in a high-throughput and selective manner.

Binding mechanism of inhibitors to p38α MAP kinase deciphered by using multiple replica Gaussian accelerated molecular dynamics and calculations of binding free energies

The p38α MAP Kinase has been an important target of drug design for treatment of inflammatory diseases and cancers. This work applies multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding mechanism of inhibitors L51, R24 and 1AU to p38α. Dynamics analyses show that inhibitor bindings exert significant effect on conformational changes of the active helix α2 and the conserved DFG loop. The rank of binding free energies calculated with MM-GBSA not only agrees well with that determined by the experimental IC50 values but also suggests that mutual compensation between the enthalpy and entropy changes can improve binding of inhibitors to p38α. The analyses of free energy landscapes indicate that the L51, R24 and 1AU bound p38α display a DFG-out conformation. The residue-based free energy decomposition method is used to evaluate contributions of separate residues to the inhibitor-p38α binding and the results imply that residues V30, V38, L74, L75, I84, T106, H107, L108, M109, L167, F169 and D168 can be utilized as efficient targets of potent inhibitors toward p38α.

Computer Simulations of the Dissociation Mechanism of Gleevec from Abl Kinase with Milestoning

Gleevec (a.k.a., imatinib) is an important anticancer (e.g., chronic myeloid leukemia) chemotherapeutic drug due to its inhibitory interaction with the Abl kinase. Here, we use atomically detailed simulations within the Milestoning framework to study the molecular dissociation mechanism of Gleevec from Abl kinase. We compute the dissociation free energy profile, the mean first passage time for unbinding, and explore the transition state ensemble of conformations. The milestones form a multidimensional network with average connectivity of about 2.93, which is significantly higher than the connectivity for a one-dimensional reaction coordinate. The free energy barrier for Gleevec dissociation is estimated to be ∼10 kcal/mol, and the exit time is ∼55 ms. We examined the transition state conformations using both, the committor and transition function. We show that near the transition state the highly conserved salt bridge K217 and E286 is transiently broken. Together with the calculated free energy profile, these calculations can advance the understanding of the molecular interaction mechanisms between Gleevec and Abl kinase and play a role in future drug design and optimization studies.

Conformational selection vs. induced fit: insights into the binding mechanisms of p38α MAP Kinase inhibitors

The conformational dynamics of a kinase's activation loop have been challenging to assess due to the activation loop's intrinsic flexibility. To directly probe the conformational equilibrium of the activation loop of mitogen-activated protein kinase p38α, we present an approach based on site-directed spin labeling, electron paramagnetic resonance (EPR) distance restraints, and multilateration. We demonstrate that the activation loop of apo p38α resides in a highly flexible equilibrium state and we reveal that binding of small molecules significantly alters this equilibrium and the populated sub-states.

Backbone resonance assignments of the catalytic and regulatory domains of Ca2+/calmodulin-dependent protein kinase 1D

The CaMK subfamily of Ser/Thr kinases are regulated by calmodulin interactions with their C-terminal regions. They are exemplified by Ca²⁺/calmodulin dependent protein kinase 1δ which is known as CaMK1D, CaMKIδ or CKLiK. CaMK1D mediates intracellular signalling downstream of Ca²⁺ influx and thereby exhibits amplifications of Ca²⁺signals and polymorphisms that have been implicated in breast cancer and diabetes. Here we report the backbone ¹H, ¹³C, ¹⁵N assignments of the 38 kDa human CaMK1D protein in its free state, including both the canonical bi-lobed kinase fold as well as the autoinhibitory and calmodulin binding domains.

Revealing the Unbinding Kinetics and Mechanism of Type I and Type II Protein Kinase Inhibitors by Local-Scaled Molecular Dynamics Simulations

Protein kinase inhibitors disrupt phosphorylation of the target kinases, which are an important class of drug for treating cancer and other diseases. Conventional structure-based design methods (such as molecular docking) focus on the static binding mode of the kinase inhibitor with its target. However, dissociation kinetic properties of a drug molecule are found to correlate with its residence time in vivo and thus have drawn the attention of drug designers in recent years. In this study, we have applied the local-scaled molecular dynamics (MD) simulation enabled in GROMACS software to explore the unbinding mechanism of a total of 41 type I and type II kinase inhibitors. Our simulation considered multiple starting configurations as well as possible protonation states of kinase inhibitors. Based on our local-scaled MD results, we discovered that the integrals of the favorable binding energy during dissociation correlated well (R² = 0.64) with the experimental dissociation rate constants of those kinase inhibitors on the entire data set. Given its accuracy and technical advantage, this method may serve as a practical option for estimating this important property in reality. Our simulation also provided a reasonable explanation of the dynamic properties of kinase and its inhibitor as well as the role of relevant water molecules in dissociation.

The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning

BACKGROUND

Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest.

METHOD

A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition.

RESULTS

The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space.

CONCLUSIONS

The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements.

GENERAL SIGNIFICANCE

The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.

p38α Mitogen-Activated Protein Kinase Is a Druggable Target in Pancreatic Adenocarcinoma

p38 mitogen-activated protein kinases are signaling molecules with major involvement in cancer. A detailed mechanistic understanding of how p38 MAPK family members function is urgently warranted for cancer targeted therapy. The conformational dynamics of the most common member of p38 MAPK family, p38α, are crucial for its function but poorly understood. Here we found that, unlike in other cancer types, p38α is significantly activated in pancreatic adenocarcinoma samples, suggesting its potential for anti-pancreatic cancer therapy. Using a state of the art supercomputer, Anton, long-timescale (39 μs) unbiased molecular dynamics simulations of p38α show that apo p38α has high structural flexibility in six regions, and reveal potential catalysis mechanism involving a “butterfly” motion. Moreover, in vitro studies show the low-selectivity of the current p38α inhibitors in both human and mouse pancreatic cancer cell lines, while computational solvent mapping identified 17 novel pockets for drug design. Taken together, our study reveals the conformational dynamics and potentially druggable pockets of p38α, which may potentiate p38α-targeting drug development and benefit pancreatic cancer patients.

How Electrostatic Coupling Enables Conformational Plasticity in a Tyrosine Kinase

Protein kinases are important cellular signaling molecules involved in cancer and a multitude of other diseases. It is well-known that inactive kinases display a remarkable conformational plasticity; however, the molecular mechanisms remain poorly understood. Conformational heterogeneity presents an opportunity but also a challenge in kinase drug discovery. The ability to predictively model various conformational states could accelerate selective inhibitor design. Here we performed a proton-coupled molecular dynamics study to explore the conformational landscape of a c-Src kinase. Starting from a completely inactive structure, the simulations captured all major types of conformational states without the use of a target structure, mutation, or bias. The simulations allowed us to test the experimental hypotheses regarding the mechanism of DFG flip, its coupling to the αC-helix movement, and the formation of regulatory spine. Perhaps the most significant finding is how key titratable residues, such as DFG-Asp, αC-Glu, and HRD-Asp, change protonation states dependent on the DFG, αC, and activation loop conformations. Our data offer direct evidence to support a long-standing hypothesis that protonation of Asp favors the DFG-out state and explain why DFG flip is also possible in simulations with deprotonated Asp. The simulations also revealed intermediate states, among which a unique DFG-out/α-C state formed as DFG-Asp is moved into a back pocket forming a salt bridge with catalytic Lys, which can be tested in selective inhibitor design. Our finding of how proton coupling enables the remarkable conformational plasticity may shift the paradigm of computational studies of kinases which assume fixed protonation states. Understanding proton-coupled conformational dynamics may hold a key to further innovation in kinase drug discovery.

Application of a Substrate-Mediated Selection with c-Src Tyrosine Kinase to a DNA-Encoded Chemical Library

As aberrant activity of protein kinases is observed in many disease states, these enzymes are common targets for therapeutics and detection of activity levels. The development of non-natural protein kinase substrates offers an approach to protein substrate competitive inhibitors, a class of kinase inhibitors with promise for improved specificity. Also, kinase activity detection approaches would benefit from substrates with improved activity and specificity. Here, we apply a substrate-mediated selection to a peptidomimetic DNA-encoded chemical library for enrichment of molecules that can be phosphorylated by the protein tyrosine kinase, c-Src. Several substrates were identified and characterized for activity. A lead compound (SrcDEL10) showed both the ability to serve as a substrate and to promote ATP hydrolysis by the kinase. In inhibition assays, compounds displayed IC_50's ranging from of 8-100 µM. NMR analysis of SrcDEL10 bound to the c-Src:ATP complex was conducted to characterize the binding mode. An ester derivative of the lead compound demonstrated cellular activity with inhibition of Src-dependent signaling in cell culture. Together, the results show the potential for substrate-mediated selections of DNA-encoded libraries to discover molecules with functions other than simple protein binding and offer a new discovery method for development of synthetic tyrosine kinase substrates.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

Potential Mean Force from Umbrella Sampling Simulations: What Can We Learn and What Is Missed?

Changes in free energy provide valuable information for molecular recognition, including both ligand-receptor binding thermodynamics and kinetics. Umbrella sampling (US), a widely used free energy calculation method, has long been used to explore the dissociation process of ligand-receptor systems and compute binding free energy. In existing publications, the binding free energy computed from the potential of mean force (PMF) with US simulation mostly yielded "ball park" values with experimental data. However, the computed PMF values are highly influenced by factors such as initial conformations and/or trajectories provided, the reaction coordinate, and sampling of conformational space in each US window. These critical factors have rarely been carefully studied. Here we used US to study the guest aspirin and 1-butanol dissociation processes of β-cyclodextrin (β-CD) and an inhibitor SB2 dissociation from a p38α mitogen-activated protein kinase (MAPK) complex. For β-CD, we used three different β-CD conformations to generate the dissociation path with US windows. For p38α, we generated the dissociation pathway by using accelerated molecular dynamics followed by conformational relaxing with short conventional MD, steered MD, and manual pulling. We found that, even for small β-CD complexes, different β-CD conformations altered the height of the PMF, but the pattern of PMF was not affected if the MD sampling in each US window was well-converged. Because changing the macrocyclic ring conformation needs to rotate dihedral angles in the ring, a bound ligand largely restrains the motion of cyclodextrin. Therefore, once a guest is in the binding site, cyclodextrin cannot freely change its initial conformation, resulting in different absolute heights of the PMF, which cannot be overcome by running excessively long MD simulations for each US window. Moreover, if the US simulations were not converged, the important barrier and minimum were missed. For ligand-protein systems, our studies also suggest that the dissociation trajectories modeled by an enhanced sampling method must maintain a natural molecular movement to avoid biased PMF plots when using US simulations.

Human p38α mitogen-activated protein kinase in the Asp168-Phe169-Gly170-in (DFG-in) state can bind allosteric inhibitor Doramapimod

Doramapimod (BIRB-796) is widely recognized as one of the most potent and selective type II inhibitors of human p38α mitogen-activated protein kinase (MAPK); however, the understanding of its binding mechanism remains incomplete. Previous studies indicated high affinity of the ligand to a so-called allosteric pocket revealed only in the 'out' state of the DFG motif (i.e. Asp168-Phe169-Gly170) when Phe169 becomes fully exposed to the solvent. The possibility of alternative binding in the DFG-in state was hypothesized, but the molecular mechanism was not known. Methods of bioinformatics, docking and long-time scale classical and accelerated molecular dynamics have been applied to study the interaction of Doramapimod with the human p38α MAPK. It was shown that Doramapimod can bind to the protein even when the Phe169 is fully buried inside the allosteric pocket and the kinase activation loop is in the DFG-in state. Orientation of the inhibitor in such a complex is significantly different from that in the known crystallographic complex formed by the kinase in the DFG-out state; however, the Doramapimod's binding is followed by the ligand-induced conformational changes, which finally improve accommodation of the inhibitor. Molecular modelling has confirmed that Doramapimod combines the features of type I and II inhibitors of p38α MAPK, i.e. can directly and indirectly compete with the ATP binding. It can be concluded that optimization of the initial binding in the DFG-in state and the final accommodation in the DFG-out state should be both considered at designing novel efficient type II inhibitors of MAPK and homologous proteins. Communicated by Ramaswamy H. Sarma.

Kinase Inhibitory Activities and Molecular Docking of a Novel Series of Anticancer Pyrazole Derivatives

A series of novel 1,3,4-triarylpyrazoles containing different heterocycles has been prepared, characterized and screened for their in vitro antiproliferative activity against HePG-2, MCF-7, PC-3, A-549 and HCT-116 cancer cell lines. The biological results revealed that compound 6 showed the highest anticancer activity so it was subjected to a kinase assay study where it reduced the activity of several protein kinases including AKT1, AKT2, BRAF V600E, EGFR, p38α and PDGFRβ at 100 μM using the radiometric or ADP-Glo assay method. Molecular docking simulation supported the initial kinase assay and suggested a common mode of interaction at the ATP-binding sites of these kinases, which demonstrates that compound 6 is a potential agent for cancer therapy deserving further research.

Direct monitoring of the conformational equilibria of the activation loop in the mitogen-activated protein kinase p38α

Conformational transitions in protein kinases are crucial for the biological function of these enzymes. Here, we characterize and assess conformational equilibria of the activation loop and the effect of small molecule inhibitors in the MAP kinase p38α. Our work experimentally revealed the existence of a two-state equilibrium for p38α while the addition of inhibitors shifts the equilibrium between these two states.

Probing Protein-Protein Interactions Using Asymmetric Labeling and Carbonyl-Carbon Selective Heteronuclear NMR Spectroscopy

Protein-protein interactions (PPIs) regulate a plethora of cellular processes and NMR spectroscopy has been a leading technique for characterizing them at the atomic resolution. Technically, however, PPIs characterization has been challenging due to multiple samples required to characterize the hot spots at the protein interface. In this paper, we review our recently developed methods that greatly simplify PPI studies, which minimize the number of samples required to fully characterize residues involved in the protein-protein binding interface. This original strategy combines asymmetric labeling of two binding partners and the carbonyl-carbon label selective (CCLS) pulse sequence element implemented into the heteronuclear single quantum correlation (¹H-15N HSQC) spectra. The CCLS scheme removes signals of the J-coupled 15N⁻13C resonances and records simultaneously two individual amide fingerprints for each binding partner. We show the application to the measurements of chemical shift correlations, residual dipolar couplings (RDCs), and paramagnetic relaxation enhancements (PRE). These experiments open an avenue for further modifications of existing experiments facilitating the NMR analysis of PPIs.

NVP-BHG712: Effects of Regioisomers on the Affinity and Selectivity toward the EPHrin Family

Erythropoietin-producing hepatocellular (EPH) receptors are transmembrane receptor tyrosine kinases. Their extracellular domains bind specifically to ephrin A/B ligands, and this binding modulates intracellular kinase activity. EPHs are key players in bidirectional intercellular signaling, controlling cell morphology, adhesion, and migration. They are increasingly recognized as cancer drug targets. We analyzed the binding of NVP-BHG712 (NVP) to EPHA2 and EPHB4. Unexpectedly, all tested commercially available NVP samples turned out to be a regioisomer (NVPiso) of the inhibitor, initially described in a Novartis patent application. They only differ by the localization of a single methyl group on either one of two adjacent nitrogen atoms. The two compounds of identical mass revealed different binding modes. Furthermore, both in vitro and in vivo experiments showed that the isomers differ in their kinase affinity and selectivity.

Dynamics of human protein kinase Aurora A linked to drug selectivity

Protein kinases are major drug targets, but the development of highly-selective inhibitors has been challenging due to the similarity of their active sites. The observation of distinct structural states of the fully-conserved Asp-Phe-Gly (DFG) loop has put the concept of conformational selection for the DFG-state at the center of kinase drug discovery. Recently, it was shown that Gleevec selectivity for the Tyr-kinase Abl was instead rooted in conformational changes after drug binding. Here, we investigate whether protein dynamics after binding is a more general paradigm for drug selectivity by characterizing the binding of several approved drugs to the Ser/Thr-kinase Aurora A. Using a combination of biophysical techniques, we propose a universal drug-binding mechanism, that rationalizes selectivity, affinity and long on-target residence time for kinase inhibitors. These new concepts, where protein dynamics in the drug-bound state plays the crucial role, can be applied to inhibitor design of targets outside the kinome.

Protein kinases are a family of enzymes found in all living organisms. These enzymes help to control many biological processes, including cell division. When particular protein kinases do not work correctly, cells may start to divide uncontrollably, which can lead to cancer. One example is the kinase Aurora A, which is over-active in many common human cancers. As a result, researchers are currently trying to design drugs that reduce the activity of Aurora A in the hope that these could form new anticancer treatments.

In general, drugs are designed to be as specific in their action as possible to reduce the risk of harmful side effects to the patient. Designing a drug that affects a single protein kinase, however, is difficult because there are hundreds of different kinases in the body, all with similar structures. Because drugs often work by binding to specific structural features, a drug that targets one protein kinase can often alter the activity of a large number of others too.

Gleevec is a successful anti-leukemia drug that specifically works on one target kinase, producing minimal side effects. It was recently discovered that the drug works through a phenomenon called ‘induced fit’. This means that after the drug binds it causes a change in the enzyme’s overall shape that alters the activity of the enzyme. The shape change is complex, and so even small structural differences can change the effect of a particular drug.

Do other drugs that target other protein kinases also produce induced fit effects? To find out, Pitsawong, Buosi, Otten, Agafonov et al. studied how three anti-cancer drugs interact with Aurora A: two drugs specifically designed to switch off Aurora A, and Gleevec (which does not target Aurora A).

The two drugs that specifically target Aurora A were thought to work by targeting one structural feature of the enzyme. However, the biochemical and biophysical experiments performed by Pitsawong et al. revealed that these drugs instead work through an induced fit effect. By contrast, Gleevec did not trigger an induced fit on Aurora A and so bound less tightly to it.

In light of these results, Pitsawong et al. suggest that future efforts to design drugs that target protein kinases should focus on exploiting the induced fit process. This will require more research into the structure of particular kinases.

Role of Molecular Interactions and Protein Rearrangement in the Dissociation Kinetics of p38α MAP Kinase Type-I/II/III Inhibitors

Understanding the governing factors of fast or slow inhibitor binding/unbinding assists in developing drugs with preferred kinetic properties. For inhibitors with the same binding affinity targeting different binding sites of the same protein, the kinetic behavior can profoundly differ. In this study, we investigated unbinding kinetics and mechanisms of fast (type-I) and slow (type-II/III) binders of p38α mitogen-activated protein kinase, where the crystal structures showed that type-I and type-II/III inhibitors bind to pockets with different conformations of the Asp-Phe-Gly (DFG) motif. The work used methods that combine conventional molecular dynamics (MD), accelerated molecular dynamics (AMD) simulations, and the newly developed pathway search guided by internal motions (PSIM) method to find dissociation pathways. The study focuses on revealing key interactions and molecular rearrangements that hinder ligand dissociation by using umbrella sampling and post-MD processing to examine changes in free energy during ligand unbinding. As anticipated, the initial dissociation steps all require breaking interactions that appeared in crystal structures of the bound complexes. Interestingly, for type-I inhibitors such as SB2, p38α keeps barrier-free conformational fluctuation in the ligand-bound complex and during ligand dissociation. In contrast, with a type-II/III inhibitor such as BIRB796, with the rearrangements of p38α in its bound state, ligand unbinding features energetically unfavorable protein-ligand concerted movement. Our results also show that the type-II/III inhibitors preferred dissociation pathways through the allosteric channel, which is consistent with an existing publication. The study suggests that the level of required protein rearrangement is one major determining factor of drug binding kinetics in p38α systems, providing useful information for development of inhibitors.

Predicting the Conformational Variability of Abl Tyrosine Kinase using Molecular Dynamics Simulations and Markov State Models

Understanding protein conformational variability remains a challenge in drug discovery. The issue arises in protein kinases, whose multiple conformational states can affect the binding of small-molecule inhibitors. To overcome this challenge, we propose a comprehensive computational framework based on Markov state models (MSMs). Our framework integrates the information from explicit-solvent molecular dynamics simulations to accurately rank-order the accessible conformational variants of a target protein. We tested the methodology using Abl kinase with a reference and blind-test set. Only half of the Abl conformational variants discovered by our approach are present in the disclosed X-ray structures. The approach successfully identified a protein conformational state not previously observed in public structures but evident in a retrospective analysis of Lilly in-house structures: the X-ray structure of Abl with WHI-P154. Using a MSM-derived model, the free energy landscape and kinetic profile of Abl was analyzed in detail highlighting opportunities for targeting the unique metastable states.

Dynamic activation and regulation of the mitogen-activated protein kinase p38

Mitogen-activated protein kinases, which include p38, are essential for cell differentiation and autophagy. The current model for p38 activation involves activation-loop phosphorylation with subsequent substrate binding leading to substrate phosphorylation. Despite extensive efforts, the molecular mechanism of activation remains unclear. Here, using NMR spectroscopy, we show how the modulation of protein dynamics across timescales activates p38. We find that activation-loop phosphorylation does not change the average conformation of p38; rather it quenches the loop ps-ns dynamics. We then show that substrate binding to nonphosphorylated and phosphorylated p38 results in uniform µs-ms backbone dynamics at catalytically essential regions and across the entire molecule, respectively. Together, these results show that phosphorylation and substrate binding flatten the energy landscape of the protein, making essential elements of allostery and activation dynamically accessible. The high degree of structural conservation among ser/thr kinases suggests that elements of this mechanism may be conserved across the kinase family.

NMR Spectroscopy to Study MAP Kinase Binding to MAP Kinase Phosphatases

NMR spectroscopy and other solution methods are increasingly being used to obtain novel insights into the mechanisms by which MAPK regulatory proteins bind and direct the activity of MAPKs. Here, we describe how interactions between the MAPK p38α and its regulatory proteins are studied using NMR spectroscopy, isothermal titration calorimetry, and small angle X-ray scattering (SAXS).

Dynamic Equilibrium of the Aurora A Kinase Activation Loop Revealed by Single-Molecule Spectroscopy

The conformation of the activation loop (T-loop) of protein kinases underlies enzymatic activity and influences the binding of small-molecule inhibitors. By using single-molecule fluorescence spectroscopy, we have determined that phosphorylated Aurora A kinase is in dynamic equilibrium between a DFG-in-like active T-loop conformation and a DFG-out-like inactive conformation, and have measured the rate constants of interconversion. Addition of the Aurora A activating protein TPX2 shifts the equilibrium towards an active T-loop conformation whereas addition of the inhibitors MLN8054 and CD532 favors an inactive T-loop. We show that Aurora A binds TPX2 and MLN8054 simultaneously and provide a new model for kinase conformational behavior. Our approach will enable conformation-specific effects to be integrated into inhibitor discovery across the kinome, and we outline some immediate consequences for structure-based drug discovery.

Unbinding Kinetics of a p38 MAP Kinase Type II Inhibitor from Metadynamics Simulations

Understanding the structural and energetic requisites of ligand binding toward its molecular target is of paramount relevance in drug design. In recent years, atomistic free energy calculations have proven to be a valid tool to complement experiments in characterizing the thermodynamic and kinetic properties of protein/ligand interaction. Here, we investigate, through a recently developed metadynamics-based protocol, the unbinding mechanism of an inhibitor of the pharmacologically relevant target p38 MAP kinase. We provide a thorough description of the ligand unbinding pathway identifying the most stable binding mode and other thermodynamically relevant poses. From our simulations, we estimated the unbinding rate as k_off = 0.020 ± 0.011 s^-1. This is in good agreement with the experimental value (k_off = 0.14 s^-1). Next, we developed a Markov state model that allowed identifying the rate-limiting step of the ligand unbinding process. Our calculations further show that the solvation of the ligand and that of the active site play crucial roles in the unbinding process. This study paves the way to investigations on the unbinding dynamics of more complex p38 inhibitors and other pharmacologically relevant inhibitors in general, demonstrating that metadynamics can be a powerful tool in designing new drugs with engineered binding/unbinding kinetics.

Elucidation of a four-site allosteric network in fibroblast growth factor receptor tyrosine kinases

Receptor tyrosine kinase (RTK) signaling is tightly regulated by protein allostery within the intracellular tyrosine kinase domains. Yet the molecular determinants of allosteric connectivity in tyrosine kinase domain are incompletely understood. By means of structural (X-ray and NMR) and functional characterization of pathogenic gain-of-function mutations affecting the FGF receptor (FGFR) tyrosine kinase domain, we elucidated a long-distance allosteric network composed of four interconnected sites termed the 'molecular brake', 'DFG latch', 'A-loop plug', and 'αC tether'. The first three sites repress the kinase from adopting an active conformation, whereas the αC tether promotes the active conformation. The skewed design of this four-site allosteric network imposes tight autoinhibition and accounts for the incomplete mimicry of the activated conformation by pathogenic mutations targeting a single site. Based on the structural similarity shared among RTKs, we propose that this allosteric model for FGFR kinases is applicable to other RTKs.

Chemical Proteomics and Structural Biology Define EPHA2 Inhibition by Clinical Kinase Drugs

The receptor tyrosine kinase EPHA2 (Ephrin type-A receptor 2) plays important roles in oncogenesis, metastasis, and treatment resistance, yet therapeutic targeting, drug discovery, or investigation of EPHA2 biology is hampered by the lack of appropriate inhibitors and structural information. Here, we used chemical proteomics to survey 235 clinical kinase inhibitors for their kinase selectivity and identified 24 drugs with submicromolar affinities for EPHA2. NMR-based conformational dynamics together with nine new cocrystal structures delineated drug-EPHA2 interactions in full detail. The combination of selectivity profiling, structure determination, and kinome wide sequence alignment allowed the development of a classification system in which amino acids in the drug binding site of EPHA2 are categorized into key, scaffold, potency, and selectivity residues. This scheme should be generally applicable in kinase drug discovery, and we anticipate that the provided information will greatly facilitate the development of selective EPHA2 inhibitors in particular and the repurposing of clinical kinase inhibitors in general.

Single versus dual-binding conformations in cellulosomal cohesin-dockerin complexes

Cohesins and dockerins are complementary interacting protein modules that form stable and highly specific receptor-ligand complexes. They play a crucial role in the assembly of cellulose-degrading multi-enzyme complexes called cellulosomes and have potential applicability in several technology areas, including biomass conversion processes. Here, we describe several exceptional properties of cohesin-dockerin complexes, including their tenacious biochemical affinity, remarkably high mechanostability and a dual-binding mode of recognition that is contrary to the conventional lock-and-key model of receptor-ligand interactions. We focus on structural aspects of the dual mode of cohesin-dockerin binding, highlighting recent single-molecule analysis techniques for its explicit characterization.

NMR Characterization of Information Flow and Allosteric Communities in the MAP Kinase p38γ

The intramolecular network structure of a protein provides valuable insights into allosteric sites and communication pathways. However, a straightforward method to comprehensively map and characterize these pathways is not currently available. Here we present an approach to characterize intramolecular network structure using NMR chemical shift perturbations. We apply the method to the mitogen activated protein kinase (MAPK) p38γ. p38γ contains allosteric sites that are conserved among eukaryotic kinases as well as unique to the MAPK family. How these regulatory sites communicate with catalytic residues is not well understood. Using our method, we observe and characterize for the first time information flux between regulatory sites through a conserved kinase infrastructure. This network is accessed, reinforced, and broken in various states of p38γ, reflecting the functional state of the protein. We demonstrate that the approach detects critical junctions in the network corresponding to biologically significant allosteric sites and pathways.

Conformational Selection and Induced Fit Mechanisms in the Binding of an Anticancer Drug to the c-Src Kinase

Understanding the conformational changes associated with the binding of small ligands to their biological targets is a fascinating and meaningful question in chemistry, biology and drug discovery. One of the most studied and important is the so-called "DFG-flip" of tyrosine kinases. The conserved three amino-acid DFG motif undergoes an "in to out" movement resulting in a particular inactive conformation to which "type II" kinase inhibitors, such as the anti-cancer drug Imatinib, bind. Despite many studies, the details of this prototypical conformational change are still debated. Here we combine various NMR experiments and surface plasmon resonance with enhanced sampling molecular dynamics simulations to shed light into the conformational dynamics associated with the binding of Imatinib to the proto-oncogene c-Src. We find that both conformational selection and induced fit play a role in the binding mechanism, reconciling opposing views held in the literature. Moreover, an external binding pose and local unfolding (cracking) of the aG helix are observed.

Fluorescence polarization-based assays for detecting compounds binding to inactive c-Jun N-terminal kinase 3 and p38α mitogen-activated protein kinase

Two fluorescein-labeled pyridinylimidazoles were synthesized and evaluated as probes for the binding affinity determination of potential kinase inhibitors to the c-Jun N-terminal kinase 3 (JNK3) and p38α mitogen-activated protein kinase (MAPK). Fluorescence polarization (FP)-based competition binding assays were developed for both enzymes using 1-(3',6'-dihydroxy-3-oxo-3H-spiro[isobenzofuran-1,9'-xanthen]-5-yl)-3-(4-((4-(4-(4-fluorophenyl)-2-(methylthio)-1H-imidazol-5-yl)pyridin-2-yl)amino)phenyl)thiourea (5) as an FP probe (JNK3: Kd = 3.0 nM; p38α MAPK: Kd = 5.7 nM). The validation of the assays with known inhibitors of JNK3 and p38α MAPK revealed that both FP assays correlate very well with inhibition data received by the activity assays. This, in addition to the viability of both FP-based binding assays for the high-throughput screening procedure, makes the assays suitable as inexpensive prescreening protocols for JNK3 and p38α MAPK inhibitors.

Allosteric inhibition of c-Met kinase in sub-microsecond molecular dynamics simulations induced by its inhibitor, tivantinib

Molecular dynamics simulations showed that conformation transition of c-Met from DFG-in to DFG-out may accomplish rapidly in the presence of tivantinib. A unique binding mode of tivantinib was found to be critical for this “DFG-flip”.

Biophysical investigation and conformational analysis of p38α kinase inhibitor doramapimod and its analogues

Doramapimod (BIRB 796) is a potent inhibitor of p38α nitrogen-activated protein kinase. By using biophysical methods, a clear correlation between kinase binding and the torsion angle θ of doramapimod analogues was found, highlighting the importance of inhibitor conformation for protein binding.

p38α MAPK and Type I Inhibitors: Binding Site Analysis and Use of Target Ensembles in Virtual Screening

Mitogen-activated protein kinase p38α plays an essential role in the regulation of pro-inflammatory signaling, and selective blockade of this kinase could be efficacious in many pathological processes. Despite considerable research efforts focused on the discovery and development of p38α MAPK inhibitors, no drug targeting this protein has been approved for clinical use so far. We herein analyze the available crystal structures of p38α MAPK in complex with ATP competitive type I inhibitors, getting insights into ATP binding site conformation and its influence on automated molecular docking results. The use of target ensembles, rather than single conformations, resulted in a performance improvement in both the ability to reproduce experimental bound conformations and the capability of mining active molecules from compound libraries. The information gathered from this study can be exploited in structure-based drug discovery programs having as the ultimate aim the identification of novel p38α MAPK type I inhibitors.

Identification of a Highly Conserved Allosteric Binding Site on Mnk1 and Mnk2

Elevated levels of phosphorylated eukaryotic initiation factor 4E (eIF4E) have been implicated in many tumor types, and mitogen activated protein kinase-interacting kinases (Mnks) are the only known kinases that phosphorylate eIF4E at Ser209. The phosphorylation of eIF4E is essential for oncogenic transformation but is of no significance to normal growth and development. Pharmacological inhibition of Mnks therefore provides a nontoxic and effective strategy for cancer therapy. However, a lack of specific Mnk inhibitors has confounded pharmacological target validation and clinical development. Herein, we report the identification of a novel series of Mnk inhibitors and their binding modes. A systematic workflow has been established to distinguish between type III and type I/II inhibitors. A selection of 66 compounds was tested for Mnk1 and Mnk2 inhibition, and 9 out of 20 active compounds showed type III interaction with an allosteric site of the proteins. Most of the type III inhibitors exhibited dual Mnk1 and Mnk2 activities and demonstrated potent antiproliferative properties against the MV4-11 acute myeloid leukemia cell line. Interestingly, ATP-/substrate-competitive inhibitors were found to be highly selective for Mnk2, with little or no activity for Mnk1. Our study suggests that Mnk1 and Mnk2 share a common structure of the allosteric inhibitory binding site but possess different structural features of the ATP catalytic domain. The findings will assist in the future design and development of Mnk targeted anticancer therapeutics.

Structural and dynamic insights into the energetics of activation loop rearrangement in FGFR1 kinase

Protein tyrosine kinases differ widely in their propensity to undergo rearrangements of the N-terminal Asp–Phe–Gly (DFG) motif of the activation loop, with some, including FGFR1 kinase, appearing refractory to this so-called ‘DFG flip'. Recent inhibitor-bound structures have unexpectedly revealed FGFR1 for the first time in a ‘DFG-out' state. Here we use conformationally selective inhibitors as chemical probes for interrogation of the structural and dynamic features that appear to govern the DFG flip in FGFR1. Our detailed structural and biophysical insights identify contributions from altered dynamics in distal elements, including the αH helix, towards the outstanding stability of the DFG-out complex with the inhibitor ponatinib. We conclude that the αC-β4 loop and ‘molecular brake' regions together impose a high energy barrier for this conformational rearrangement, and that this may have significance for maintaining autoinhibition in the non-phosphorylated basal state of FGFR1.

Receptor tyrosine kinases are key mediators of cell proliferation that have been implicated in several disease states for which they represent promising drug targets. Here the authors determine the thermodynamic basis for the low propensity of FGFR1 to access the DFG-Phe-out conformation required to bind type-II inhibitors.

Kinetic Cooperativity in Human Pancreatic Glucokinase Originates from Millisecond Dynamics of the Small Domain

The hallmark of glucokinase (GCK), which catalyzes the phosphorylation of glucose during glycolysis, is its kinetic cooperativity, whose understanding at atomic detail has remained open since its discovery over 40 years ago. Herein, by using kinetic CPMG NMR spectroscopic data for 17 isoleucine side chains distributed over all parts of GCK, we show that the origin of kinetic cooperativity is rooted in intramolecular protein dynamics. Residues of glucose-free GCK located in the small domain displayed distinct exchange behavior involving multiple conformers that are substantially populated (p>17 %) with a kex  value of 509±51 s(-1) , whereas in the glucose-bound form these exchange processes were quenched. This exchange behavior directly competes with the enzymatic turnover rate at physiological glucose concentrations, thereby generating the sigmoidal rate dependence that defines kinetic cooperativity.