Cyclic AMP analog blocks kinase activation by stabilizing inactive conformation: conformational selection highlights a new concept in allosteric inhibitor design

Multiplicity of Regulatory Subunit Conformations Defines Structural Ensemble of Reset Protein Kinase A Holoenzyme

How protein kinase A (PKA) is reset to a basal state following 3'5'-cyclic adenosine monophosphate (cAMP)-mediated activation is unknown. Here we describe the mechanism of cAMP-PKA type I signal termination leading to a reset of PKA by holoenzyme formation through the obligatory action of phosphodiesterases (PDEs). We report a catalytic subunit (Cα)-assisted mechanism for the reset of type I PKA and describe for the first time multiple structures of the reset PKA holoenzyme (RIα2:Cα2) that capture an ensemble of multiple conformational end-states through integrative electron microscopy and structural mass spectrometry approaches. Together these complementary methods highlight the large conformational dynamics of the regulatory subunit (RIα) within the tetrameric reset PKA holoenzyme. The cAMP-free reset PKA holoenzyme adopts multiple distinct conformations of RIα with contributions from the N-terminal linker and CNB-B dynamics. Our findings highlight the interplay between RIα, Cα, and PDEs (PDE8) in cAMP-PKA signalosomes to offer a new paradigm for PDE-mediated regulation of cAMP-PKA signaling.

Exploring the Conformational Landscape and Stability of Aurora A Using Ion-Mobility Mass Spectrometry and Molecular Modeling

Protein kinase inhibitors are highly effective in treating diseases driven by aberrant kinase signaling and as chemical tools to help dissect the cellular roles of kinase signaling complexes. Evaluating the effects of binding of small molecule inhibitors on kinase conformational dynamics can assist in understanding both inhibition and resistance mechanisms. Using gas-phase ion-mobility mass spectrometry (IM-MS), we characterize changes in the conformational landscape and stability of the protein kinase Aurora A (Aur A) driven by binding of the physiological activator TPX2 or small molecule inhibition. Aided by molecular modeling, we establish three major conformations, the relative abundances of which were dependent on the Aur A activation status: one highly populated compact conformer similar to that observed in most crystal structures, a second highly populated conformer possessing a more open structure infrequently found in crystal structures, and an additional low-abundance conformer not currently represented in the protein databank. Notably, inhibitor binding induces more compact configurations of Aur A, as adopted by the unbound enzyme, with both IM-MS and modeling revealing inhibitor-mediated stabilization of active Aur A.

Depolarization induces nociceptor sensitization by CaV1.2-mediated PKA-II activation

Depolarization drives neuronal plasticity. However, whether depolarization drives sensitization of peripheral nociceptive neurons remains elusive. By high-content screening (HCS) microscopy, we revealed that depolarization of cultured sensory neurons rapidly activates protein kinase A type II (PKA-II) in nociceptors by calcium influx through Ca_V1.2 channels. This effect was modulated by calpains but insensitive to inhibitors of cAMP formation, including opioids. In turn, PKA-II phosphorylated Ser1928 in the distal C terminus of Ca_V1.2, thereby increasing channel gating, whereas dephosphorylation of Ser1928 involved the phosphatase calcineurin. Patch-clamp and behavioral experiments confirmed that depolarization leads to calcium- and PKA-dependent sensitization of calcium currents ex vivo and local peripheral hyperalgesia in the skin in vivo. Our data suggest a local activity-driven feed-forward mechanism that selectively translates strong depolarization into further activity and thereby facilitates hypersensitivity of nociceptor terminals by a mechanism inaccessible to opioids.

Integrated Phosphoproteomics for Identifying Substrates of Human Protein Kinase A (PRKACA) and Its Oncogenic Mutant DNAJB1-PRKACA

The DNAJB1-PRKACA fusion is the signature genetic event of fibrolamellar hepatocellular carcinoma (FL-HCC), a rare but lethal liver cancer that primarily affects adolescents and young adults. A deletion fuses the first exon of the HSP40 gene (DNAJB1), with exons 2-10 of protein kinase A (PRKACA), producing the chimeric kinase DNAJB1-PKA_ca (J-PKA_ca). The HSP40 portion's scaffolding/chaperone function has been implicated in redirecting substrate recognition to upregulate oncogenic pathways, but the direct substrates of this fusion are not fully known. We integrated cell-based and in vitro phosphoproteomics to identify substrates targeted directly by PKA and J-PKA_ca, comparing phosphoproteome profiles from cells with in vitro rephosphorylation of peptides and proteins from lysates using recombinant enzymes. We identified a subset of phosphorylation sites in both cell-based and in vitro experiments, as well as altered pathways and proteins consistent with observations from related studies. We also treated cells with PKA inhibitors that function by two different mechanisms (rpcAMPs and PKI) and examined phosphoproteome profiles, finding some substrates that persisted in the presence of inhibitors and revealing differences between WT and chimera. Overall, these results provide potential insights into J-PKA_ca's oncogenic activity in a complex cellular system and may provide candidate targets for therapeutic follow-up.

A Review on Computational Approaches for Analyzing Hydrogen- Deuterium (H/D) Exchange of Proteins

Native state Hydrogen-Deuterium (H/D) exchange method has been used to study the structures and the unfolding pathways for quite a number of proteins. The H/D exchange method is generally monitored using nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) techniques. NMR-assisted H/D exchange methods primarily monitor the residue level fluctuation of proteins, whereas MS-assisted H/D exchange methods analyze multifold ensemble conformations of proteins. In this connection, quite a large number of computational tools and algorithms have been developed for processing and analyzing huge amount of the H/D exchange data generated from these techniques. In this review, most of the freely available computational tools associated with the H/D exchange of proteins have been comprehensively reviewed and scopes to improve/ develop novel computational approaches for analyzing the H/D exchange data of proteins have also been brought into fore.

Noncanonical protein kinase A activation by oligomerization of regulatory subunits as revealed by inherited Carney complex mutations

Familial mutations of the protein kinase A (PKA) R1α regulatory subunit lead to a generalized predisposition for a wide range of tumors, from pituitary adenomas to pancreatic and liver cancers, commonly referred to as Carney complex (CNC). CNC mutations are known to cause overactivation of PKA, but the molecular mechanisms underlying such kinase overactivity are not fully understood in the context of the canonical cAMP-dependent activation of PKA. Here, we show that oligomerization-induced sequestration of R1α from the catalytic subunit of PKA (C) is a viable mechanism of PKA activation that can explain the CNC phenotype. Our investigations focus on comparative analyses at the level of structure, unfolding, aggregation, and kinase inhibition profiles of wild-type (wt) PKA R1α, the A211D and G287W CNC mutants, as well as the cognate acrodysostosis type 1 (ACRDYS1) mutations A211T and G287E. The latter exhibit a phenotype opposite to CNC with suboptimal PKA activation compared with wt. Overall, our results show that CNC mutations not only perturb the classical cAMP-dependent allosteric activation pathway of PKA, but also amplify significantly more than the cognate ACRDYS1 mutations nonclassical and previously unappreciated activation pathways, such as oligomerization-induced losses of the PKA R1α inhibitory function.

Adenylate control in cAMP signaling: implications for adaptation in signalosomes

In cAMP-Protein Kinase A (PKA) signaling, A-kinase anchoring protein scaffolds assemble PKA in close proximity to phosphodiesterases (PDE), kinase-substrates to form signaling islands or 'signalosomes'. In its basal state, inactive PKA holoenzyme (R2:C2) is activated by binding of cAMP to regulatory (R)-subunits leading to dissociation of active catalytic (C)-subunits. PDEs hydrolyze cAMP-bound to the R-subunits to generate 5'-AMP for termination and resetting the cAMP signaling. Mechanistic basis for cAMP signaling has been derived primarily by focusing on the proteins in isolation. Here, we set out to simulate cAMP signaling activation-termination cycles in a signalosome-like environment with PDEs and PKA subunits in close proximity to each other. Using a combination of fluorescence polarization and amide hydrogen exchange mass spectrometry with regulatory (RIα), C-subunit (Cα) and PDE8 catalytic domain, we have tracked movement of cAMP through activation-termination cycles. cAMP signaling operates as a continuum of four phases: (1) Activation and dissociation of PKA into R- and C-subunits by cAMP and facilitated by substrate (2) PDE recruitment to R-subunits (3) Hydrolysis of cAMP to 5'-AMP (4) Reassociation of C-subunit to 5'-AMP-bound-RIα in the presence of excess ATP to reset cAMP signaling to form the inactive PKA holoenzyme. Our results demonstrate that 5'-AMP is not merely a passive hydrolysis end-product of PDE action. A 'ligand-free' state R subunit does not exist in signalosomes as previously assumed. Instead the R-subunit toggles between cAMP- or 5'-AMP bound forms. This highlights, for the first time, the importance of 5'-AMP in promoting adaptation and uncovers adenylate control in cAMP signaling.

Allosteric inhibition explained through conformational ensembles sampling distinct "mixed" states

Allosteric modulation provides an effective avenue for selective and potent enzyme inhibition. Here, we summarize and critically discuss recent advances on the mechanisms of allosteric partial agonists for three representative signalling enzymes activated by cyclic nucleotides: the cAMP-dependent protein kinase (PKA), the cGMP-dependent protein kinase (PKG), and the exchange protein activated by cAMP (EPAC). The comparative analysis of partial agonism in PKA, PKG and EPAC reveals a common emerging theme, i.e. the sampling of distinct “mixed” conformational states, either within a single domain or between distinct domains. Here, we show how such “mixed” states play a crucial role in explaining the observed functional response, i.e. partial agonism and allosteric pluripotency, as well as in maximizing inhibition while minimizing potency losses. In addition, by combining Nuclear Magnetic Resonance (NMR), Molecular Dynamics (MD) simulations and Ensemble Allosteric Modeling (EAM), we also show how to map the free-energy landscape of conformational ensembles containing “mixed” states. By discussing selected case studies, we illustrate how MD simulations and EAM complement NMR to quantitatively relate protein dynamics to function. The resulting NMR- and MD-based EAMs are anticipated to inform not only the design of new generations of highly selective allosteric inhibitors, but also the choice of multidrug combinations.

Synergistic Allostery in Multiligand-Protein Interactions

Amide hydrogen-deuterium exchange mass spectrometry is powerful for describing combinatorial coupling effects of a cooperative ligand pair binding at noncontiguous sites: adenosine at the ATP-pocket and a docking peptide (PIFtide) at the PIF-pocket, on a model protein kinase PDK1. Binding of two ligands to PDK1 reveal multiple hotspots of synergistic allostery with cumulative effects greater than the sum of individual effects mediated by each ligand. We quantified this synergism and ranked these hotspots using a difference in deuteration-based approach, which showed that the strongest synergistic effects were observed at three of the critical catalytic loci of kinases: the αB-αC helices, and HRD-motif loop, and DFG-motif. Additionally, we observed weaker synergistic effects at a distal GHI-subdomain locus. Synergistic changes in deuterium exchange observed at a distal site but not at the intermediate sites of the large lobe of the kinase reveals allosteric propagation in proteins to operate through two modes. Direct electrostatic interactions between polar and charged amino acids that mediate targeted relay of allosteric signals, and diffused relay of allosteric signals through soft matter-like hydrophobic core amino acids. Furthermore, we provide evidence that the conserved β-3 strand lysine of protein kinases (Lys111 of PDK1) functions as an integrator node to coordinate allosteric coupling of the two ligand-binding sites. It maintains indirect interactions with the ATP-pocket and mediates a critical salt bridge with a glutamate (Glu130) of αC helix, which is conserved across all kinases. In summary, allosteric propagation in cooperative, dual-liganded enzyme targets is bidirectional and synergistic and offers a strategy for combinatorial drug development.

A substitution in cGMP-dependent protein kinase 1 associated with aortic disease induces an active conformation in the absence of cGMP

Type 1 cGMP-dependent protein kinases (PKGs) play important roles in human cardiovascular physiology, regulating vascular tone and smooth-muscle cell phenotype. A mutation in the human PRKG1 gene encoding cGMP-dependent protein kinase 1 (PKG1) leads to thoracic aortic aneurysms and dissections. The mutation causes an arginine-to-glutamine (RQ) substitution within the first cGMP-binding pocket in PKG1. This substitution disrupts cGMP binding to the pocket, but it also unexpectedly causes PKG1 to have high activity in the absence of cGMP via an unknown mechanism. Here, we identified the molecular mechanism whereby the RQ mutation increases basal kinase activity in the human PKG1α and PKG1β isoforms. Although we found that the RQ substitution (R177Q in PKG1α and R192Q in PKG1β) increases PKG1α and PKG1β autophosphorylation in vitro, we did not detect increased autophosphorylation of the PKG1α or PKG1β RQ variant isolated from transiently transfected 293T cells, indicating that increased basal activity of the RQ variants in cells was not driven by PKG1 autophosphorylation. Replacement of Arg-177 in PKG1α with alanine or methionine also increased basal activity. PKG1 exists as a parallel homodimer linked by an N-terminal leucine zipper, and we show that the WT chain in WT-RQ heterodimers partly reduces basal activity of the RQ chain. Using hydrogen/deuterium-exchange MS, we found that the RQ substitution causes PKG1β to adopt an active conformation in the absence of cGMP, similar to that of cGMP-bound WT enzyme. We conclude that the RQ substitution in PKG1 increases its basal activity by disrupting the formation of an inactive conformation.

Allosteric pluripotency as revealed by protein kinase A

The functional response of a signaling system to an allosteric stimulus often depends on subcellular conditions, a phenomenon known as pluripotent allostery. For example, a single allosteric modulator, Rp-cAMPS, of the prototypical protein kinase A (PKA) switches from antagonist to agonist depending on MgATP levels. However, the mechanism underlying such pluripotent allostery has remained elusive for decades. Using nuclear magnetic resonance spectroscopy, ensemble models, kinase assays, and molecular dynamics simulations, we show that allosteric pluripotency arises from surprisingly divergent responses of highly homologous tandem domains. The differential responses perturb domain-domain interactions and remodel the free-energy landscape of inhibitory excited states sampled by the regulatory subunit of PKA. The resulting activation threshold values are comparable to the effective free energy of regulatory and catalytic subunit binding, which depends on metabolites, substrates, and mutations, explaining pluripotent allostery and warranting a general redefinition of allosteric targets to include specific subcellular environments.

Recent Advances in EPAC-Targeted Therapies: A Biophysical Perspective

The universal second messenger cAMP regulates diverse intracellular processes by interacting with ubiquitously expressed proteins, such as Protein Kinase A (PKA) and the Exchange Protein directly Activated by cAMP (EPAC). EPAC is implicated in multiple pathologies, thus several EPAC-specific inhibitors have been identified in recent years. However, the mechanisms and molecular interactions underlying the EPAC inhibition elicited by such compounds are still poorly understood. Additionally, being hydrophobic low molecular weight species, EPAC-specific inhibitors are prone to forming colloidal aggregates, which result in non-specific aggregation-based inhibition (ABI) in aqueous systems. Here, we review from a biophysical perspective the molecular basis of the specific and non-specific interactions of two EPAC antagonists-CE3F4R, a non-competitive inhibitor, and ESI-09, a competitive inhibitor of EPAC. Additionally, we discuss the value of common ABI attenuators (e.g., TX and HSA) to reduce false positives at the expense of introducing false negatives when screening aggregation-prone compounds. We hope this review provides the EPAC community effective criteria to evaluate similar compounds, aiding in the optimization of existing drug leads, and informing the development of the next generation of EPAC-specific inhibitors.

Advances in MS Based Strategies for Probing Ligand-Target Interactions: Focus on Soft Ionization Mass Spectrometric Techniques

The non-covalent interactions between small drug molecules and disease-related proteins (ligand-target interactions) mediate various pharmacological processes in the treatment of different diseases. The development of the analytical methods to assess those interactions, including binding sites, binding energies, stoichiometry and association-dissociation constants, could assist in clarifying the mechanisms of action, precise treatment of targeted diseases as well as the targeted drug discovery. For the last decades, mass spectrometry (MS) has been recognized as a powerful tool to study the non-covalent interactions of the ligand-target complexes with the characteristics of high sensitivity, high-resolution, and high-throughput. Soft ionization mass spectrometry, especially the electrospray mass spectrometry (ESI-MS) and matrix assisted laser desorption ionization mass spectrometry (MALDI-MS), could achieve the complete transformation of the target analytes into the gas phase, and subsequent detection of the small drug molecules and disease-related protein complexes, and has exerted great advantages for studying the drug ligands-protein targets interactions, even in case of identifying active components as drug ligands from crude extracts of medicinal plants. Despite of other analytical techniques for this purpose, such as the NMR and X-ray crystallography, this review highlights the principles, research hotspots and recent applications of the soft ionization mass spectrometry and its hyphenated techniques, including hydrogen-deuterium exchange mass spectrometry (HDX-MS), chemical cross-linking mass spectrometry (CX-MS), and ion mobility spectrometry mass spectrometry (IMS-MS), in the study of the non-covalent interactions between small drug molecules and disease-related proteins.

Structures of the PKA RIα Holoenzyme with the FLHCC Driver J-PKAcα or Wild-Type PKAcα

Fibrolamellar hepatocellular carcinoma (FLHCC) is driven by J-PKAcα, a kinase fusion chimera of the J domain of DnaJB1 with PKAcα, the catalytic subunit of protein kinase A (PKA). Here we report the crystal structures of the chimeric fusion RIα2:J-PKAcα2 holoenzyme formed by J-PKAcα and the PKA regulatory (R) subunit RIα, and the wild-type (WT) RIα2:PKAcα2 holoenzyme. The chimeric and WT RIα holoenzymes have quaternary structures different from the previously solved WT RIβ and RIIβ holoenzymes. The WT RIα holoenzyme showed the same configuration as the chimeric RIα2:J-PKAcα2 holoenzyme and a distinct second conformation. The J domains are positioned away from the symmetrical interface between the two RIα:J-PKAcα heterodimers in the chimeric fusion holoenzyme and are highly dynamic. The structural and dynamic features of these holoenzymes enhance our understanding of the fusion chimera protein J-PKAcα that drives FLHCC as well as the isoform specificity of PKA.

Mechanism of Selective Enzyme Inhibition through Uncompetitive Regulation of an Allosteric Agonist

Classical uncompetitive inhibitors are potent pharmacological modulators of enzyme function. Since they selectively target enzyme-substrate complexes (E:S), their inhibitory potency is amplified by increasing substrate concentrations. Recently, an unconventional uncompetitive inhibitor, called CE3F4R, was discovered for the exchange protein activated by cAMP isoform 1 (EPAC1). Unlike conventional uncompetitive inhibitors, CE3F4R is uncompetitive with respect to an allosteric effector, cAMP, as opposed to the substrate (i.e., CE3F4R targets the E:cAMP rather than the E:S complex). However, the mechanism of CE3F4R as an uncompetitive inhibitor is currently unknown. Here, we elucidate the mechanism of CE3F4R's action using NMR spectroscopy. Due to limited solubility and line broadening, which pose major challenges for traditional structural determination approaches, we resorted to a combination of protein- and ligand-based NMR experiments to comparatively analyze EPAC mutations, inhibitor analogs, and cyclic nucleotide derivatives that trap EPAC at different stages of activation. We discovered that CE3F4R binds within the EPAC cAMP-binding domain (CBD) at a subdomain interface distinct from the cAMP binding site, acting as a wedge that stabilizes a cAMP-bound mixed-intermediate. The mixed-intermediate includes attributes of both the apo/inactive and cAMP-bound/active states. In particular, the intermediate targeted by CE3F4R traps a CBD's hinge helix in its inactive conformation, locking EPAC into a closed domain topology that restricts substrate access to the catalytic domain. The proposed mechanism of action also explains the isoform selectivity of CE3F4R in terms of a single EPAC1 versus EPAC2 amino acid difference that destabilizes the active conformation of the hinge helix.

Reversing allosteric communication: From detecting allosteric sites to inducing and tuning targeted allosteric response

The omnipresence of allosteric regulation together with the fundamental role of structural dynamics in this phenomenon have initiated a great interest to the detection of regulatory exosites and design of corresponding effectors. However, despite a general consensus on the key role of dynamics most of the earlier efforts on the prediction of allosteric sites are heavily crippled by the static nature of the underlying methods, which are either structure-based approaches seeking for deep surface pockets typical for “traditional” orthosteric drugs or sequence-based techniques exploiting the conservation of protein sequences. Because of the critical role of global protein dynamics in allosteric signaling, we investigate the hypothesis of reversibility in allosteric communication, according to which allosteric sites can be detected via the perturbation of the functional sites. The reversibility is tested here using our structure-based perturbation model of allostery, which allows one to analyze the causality and energetics of allosteric communication. We validate the “reverse perturbation” hypothesis and its predictive power on a set of classical allosteric proteins, then, on the independent extended benchmark set. We also show that, in addition to known allosteric sites, the perturbation of the functional sites unravels rather extended protein regions, which can host latent regulatory exosites. These protein parts that are dynamically coupled with functional sites can also be used for inducing and tuning allosteric communication, and an exhaustive exploration of the per-residue contributions to allosteric effects can eventually lead to the optimal modulation of protein activity. The site-effector interactions necessary for a specific mode and level of allosteric communication can be fine-tuned by adjusting the site’s structure to an available effector molecule and by the design or selection of an appropriate ligand.

Recent advances in the development of allosteric drugs allow one to fully appreciate the sheer power of allosteric effectors in the avoiding toxicity, receptor desensitization and modulatory rather than on/off mode of action, compared to the traditional orthosteric compounds. The detection of allosteric sites is one of the major challenges in the quest for allosteric drugs. This work proposes a “reverse perturbation” approach for identifying allosteric sites as a result of a perturbation applied to the functional ones. We show that according to the traditional Monod-Changeux-Jacob’s definition of allostery, considering non-overlapping regulatory and functional sites is a critical prerequisite for the successful detection of allosteric sites. Using the reverse perturbation method, it is possible to determine wide protein regions with a potential to induce an allosteric response and to adjust its strength. Further studies on inducing and fine-tuning of allosteric signalling seem to be of a great importance for efficient design of non-orthosteric ligands in the development of novel drugs.

PKA-RII subunit phosphorylation precedes activation by cAMP and regulates activity termination

Activity of endogenous protein kinase A (PKA) could never be analyzed directly in the cellular environment. Isensee et al. used antibodies to quantify conformational changes leading to an open conformation of endogenous PKA-II holoenzymes, which allowed them to analyze and model its activation cycle in primary sensory neurons.

Type II isoforms of cyclic adenosine monophosphate (cAMP)–dependent protein kinase A (PKA-II) contain a phosphorylatable epitope within the inhibitory domain of RII subunits (pRII) with still unclear function. In vitro, RII phosphorylation occurs in the absence of cAMP, whereas staining of cells with pRII-specific antibodies revealed a cAMP-dependent pattern. In sensory neurons, we found that increased pRII immunoreactivity reflects increased accessibility of the already phosphorylated RII epitope during cAMP-induced opening of the tetrameric RII₂:C₂ holoenzyme. Accordingly, induction of pRII by cAMP was sensitive to novel inhibitors of dissociation, whereas blocking catalytic activity was ineffective. Also in vitro, cAMP increased the binding of pRII antibodies to RII₂:C₂ holoenzymes. Identification of an antibody specific for the glycine-rich loop of catalytic subunits facing the pRII-epitope confirmed activity-dependent binding with similar kinetics, proving that the reassociation is rapid and precisely controlled. Mechanistic modeling further supported that RII phosphorylation precedes cAMP binding and controls the inactivation by modulating the reassociation involving the coordinated action of phosphodiesterases and phosphatases.

Neutron Crystallography Detects Differences in Protein Dynamics: Structure of the PKG II Cyclic Nucleotide Binding Domain in Complex with an Activator

As one of the main receptors of a second messenger, cGMP, cGMP-dependent protein kinase (PKG) isoforms I and II regulate distinct physiological processes. The design of isoform-specific activators is thus of great biomedical importance and requires detailed structural information about PKG isoforms bound with activators, including accurate positions of hydrogen atoms and a description of the hydrogen bonding and water architecture. Here, we determined a 2.2 Å room-temperature joint X-ray/neutron (XN) structure of the human PKG II carboxyl cyclic nucleotide binding (CNB-B) domain bound with a potent PKG II activator, 8-pCPT-cGMP. The XN structure directly visualizes intermolecular interactions and reveals changes in hydrogen bonding patterns upon comparison to the X-ray structure determined at cryo-temperatures. Comparative analysis of the backbone hydrogen/deuterium exchange patterns in PKG II:8-pCPT-cGMP and previously reported PKG Iβ:cGMP XN structures suggests that the ability of these agonists to activate PKG is related to how effectively they quench dynamics of the cyclic nucleotide binding pocket and the surrounding regions.

A novel allosteric site in casein kinase 2α discovered using combining bioinformatics and biochemistry methods

Casein kinase 2 (CK2) is a highly pleiotropic serine-threonine kinase, which catalyzed phosphorylation of more than 300 proteins that are implicated in regulation of many cellular functions, such as signal transduction, transcriptional control, apoptosis and the cell cycle. On the other hand, CK2 is abnormally elevated in a variety of tumors, and is considered as a promising therapeutic target. The currently available ATP-competitive CK2 inhibitors, however, lack selectivity, which has impeded their development in cancer therapy. Because allosteric inhibitors can avoid the shortcomings of conventional kinase inhibitors, this study was aimed to discover a new allosteric site in CK2α and to investigate the effects of mutations in this site on the activity of CK2α. Using Allosite based on protein dynamics and structural alignment, we predicted a new allosteric site that was partly located in the αC helix of CK2α. Five residues exposed on the surface of this site were mutated to validate the prediction. Kinetic analyses were performed using a luminescent ADP detection assay by varying the concentrations of a peptide substrate, and the results showed that the mutations I78C and I78W decreased CK2α activity, whereas V31R, K75E, I82C and P109C increased CK2α activity. Potential allosteric pathways were identified using the Monte Carlo path generation approach, and the results of these predicted allosteric pathways were consistent with the mutation analysis. Multiple sequence alignments of CK2α with the other kinases in the family were conducted using the ClustalX method, which revealed the diversity of the residues in the site. In conclusion, we identified a new allosteric site in CK2α that can be altered to modulate the activity of the kinase. Because of the high diversity of the residues in the site, the site can be targeted using rational drug design of specific CK2α inhibitors for biological relevance.

Molecular Simulations Reveal an Unresolved Conformation of the Type IA Protein Kinase A Regulatory Subunit and Suggest Its Role in the cAMP Regulatory Mechanism

We identify a previously unresolved, unrecognized, and highly stable conformation of the protein kinase A (PKA) regulatory subunit RIα. This conformation, which we term the "Flipback" structure, bridges conflicting characteristics in crystallographic structures and solution experiments of the PKA RIα heterotetramer. Our simulations reveal a hinge residue, G235, in the B/C helix that is conserved through all isoforms of RI. Brownian dynamics simulations suggest that the Flipback conformation plays a role in cAMP association to the A domain of the R subunit.

Dissecting Orthosteric Contacts for a Reverse-Fragment-Based Ligand Design

Orthosteric sites on proteins are formed typically from noncontiguous interacting sites in three-dimensional space where the composite binding interaction of a biological ligand is mediated by multiple synergistic interactions of its constituent functional groups. Through these multiple interactions, ligands stabilize both the ligand binding site and the local secondary structure. However, relative energetic contributions of the individual contacts in these protein-ligand interactions are difficult to resolve. Deconvolution of the contributions of these various functional groups in natural inhibitors/ligand would greatly aid in iterative fragment-based drug discovery (FBDD). In this study, we describe an approach of progressive unfolding of a target protein using a gradient of denaturant urea to reveal the individual energetic contributions of various ligand-functional groups to the affinity of the entire ligand. Through calibrated unfolding of two protein-ligand systems: cAMP-bound regulatory subunit of Protein Kinase A (RIα) and IBMX-bound phosphodiesterase8 (PDE8), monitored by amide hydrogen-deuterium exchange mass spectrometry, we show progressive disruption of individual orthosteric contacts in the ligand binding sites, allowing us to rank the energetic contributions of these individual interactions. In the two cAMP-binding sites of RIα, exocyclic phosphate oxygens of cAMP were identified to mediate stronger interactions than ribose 2'-OH in both the RIα-cAMP binding interfaces. Further, we have also ranked the relative contributions of the different functional groups of IBMX based on their interactions with the orthosteric residues of PDE8. This strategy for deconstruction of individual binding sites and identification of the strongest functional group interaction in enzyme orthosteric sites offers a rational starting point for FBDD.

Electrostatic Switch Function in the Mechanism of Protein Kinase A I Activation: Results of the Molecular Dynamics Simulation

We used molecular dynamics to find the average path of the A-domain H → B conformational transition in protein kinase A Iα. We obtained thirteen productive trajectories and processed them sequentially using factor and cross-correlation analyses. The conformational transition is presented as partly deterministic sequence of six events. Event B represents H → B transition of the phosphate binding cassette. Main participants of this event form electrostatic switch cAMP(O6)–A202(N-H)–G199(C=O). Through this switch, cAMP transmits information about its binding to hydrophobic switch L203–Y229 and thus triggers conformational transition of A-domain. Events C and D consist in N3A-motif displacement towards phosphate binding cassette and B/C-helix rotation. Event E involves an increase in interaction energy between Y229 and β-subdomain. Taken together, events B, E, and D correspond to the hinge movement towards β-barrel. Transition of B/C-helix turn (a.a. 229–234) from α-form to π-form accounts for event F. Event G implies that π-helical turn is replaced by kink. Emerging in the resulting conformation, electrostatic interaction R241–E200 facilitates kink formation. The obtained data on the mechanism of cAMP-dependent activation of PKA Iα may contribute to new approaches to designing pharmaceuticals based on cAMP analogs.

Crystal structure of cGMP-dependent protein kinase Iβ cyclic nucleotide-binding-B domain : Rp-cGMPS complex reveals an apo-like, inactive conformation

The R-diastereomer of phosphorothioate analogs of cGMP, Rp-cGMPS, is one of few known inhibitors of cGMP-dependent protein kinase I (PKG I); however, its mechanism of inhibition is currently not fully understood. Here, we determined the crystal structure of the PKG Iβ cyclic nucleotide-binding domain (PKG Iβ CNB-B), considered a 'gatekeeper' for cGMP activation, bound to Rp-cGMPS at 1.3 Å. Our structural and NMR data show that PKG Iβ CNB-B bound to Rp-cGMPS displays an apo-like structure with its helical domain in an open conformation. Comparison with the cAMP-dependent protein kinase regulatory subunit (PKA RIα) showed that this conformation resembles the catalytic subunit-bound inhibited state of PKA RIα more closely than the apo or Rp-cAMPS-bound conformations. These results suggest that Rp-cGMPS inhibits PKG I by stabilizing the inactive conformation of CNB-B.

β-subunit myristoylation functions as an energy sensor by modulating the dynamics of AMP-activated Protein Kinase

The heterotrimeric AMP-activated protein kinase (AMPK), consisting of α, β and γ subunits, is a stress-sensing enzyme that is activated by phosphorylation of its activation loop in response to increases in cellular AMP. N-terminal myristoylation of the β-subunit has been shown to suppress Thr172 phosphorylation, keeping AMPK in an inactive state. Here we use amide hydrogen-deuterium exchange mass spectrometry (HDX-MS) to investigate the structural and dynamic properties of the mammalian myristoylated and non-myristoylated inactivated AMPK (D139A) in the presence and absence of nucleotides. HDX MS data suggests that the myristoyl group binds near the first helix of the C-terminal lobe of the kinase domain similar to other kinases. Our data, however, also shows that ATP.Mg²⁺ results in a global stabilization of myristoylated, but not non-myristoylated AMPK, and most notably for peptides of the activation loop of the α-kinase domain, the autoinhibitory sequence (AIS) and the βCBM. AMP does not have that effect and HDX measurements for myristoylated and non-myristoylated AMPK in the presence of AMP are similar. These differences in dynamics may account for a reduced basal rate of phosphorylation of Thr172 in myristoylated AMPK in skeletal muscle where endogenous ATP concentrations are very high.

Ligand-mediated changes in conformational dynamics of NpmA: implications for ribosomal interactions

Aminoglycosides are broad-spectrum antibiotics that bind to the 30S ribosomal subunit (30S) of bacteria and disrupt protein translation. NpmA, a structurally well-characterized methyltransferase identified in an E. coli clinical isolate, catalyzes methylation of 30S at A1408 of the 16S rRNA and confers aminoglycoside resistance. Using sucrose cushion centrifugation and isothermal titration calorimetry, we first confirmed the binding between NpmA and 30S. Next, we performed amide Hydrogen/Deuterium Exchange Mass Spectrometry (HDXMS) of apo NpmA and in the presence and absence of SAM/SAH. We observed that ligand binding resulted in time-dependent differences in deuterium exchange not only at the ligand-binding pocket (D25–D55 and A86–E112) but also in distal regions (F62-F82 and Y113-S144) of NpmA. These results provide insights into methylation group donor cofactor-mediated allostery in NpmA in the ligand-bound states, which could not be observed in the static endpoint crystal structures. We predict that the two distal sites in NpmA form part of the allosteric sites that importantly are part of the main 16S rRNA binding interface. Thus HDXMS helped uncover allosteric communication relays that couple SAM/SAH binding sites with the ribosome-binding site. This highlights how HDXMS together with X-ray crystallography can provide important allosteric insights in protein-ligand complexes.

Disrupting the Allosteric Interaction between the Plasmodium falciparum cAMP-dependent Kinase and Its Regulatory Subunit

The ubiquitous second messenger cAMP mediates signal transduction processes in the malarial parasite that regulate host erythrocyte invasion and the proliferation of merozoites. In Plasmodium falciparum, the central receptor for cAMP is the single regulatory subunit (R) of protein kinase A (PKA). To aid the development of compounds that can selectively dysregulate parasite PKA signaling, we solved the structure of the PKA regulatory subunit in complex with cAMP and a related analogue that displays antimalarial activity, (S_p)-2-Cl-cAMPS. Prior to signaling, PKA-R holds the kinase's catalytic subunit (C) in an inactive state by exerting an allosteric inhibitory effect. When two cAMP molecules bind to PKA-R, they stabilize a structural conformation that facilitates its dissociation, freeing PKA-C to phosphorylate downstream substrates such as apical membrane antigen 1. Although PKA activity was known to be necessary for erythrocytic proliferation, we show that uncontrolled induction of PKA activity using membrane-permeable agonists is equally disruptive to growth.

Predicting Allosteric Effects from Orthosteric Binding in Hsp90-Ligand Interactions: Implications for Fragment-Based Drug Design

A key question in mapping dynamics of protein-ligand interactions is to distinguish changes at binding sites from those associated with long range conformational changes upon binding at distal sites. This assumes a greater challenge when considering the interactions of low affinity ligands (dissociation constants, K_D, in the μM range or lower). Amide hydrogen deuterium Exchange mass spectrometry (HDXMS) is a robust method that can provide both structural insights and dynamics information on both high affinity and transient protein-ligand interactions. In this study, an application of HDXMS for probing the dynamics of low affinity ligands to proteins is described using the N-terminal ATPase domain of Hsp90. Comparison of Hsp90 dynamics between high affinity natural inhibitors (K_D ~ nM) and fragment compounds reveal that HDXMS is highly sensitive in mapping the interactions of both high and low affinity ligands. HDXMS reports on changes that reflect both orthosteric effects and allosteric changes accompanying binding. Orthosteric sites can be identified by overlaying HDXMS onto structural information of protein-ligand complexes. Regions distal to orthosteric sites indicate long range conformational changes with implications for allostery. HDXMS, thus finds powerful utility as a high throughput method for compound library screening to identify binding sites and describe allostery with important implications for fragment-based ligand discovery (FBLD).

Ligand interactions with proteins result in broad changes that are propagated throughout the target proteins, across space and time. These changes can be broadly classified into: orthosteric effects at the ligand binding site and allosteric changes at distal sites. These allosteric changes are difficult to localize and distinguish from binding interactions. In this study, we describe the application of amide hydrogen/deuterium exchange mass-spectrometry (HDXMS) to differentiate between changes occurring at the binding site and at distal allosteric sites by combining HDXMS with X-ray crystallography. Every ligand or a fragment mediates distinct contacts and results in changes in deuterium uptake across the protein. By comparing with orthosteric structural information, it is possible to identify long-range changes (action at a distance) due to the ligands. An important application of HDXMS is that it can identify subtle changes in protein dynamics that cannot be picked up by quantitative screens of protein-ligand interactions or crystal structures. This gives us the ability to describe ligand binding based on the response from different regions in the proteins. Thus it provides us with the potential to accurately measure and compare changes in dynamics upon binding different ligands and fragments, which is greatly valuable in fragment-based ligand design.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.

Tandem Domains with Tuned Interactions Are a Powerful Biological Design Principle

Allosteric effects of mutations, ligand binding, or post-translational modifications on protein function occur through changes to the protein's shape, or conformation. In a cell, there are many copies of the same protein, all experiencing these perturbations in a dynamic fashion and fluctuating through different conformations and activity states. According to the "conformational selection and population shift" theory, ligand binding selects a particular conformation. This perturbs the ensemble and induces a population shift. In a new PLOS Biology paper, Melacini and colleagues describe a novel model of protein regulation, the "Double-Conformational Selection Model", which demonstrates how two tandem ligand-binding domains interact to regulate protein function. Here we explain how tandem domains with tuned interactions-but not single domains-can provide a blueprint for sensitive activation sensors within a narrow window of ligand concentration, thereby promoting signaling control.

Limited Proteolysis Combined with Stable Isotope Labeling Reveals Conformational Changes in Protein (Pseudo)kinases upon Binding Small Molecules

Likely due to conformational rearrangements, small molecule inhibitors may stabilize the active conformation of protein kinases and paradoxically promote tumorigenesis. We combined limited proteolysis with stable isotope labeling MS to monitor protein conformational changes upon binding of small molecules. Applying this method to the human serine/threonine kinase B-Raf, frequently mutated in cancer, we found that binding of ATP or its nonhydrolyzable analogue AMP-PNP, but not ADP, stabilized the structure of both B-Raf(WT) and B-Raf(V600E). The ATP-competitive type I B-Raf inhibitor vemurafenib and the type II inhibitor sorafenib stabilized the kinase domain (KD) but had distinct effects on the Ras-binding domain. Stabilization of the B-Raf(WT) KD was confirmed by hydrogen/deuterium exchange MS and molecular dynamics simulations. Our results are further supported by cellular assays in which we assessed cell viability and phosphorylation profiles in cells expressing B-Raf(WT) or B-Raf(V600E) in response to vemurafenib or sorafenib. Our data indicate that an overall stabilization of the B-Raf structure by specific inhibitors activates MAPK signaling and increases cell survival, helping to explain clinical treatment failure. We also applied our method to monitor conformational changes upon nucleotide binding of the pseudokinase KSR1, which holds high potential for inhibition in human diseases.

Measurement of State-Specific Association Constants in Allosteric Sensors through Molecular Stapling and NMR

Allostery is a ubiquitous mechanism to control biological function and arises from the coupling of inhibitory and binding equilibria. The extent of coupling reflects the inactive vs active state selectivity of the allosteric effector. Hence, dissecting allosteric determinants requires quantification of state-specific association constants. However, observed association constants are typically population-averages, reporting on overall affinities but not on allosteric coupling. Here we propose a general method to measure state-specific association constants in allosteric sensors based on three key elements, i.e., state-selective molecular stapling through disulfide bridges, competition binding saturation transfer experiments and chemical shift correlation analyses to gauge state populations. The proposed approach was applied to the prototypical cyclic adenosine monophosphate (cAMP)-dependent protein kinase (PKA-RIα), for which the structures of the inactive and active states are available, as needed to design the state-selective disulfide bridges. Surprisingly, the PKA-RIα state-specific association constants are comparable to those of a structurally homologous domain with ∼10(3)-fold lower cAMP-affinity, suggesting that the affinity difference arises primarily from changes in the position of the dynamic apo inhibitory equilibrium.

Parallel Allostery by cAMP and PDE Coordinates Activation and Termination Phases in cAMP Signaling

The second messenger molecule cAMP regulates the activation phase of the cAMP signaling pathway through high-affinity interactions with the cytosolic cAMP receptor, the protein kinase A regulatory subunit (PKAR). Phosphodiesterases (PDEs) are enzymes responsible for catalyzing hydrolysis of cAMP to 5' AMP. It was recently shown that PDEs interact with PKAR to initiate the termination phase of the cAMP signaling pathway. While the steps in the activation phase are well understood, steps in the termination pathway are unknown. Specifically, the binding and allosteric networks that regulate the dynamic interplay between PKAR, PDE, and cAMP are unclear. In this study, PKAR and PDE from Dictyostelium discoideum (RD and RegA, respectively) were used as a model system to monitor complex formation in the presence and absence of cAMP. Amide hydrogen/deuterium exchange mass spectrometry was used to monitor slow conformational transitions in RD, using disordered regions as conformational probes. Our results reveal that RD regulates its interactions with cAMP and RegA at distinct loci by undergoing slow conformational transitions between two metastable states. In the presence of cAMP, RD and RegA form a stable ternary complex, while in the absence of cAMP they maintain transient interactions. RegA and cAMP each bind at orthogonal sites on RD with resultant contrasting effects on its dynamics through parallel allosteric relays at multiple important loci. RD thus serves as an integrative node in cAMP termination by coordinating multiple allosteric relays and governing the output signal response.

Allostery through the computational microscope: cAMP activation of a canonical signalling domain

Ligand-induced protein allostery plays a central role in modulating cellular signalling pathways. Here using the conserved cyclic nucleotide-binding domain of protein kinase A's (PKA) regulatory subunit as a prototype signalling unit, we combine long-timescale, all-atom molecular dynamics simulations with Markov state models to elucidate the conformational ensembles of PKA's cyclic nucleotide-binding domain A for the cAMP-free (apo) and cAMP-bound states. We find that both systems exhibit shallow free-energy landscapes that link functional states through multiple transition pathways. This observation suggests conformational selection as the general mechanism of allostery in this canonical signalling domain. Further, we expose the propagation of the allosteric signal through key structural motifs in the cyclic nucleotide-binding domain and explore the role of kinetics in its function. Our approach integrates disparate lines of experimental data into one cohesive framework to understand structure, dynamics and function in complex biological systems.

Probing allosteric mechanisms using native mass spectrometry

Native mass spectrometry (MS) and ion mobility MS provide a way to discriminate between various allosteric mechanisms that cannot be distinguished using ensemble measurements of ligand binding in bulk protein solutions. Native MS, which yields mass measurements of intact assemblies, can be used to determine the values of ligand binding constants of multimeric allosteric proteins, thereby providing a way to distinguish, for example, between concerted and sequential allosteric models. Native MS can also be employed to study cooperativity owing to ligand-modulated protein oligomerization. The rotationally averaged cross-section areas of complexes obtained by ion mobility MS can be used to distinguish between induced fit and conformational selection. Native MS and its allied techniques are, therefore, becoming increasingly powerful tools for dissecting allosteric mechanisms.

Noncognate DNA damage prevents the formation of the active conformation of the Y-family DNA polymerases DinB and DNA polymerase κ

Y-family DNA polymerases are specialized to copy damaged DNA, and are associated with increased mutagenesis, owing to their low fidelity. It is believed that the mechanism of nucleotide selection by Y-family DNA polymerases involves conformational changes preceding nucleotidyl transfer, but there is limited experimental evidence for such structural changes. In particular, nucleotide-induced conformational changes in bacterial or eukaryotic Y-family DNA polymerases have, to date, not been extensively characterized. Using hydrogen-deuterium exchange mass spectrometry, we demonstrate here that the Escherichia coli Y-family DNA polymerase DinB and its human ortholog DNA polymerase κ undergo a conserved nucleotide-induced conformational change in the presence of undamaged DNA and the correct incoming nucleotide. Notably, this holds true for damaged DNA containing N(2) -furfuryl-deoxyguanosine, which is efficiently copied by these two polymerases, but not for damaged DNA containing the major groove modification O(6) -methyl-deoxyguanosine, which is a poor substrate. Our observations suggest that DinB and DNA polymerase κ utilize a common mechanism for nucleotide selection involving a conserved prechemical conformational transition promoted by the correct nucleotide and only preferred DNA substrates.

Cyclic nucleotide binding and structural changes in the isolated GAF domain of Anabaena adenylyl cyclase, CyaB2

GAF domains are a large family of regulatory domains, and a subset are found associated with enzymes involved in cyclic nucleotide (cNMP) metabolism such as adenylyl cyclases and phosphodiesterases. CyaB2, an adenylyl cyclase from Anabaena, contains two GAF domains in tandem at the N-terminus and an adenylyl cyclase domain at the C-terminus. Cyclic AMP, but not cGMP, binding to the GAF domains of CyaB2 increases the activity of the cyclase domain leading to enhanced synthesis of cAMP. Here we show that the isolated GAFb domain of CyaB2 can bind both cAMP and cGMP, and enhanced specificity for cAMP is observed only when both the GAFa and the GAFb domains are present in tandem (GAFab domain). In silico docking and mutational analysis identified distinct residues important for interaction with either cAMP or cGMP in the GAFb domain. Structural changes associated with ligand binding to the GAF domains could not be detected by bioluminescence resonance energy transfer (BRET) experiments. However, amide hydrogen-deuterium exchange mass spectrometry (HDXMS) experiments provided insights into the structural basis for cAMP-induced allosteric regulation of the GAF domains, and differences in the changes induced by cAMP and cGMP binding to the GAF domain. Thus, our findings could allow the development of molecules that modulate the allosteric regulation by GAF domains present in pharmacologically relevant proteins.

Identification of novel cyclic nucleotide binding proteins in Trypanosoma cruzi

Cyclic AMP has been implicated as second messenger in a wide range of cellular processes. In the protozoan parasite Trypanosoma cruzi, cAMP is involved in the development of the parasite's life cycle. While cAMP effectors have been widely studied in other eukaryotic cells, little is known about cAMP's mechanism of action in T. cruzi. To date, only a cAMP-dependent protein kinase A (PKA) has been cloned and characterised in this parasite; however experimental evidence indicates the existence of cAMP-dependent, PKA-independent events. In order to identify new cAMP binding proteins as potential cAMP effectors, we carried out in silico studies using the predicted T. cruzi proteome. Using a combination of search methods 27 proteins with putative cNMP binding domains (CBDs) were identified. Phylogenetic analysis of the CBDs presented a homogeneous distribution, with sequences segregated into two main branches: one containing kinases-like proteins and the other gathering hypothetical proteins with different function or no other known. Comparative modelling of the strongest candidates provides support for the hypothesis that these proteins may give rise to structurally viable cyclic nucleotide binding domains. Pull-down and nucleotide displacement assays strongly suggest that TcCLB.508523.80 could bind cAMP and eventually be a new putative PKA-independent cAMP effector in T. cruzi.

Active site coupling in PDE:PKA complexes promotes resetting of mammalian cAMP signaling

Cyclic 3'5' adenosine monophosphate (cAMP)-dependent-protein kinase (PKA) signaling is a fundamental regulatory pathway for mediating cellular responses to hormonal stimuli. The pathway is activated by high-affinity association of cAMP with the regulatory subunit of PKA and signal termination is achieved upon cAMP dissociation from PKA. Although steps in the activation phase are well understood, little is known on how signal termination/resetting occurs. Due to the high affinity of cAMP to PKA (KD ∼ low nM), bound cAMP does not readily dissociate from PKA, thus begging the question of how tightly bound cAMP is released from PKA to reset its signaling state to respond to subsequent stimuli. It has been recently shown that phosphodiesterases (PDEs) can catalyze dissociation of bound cAMP and thereby play an active role in cAMP signal desensitization/termination. This is achieved through direct interactions with the regulatory subunit of PKA, thereby facilitating cAMP dissociation and hydrolysis. In this study, we have mapped direct interactions between a specific cyclic nucleotide phosphodiesterase (PDE8A) and a PKA regulatory subunit (RIα isoform) in mammalian cAMP signaling, by a combination of amide hydrogen/deuterium exchange mass spectrometry, peptide array, and computational docking. The interaction interface of the PDE8A:RIα complex, probed by peptide array and hydrogen/deuterium exchange mass spectrometry, brings together regions spanning the phosphodiesterase active site and cAMP-binding sites of RIα. Computational docking combined with amide hydrogen/deuterium exchange mass spectrometry provided a model for parallel dissociation of bound cAMP from the two tandem cAMP-binding domains of RIα. Active site coupling suggests a role for substrate channeling in the PDE-dependent dissociation and hydrolysis of cAMP bound to PKA. This is the first instance, to our knowledge, of PDEs directly interacting with a cAMP-receptor protein in a mammalian system, and highlights an entirely new class of binding partners for RIα. This study also highlights applications of structural mass spectrometry combined with computational docking for mapping dynamics in transient signaling protein complexes. Together, these results present a novel and critical role for phosphodiesterases in moderating local concentrations of cAMP in microdomains and signal resetting.

A tool set to map allosteric networks through the NMR chemical shift covariance analysis

Allostery is an essential regulatory mechanism of biological function. Allosteric sites are also pharmacologically relevant as they are often targeted with higher selectivity than orthosteric sites. However, a comprehensive map of allosteric sites poses experimental challenges because allostery is driven not only by structural changes, but also by modulations in dynamics that typically remain elusive to classical structure determination methods. An avenue to overcome these challenges is provided by the NMR chemical shift covariance analysis (CHESCA), as chemical shifts are exquisitely sensitive to redistributions in dynamic conformational ensembles. Here, we propose a set of complementary CHESCA algorithms designed to reliably detect allosteric networks with minimal occurrences of false positives or negatives. The proposed CHESCA toolset was tested for two allosteric proteins (PKA and EPAC) and is expected to complement traditional comparative structural analyses in the comprehensive identification of functionally relevant allosteric sites, including those in otherwise elusive partially unstructured regions.

Using Markov state models to develop a mechanistic understanding of protein kinase A regulatory subunit RIα activation in response to cAMP binding

Protein kinase A (PKA) holoenzyme consists of two catalytic (C) subunits and a regulatory (R) subunit dimer (R(2)C(2)). The kinase is activated by the binding of cAMPs to the two cyclic nucleotide binding domains (CBDs), A and B, on each R-subunit. Despite extensive study, details of the allosteric mechanisms underlying the cooperativity of holoenzyme activation remain unclear. Several Markov state models of PKA-RIα were developed to test competing theories of activation for the R(2)C(2) complex. We found that CBD-B plays an essential role in R-C interaction and promotes the release of the first C-subunit prior to the binding to CBD-A. This favors a conformational selection mechanism for release of the first C-subunit of PKA. However, the release of the second C-subunit requires all four cAMP sites to be occupied. These analyses elucidate R-C heterodimer interactions in the cooperative activation of PKA and cAMP binding and represent a new mechanistic model of R(2)C(2) PKA-RIα activation.

A mechanism for the auto-inhibition of hyperpolarization-activated cyclic nucleotide-gated (HCN) channel opening and its relief by cAMP

Hyperpolarization-activated cyclic nucleotide-gated (HCN) ion channels control neuronal and cardiac electrical rhythmicity. There are four homologous isoforms (HCN1-4) sharing a common multidomain architecture that includes an N-terminal transmembrane tetrameric ion channel followed by a cytoplasmic "C-linker," which connects a more distal cAMP-binding domain (CBD) to the inner pore. Channel opening is primarily stimulated by transmembrane elements that sense membrane hyperpolarization, although cAMP reduces the voltage required for HCN activation by promoting tetramerization of the intracellular C-linker, which in turn relieves auto-inhibition of the inner pore gate. Although binding of cAMP has been proposed to relieve auto-inhibition by affecting the structure of the C-linker and CBD, the nature and extent of these cAMP-dependent changes remain limitedly explored. Here, we used NMR to probe the changes caused by the binding of cAMP and of cCMP, a partial agonist, to the apo-CBD of HCN4. Our data indicate that the CBD exists in a dynamic two-state equilibrium, whose position as gauged by NMR chemical shifts correlates with the V½ voltage measured through electrophysiology. In the absence of cAMP, the most populated CBD state leads to steric clashes with the activated or "tetrameric" C-linker, which becomes energetically unfavored. The steric clashes of the apo tetramer are eliminated either by cAMP binding, which selects for a CBD state devoid of steric clashes with the tetrameric C-linker and facilitates channel opening, or by a transition of apo-HCN to monomers or dimer of dimers, in which the C-linker becomes less structured, and channel opening is not facilitated.

Unraveling structural mechanisms of allosteric drug action

Orthosteric drugs block the active site to obstruct function; allosteric drugs modify the population of the active state, to modulate function. Available data lead us to propose that allosteric drugs can constitute anchors and drivers. The anchor docks into an allosteric pocket. The conformation with which it interacts is unchanged during the transition between the inactive and active states. The anchor provides the foundation that allows the driver to exert a 'pull' and/or 'push' action that shifts the receptor population from the inactive to the active state. The presence or absence of driver atom in an allosteric drug can exert opposite agonism. We map a strategy for driver identification and expect the allosteric trigger concept to transform agonist/antagonist drug discovery.

Role of arginine 209 in the conformational transition of the protein kinase A regulatory subunit RIα A-domain

Using the combination of molecular dynamics (MD) simulations and geometric clustering we analyzed the role of arginine at 209 position in the transition of protein kinase A Iα (PKA Iα) regulatory subunit A-domain from H- to B-conformation and stabilization of the latter. The mechanism underlying the role of the residue at position 209 in the realization of B-conformation includes: (1) possibility to bind the ligand tightly (if transition happens in the presence of cAMP), (2) capability to hold β2β3-loop in the correct conformation, (3) tendency of residue at 209 position to stabilize B-conformation in the absence and in presence of the ligand. In terms of the effect produced on transition of A-domain from H- to B-conformation in the presence of cAMP, mutational substitutions for R209 can be arranged in the following order: Glu(Gly)>Lys>Ile. In the absence of cAMP the order is different Lys>Gly>Glu>Ile. Thus, our results allow us to presume that the role of arginine at 209 position can be important though not crucial.

HDX-MS takes centre stage at unravelling kinase dynamics

In recent years, HDX-MS (hydrogen–deuterium exchange coupled to MS) on biomolecules has evolved from a niche technique to a powerful method in the investigation of protein dynamics. Protein kinases, in particular, represent an area of active study using this technique owing to their well-characterized protein structures and their relevance to diseases such as cancer, immune disorders and neurodegenerative defects. In the present review, we describe how HDX-MS has revealed important dynamic properties of protein kinases and provided insight into the mechanisms of drug binding.

PKA RIα homodimer structure reveals an intermolecular interface with implications for cooperative cAMP binding and Carney complex disease

The regulatory (R) subunit is the cAMP receptor of protein kinase A. Following cAMP binding, the inactive PKA holoenzyme complex separates into two active catalytic (C) subunits and a cAMP-bound R dimer. Thus far, only monomeric R structures have been solved, which fell short in explaining differences of cAMP binding for the full-length protein as compared to the truncated R subunits. Here we solved a full-length R-dimer structure that reflects the biologically relevant conformation, and this structure agrees well with small angle X-ray scattering. An isoform-specific interface is revealed between the protomers. This interface acts as an intermolecular sensor for cAMP and explains the cooperative character of cAMP binding to the RIα dimer. Mutagenesis of residues on this interface not only leads to structural and biochemical changes, but is also linked to Carney complex disease.