Probing ligand selectivity in pathogens.

Why does protein kinase A respond to purine nucleosides in certain pathogens, but not to the cyclic nucleotides that activate this kinase in most other organisms?

Fractionation factors reveal hidden frustration in an ancient allosteric module.

Protein kinase G (PKG) is an essential regulator of eukaryotic cyclic guanosine monophosphate (cGMP)-dependent intracellular signaling, controlling pathways that are often distinct from those regulated by cyclic adenosine monophosphate (cAMP). Specifically, the C-terminal cyclic-nucleotide-binding domain (CNB-B) of PKG has emerged as a critical module to control allostery and cGMP-selectivity in PKG. While key contributions to the cGMP-versus-cAMP selectivity of CNB-B were previously assessed, only limited knowledge is currently available on how cyclic nucleotide binding rewires the network of hydrogen bonds in CNB-B, and how such rewiring contributes to allostery and cGMP selectivity. To address this gap, we extend the comparative analysis of apo, cAMP- and cGMP-bound CNB-B to H/D fractionation factors (FFs), which are well-suited for assessing backbone hydrogen-bond strengths within proteins. Apo-vs-bound comparisons inform of perturbations arising from both binding and allostery, while cGMP-bound vs cAMP-bound comparisons inform of perturbations that are purely allosteric. The comparative FF analyses of the bound states revealed mixed patterns of hydrogen-bond strengthening and weakening, pointing to inherent frustration, whereby not all hydrogen bonds can be simultaneously stabilized. Interestingly, contrary to expectations, these patterns include a weakening of hydrogen bonds not only within critical recognition and allosteric elements of CNB-B, but also within elements known to undergo rigid-body movement upon cyclic nucleotide binding. These results suggest that frustration may contribute to the reversibility of allosteric conformational shifts by avoiding over-rigidification that may otherwise trap CNB-B in its active state. Considering that PKG CNB-B serves as a prototype for allosteric conformational switches, similar concepts may be applicable to allosteric domains in general.

An auto-inhibited state of protein kinase G and implications for selective activation.

Cyclic GMP-dependent protein kinases (PKGs) are key mediators of the nitric oxide/cyclic guanosine monophosphate (cGMP) signaling pathway that regulates biological functions as diverse as smooth muscle contraction, cardiac function, and axon guidance. Understanding how cGMP differentially triggers mammalian PKG isoforms could lead to new therapeutics that inhibit or activate PKGs, complementing drugs that target nitric oxide synthases and cyclic nucleotide phosphodiesterases in this signaling axis. Alternate splicing of PRKG1 transcripts confers distinct leucine zippers, linkers, and auto-inhibitory (AI) pseudo-substrate sequences to PKG Iα and Iß that result in isoform-specific activation properties, but the mechanism of enzyme auto-inhibition and its alleviation by cGMP is not well understood. Here, we present a crystal structure of PKG Iß in which the AI sequence and the cyclic nucleotide-binding (CNB) domains are bound to the catalytic domain, providing a snapshot of the auto-inhibited state. Specific contacts between the PKG Iß AI sequence and the enzyme active site help explain isoform-specific activation constants and the effects of phosphorylation in the linker. We also present a crystal structure of a PKG I CNB domain with an activating mutation linked to Thoracic Aortic Aneurysms and Dissections. Similarity of this structure to wildtype cGMP-bound domains and differences with the auto-inhibited enzyme provide a mechanistic basis for constitutive activation. We show that PKG Iß auto-inhibition is mediated by contacts within each monomer of the native full-length dimeric protein, and using the available structural and biochemical data we develop a model for the regulation and cooperative activation of PKGs.

Identification of core allosteric sites through temperature- and nucleus-invariant chemical shift covariance.

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or in the design of allosteric drugs. Here, we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. T- and CLASS-CHESCAs, as well as complete-linkage CHESCA, were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, T-, CLASS-, and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.

Allosteric pluripotency: challenges and opportunities.

Allosteric pluripotency arises when the functional response of an allosteric receptor to an allosteric stimulus depends on additional allosteric modulators. Here, we discuss allosteric pluripotency as observed in the prototypical Protein Kinase A (PKA) as well as in other signaling systems, from typical multidomain signaling proteins to bacterial enzymes. We identify key drivers of pluripotent allostery and illustrate how hypothesizing allosteric pluripotency may solve apparent discrepancies currently present in the literature regarding the dual nature of known allosteric modulators. We also outline the implications of allosteric pluripotency for cellular signaling and allosteric drug design, and analyze the challenges and opportunities opened by the pluripotent nature of allostery.

Non-Canonical Allostery in Cyclic Nucleotide Dependent Kinases.

The cAMP- and cGMP-dependent protein kinases (PKA and PKG) are canonically activated by the corresponding cyclic nucleotides. However, both systems are also sensitive to a wide range of non-canonical allosteric effectors, such as reactive oxygen species, which induce the formation of regulatory inter- and intra-molecular disulfide bridges, and disease-related mutations (DRMs). Here, we present a combined analysis of representative non-canonical allosteric effectors for PKA and PKG, and we identify common molecular mechanisms underlying non-canonical allostery in these kinases, from shifts in dynamical regulatory equilibria to modulation of inter-protomer interactions. In addition, mutations may also drive oligomerization beyond dimerization, and possibly phase transitions, causing loss of kinase inhibitory function and amplifying the allosteric effects of DRMs. Hence non-canonical allosteric stimuli often result in constitutive kinase activation underlying either physiological control of downstream signaling pathways or pathological outcomes, from aortic aneurisms to cancer predisposition. Overall, PKA and PKG emerge as "pan-sensors" going well beyond canonical cyclic nucleotide activation, revealing their versatile roles as central signaling hubs.

Divergent allostery reveals critical differences between structurally homologous regulatory domains of Plasmodium falciparum and human protein kinase G.

Malaria is a life-threatening infectious disease primarily caused by the Plasmodium falciparum parasite. The increasing resistance to current antimalarial drugs and their side effects has led to an urgent need for novel malaria drug targets, such as the P. falciparum cGMP-dependent protein kinase (pfPKG). However, PKG plays an essential regulatory role also in the human host. Human cGMP-dependent protein kinase (hPKG) and pfPKG are controlled by structurally homologous cGMP-binding domains (CBDs). Here, we show that despite the structural similarities between the essential CBDs in pfPKG and hPKG, their respective allosteric networks differ significantly. Through comparative analyses of chemical shift covariance analyses, molecular dynamics simulations, and backbone internal dynamics measurements, we found that conserved allosteric elements within the essential CBDs are wired differently in pfPKG and hPKG to implement cGMP-dependent kinase activation. Such pfPKG versus hPKG rewiring of allosteric networks was unexpected because of the structural similarity between the two essential CBDs. Yet, such finding provides crucial information on which elements to target for selective inhibition of pfPKG versus hPKG, which may potentially reduce undesired side effects in malaria treatments.

State-selective frustration as a key driver of allosteric pluripotency.

Allosteric pluripotency arises when an allosteric effector switches from agonist to antagonist depending on the experimental conditions. For example, the Rp-cAMPS ligand of Protein Kinase A (PKA) switches from agonist to antagonist as the MgATP concentration increases and/or the kinase substrate affinity or concentration decreases. Understanding allosteric pluripotency is essential to design effective allosteric therapeutics with minimal side effects. Allosteric pluripotency of PKA arises from divergent allosteric responses of two homologous tandem cAMP-binding domains, resulting in a free energy landscape for the Rp-cAMPS-bound PKA regulatory subunit R1a in which the ground state is kinase inhibition-incompetent and the kinase inhibition-competent state is excited. The magnitude of the free energy difference between the ground non-inhibitory and excited inhibitory states (ΔG R,Gap) relative to the effective free energy of R1a binding to the catalytic subunit of PKA (ΔG R:C) dictates whether the antagonism-to-agonism switch occurs. However, the key drivers of ΔG R,Gap are not fully understood. Here, by analyzing an R1a mutant that selectively silences allosteric pluripotency, we show that a major determinant of ΔG R,Gap unexpectedly arises from state-selective frustration in the ground inhibition-incompetent state of Rp-cAMPS-bound R1a. Such frustration is caused by steric clashes between the phosphate-binding cassette and the helices preceding the lid, which interact with the phosphate and base of Rp-cAMPS, respectively. These clashes are absent in the excited inhibitory state, thus reducing the ΔG R,Gap to values comparable to ΔG R:C, as needed for allosteric pluripotency to occur. The resulting model of allosteric pluripotency is anticipated to assist the design of effective allosteric modulators.

Mutual Protein-Ligand Conformational Selection Drives cGMP vs. cAMP Selectivity in Protein Kinase G.

Protein kinase G (PKG) is a major receptor of cGMP, and controls signaling pathways distinct from those regulated by cAMP. However, the contributions of the two substituents that differentiate cGMP from cAMP (i.e. 6-oxo and 2-NH2) to the cGMP-versus-cAMP selectivity of PKG remain unclear. Here, using NMR to map how binding affinity and dynamics of the protein and ligand vary along a ligand double-substitution cycle, we show that the contributions of the two substituents to binding affinity are surprisingly non-additive. Such non-additivity stems primarily from mutual protein-ligand conformational selection, whereby not only does the ligand select for a preferred protein conformation upon binding, but also, the protein selects for a preferred ligand conformation. The 6-oxo substituent mainly controls the conformational equilibrium of the bound protein, while the 2-NH2 substituent primarily controls the conformational equilibrium of the unbound ligand (i.e. syn versus anti). Therefore, understanding the conformational dynamics of both the protein and ligand is essential to explain the cGMP-versus-cAMP selectivity of PKG.

Dynamical Basis of Allosteric Activation for the Plasmodium falciparum Protein Kinase G.

The Plasmodium falciparum cGMP-dependent protein kinase (PfPKG) is required for the progression of the Plasmodium's life cycle and is therefore a promising malaria drug target. PfPKG includes four cGMP-binding domains (CBD-A to CBD-D). CBD-D plays a crucial role in PfPKG regulation as it is the primary determinant for the inhibition and cGMP-dependent activation of the catalytic domain. Hence, it is critical to understand how CBD-D is allosterically regulated by cGMP. Although the apo versus holo conformational changes of CBD-D have been reported, information on the intermediates of the activation pathway is currently lacking. Here, we employed molecular dynamics simulations to model four key states along the thermodynamic cycle for the cGMP-dependent activation of the PfPKG CBD-D domain. The simulations were compared to NMR data, and they revealed that the PfPKG CBD-D activation pathway samples a compact intermediate in which the N- and C-terminal helices approach the central ß-barrel. In addition, by comparing the cGMP-bound active and inactive states, the essential binding interactions that differentiate these states were identified. The identification of structural and dynamical features unique to the cGMP-bound inactive state provides a promising basis to design PfPKG-selective allosteric inhibitors as a viable treatment for malaria.

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.

Allosteric Mechanisms of Nonadditive Substituent Contributions to Protein-Ligand Binding.

Quantifying chemical substituent contributions to ligand-binding free energies is challenging due to nonadditive effects. Protein allostery is a frequent cause of nonadditivity, but the underlying allosteric mechanisms often remain elusive. Here, we propose a general NMR-based approach to elucidate such mechanisms and we apply it to the HCN4 ion channel, whose cAMP-binding domain is an archetypal conformational switch. Using NMR, we show that nonadditivity arises not only from concerted conformational transitions, but also from conformer-specific effects, such as steric frustration. Our results explain how affinity-reducing functional groups may lead to affinity gains if combined. Surprisingly, our approach also reveals that nonadditivity depends markedly on the receptor conformation. It is negligible for the inhibited state but highly significant for the active state, opening new opportunities to tune potency and agonism of allosteric effectors.

The structure of the apo cAMP-binding domain of HCN4 - a stepping stone toward understanding the cAMP-dependent modulation of the hyperpolarization-activated cyclic-nucleotide-gated ion channels.

The hyperpolarization-activated cyclic-nucleotide-gated (HCN) ion channels control nerve impulse transmission and cardiac pacemaker activity. The modulation by cAMP is critical for the regulatory function of HCN in both neurons and cardiomyocytes, but the underlying mechanism is not fully understood. Here, we show how the structure of the apo cAMP-binding domain of the HCN4 isoform has contributed to a model for the cAMP-dependent modulation of the HCN ion-channel. This model recapitulates the structural and dynamical changes that occur along the thermodynamic cycle arising from the coupling of cAMP-binding and HCN self-association equilibria. The proposed model addresses some of the questions previously open about the auto-inhibition of HCN and its cAMP-induced activation, while opening new opportunities for selectively targeting HCN through allosteric ligands. A remaining challenge is the investigation of HCN dimers and their regulatory role. Overcoming this challenge will require the integration of crystallography, cryo electron microscopy, NMR, electrophysiology and simulations.

Role of Dimers in the cAMP-Dependent Activation of Hyperpolarization-Activated Cyclic-Nucleotide-Modulated (HCN) Ion Channels.

Hyperpolarization-activated cyclic-nucleotide-modulated (HCN) ion channels control rhythmicity in neurons and cardiomyocytes. Cyclic AMP (cAMP) modulates HCN activity through the cAMP-induced formation of a tetrameric gating ring spanning the intracellular region (IR) of HCN. Although evidence from confocal patch-clamp fluorometry indicates that the cAMP-dependent gating of HCN occurs through a dimer of dimers, the structural and dynamical basis of cAMP allostery in HCN dimers has so far remained elusive. Thus, here we examine how dimers influence IR structural dynamics, and the role that such structural dynamics play in HCN allostery. To this end, we performed molecular dynamics (MD) simulations of HCN4 IR dimers in their fully apo, fully holo, and partially cAMP-bound states, resulting in a total simulated time of 1.2 µs. Comparative analyses of these MD trajectories, as well as previous monomer and tetramer simulations utilized as benchmarks for comparison, reveal that dimers markedly sensitize the HCN IR to cAMP-modulated allostery. Our results indicate that dimerization fine-tunes the IR dynamics to enhance, relative to both monomers and tetramers, the allosteric intra- and interprotomer coupling between the cAMP-binding domain and tetramerization domain components of the IR. The resulting allosteric model provides a viable rationalization of electrophysiological data on the role of IR dimers in HCN activation.

Functional dynamics in cyclic nucleotide signaling and amyloid inhibition.

It is now established that understanding the molecular basis of biological function requires atomic resolution maps of both structure and dynamics. Here, we review several illustrative examples of functional dynamics selected from our work on cyclic nucleotide signaling and amyloid inhibition. Although fundamentally diverse, a central aspect common to both fields is that function can only be rationalized by considering dynamic equilibria between distinct states of the accessible free energy landscape. The dynamic exchange between ground and excited states of signaling proteins is essential to explain auto-inhibition and allosteric activation. The dynamic exchange between non-toxic monomeric species and toxic oligomers of amyloidogenic proteins provides a foundation to understand amyloid inhibition. NMR ideally probes both types of dynamic exchange at atomic resolution. Specifically, we will show how NMR was utilized to reveal the dynamical basis of cyclic nucleotide affinity, selectivity, agonism and antagonism in multiple eukaryotic cAMP and cGMP receptors. We will also illustrate how NMR revealed the mechanism of action of plasma proteins that act as extracellular chaperones and inhibit the self-association of the prototypical amyloidogenic Aß peptide. The examples outlined in this review illustrate the widespread implications of functional dynamics and the power of NMR as an indispensable tool in molecular pharmacology and pathology.

Molecular Mechanism for the (-)-Epigallocatechin Gallate-Induced Toxic to Nontoxic Remodeling of Aß Oligomers.

(-)-Epigallocatechin gallate (EGCG) effectively reduces the cytotoxicity of the Alzheimer's disease ß-amyloid peptide (Aß) by remodeling seeding-competent Aß oligomers into off-pathway seeding-incompetent Aß assemblies. However, the mechanism of EGCG-induced remodeling is not fully understood. Here we combine 15N and 1H dark-state exchange saturation transfer (DEST), relaxation, and chemical shift projection NMR analyses with fluorescence, dynamic light scattering, and electron microscopy to elucidate how EGCG remodels Aß oligomers. We show that the remodeling adheres to a Hill-Scatchard model whereby the Aß(1-40) self-association occurs cooperatively and generates Aß(1-40) oligomers with multiple independent binding sites for EGCG with a Kd ∼10-fold lower than that for the Aß(1-40) monomers. Upon binding to EGCG, the Aß(1-40) oligomers become less solvent exposed, and the ß-regions, which are involved in direct monomer-protofibril contacts in the absence of EGCG, undergo a direct-to-tethered contact shift. This switch toward less engaged monomer-protofibril contacts explains the seeding incompetency observed upon EGCG remodeling and suggests that EGCG interferes with secondary nucleation events known to generate toxic Aß assemblies. Unexpectedly, the N-terminal residues experience an opposite EGCG-induced shift from tethered to direct contacts, explaining why EGCG remodeling occurs without release of Aß(1-40) monomers. We also show that upon binding Aß(1-40) oligomers the relative positions of the EGCG B and D rings change with respect to that of ring A. These distinct structural changes occurring in both Aß(1-40) oligomers and EGCG during remodeling offer a foundation for understanding the molecular mechanism of EGCG as a neurotoxicity inhibitor. Furthermore, the results reported here illustrate the effectiveness of DEST-based NMR approaches in investigating the mechanism of low-molecular-weight amyloid inhibitors.

Regulation of HCN Ion Channels by Non-canonical Cyclic Nucleotides.

The hyperpolarization-activated cyclic-nucleotide-modulated (HCN) proteins are cAMP-regulated ion channels that play a key role in nerve impulse transmission and heart rate modulation in neuronal and cardiac cells, respectively. Although they are regulated primarily by cAMP, other cyclic nucleotides such as cGMP, cCMP, and cUMP serve as partial agonists for the HCN2 and HCN4 isoforms. By competing with cAMP for binding, these non-canonical ligands alter ion channel gating, and in turn, modulate the cAMP-dependent activation profiles. The partial activation of non-canonical cyclic nucleotides can be rationalized by either a partial reversal of a two-state inactive/active conformational equilibrium, or by sampling of a third conformational state with partial activity. Furthermore, different mechanisms and degrees of activation have been observed upon binding of non-canonical cyclic nucleotides to HCN2 versus HCN4, suggesting that these ligands control HCN ion channels in an isoform-specific manner. While more work remains to be done to achieve a complete understanding of ion channel modulation by non-canonical cyclic nucleotides, it is already clear that such knowledge will ultimately prove invaluable in achieving a more complete understanding of ion channel signaling in vivo, as well as in the development of therapeutics designed to selectively modulate ion channel gating.

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.

Structural Basis of Tonic Inhibition by Dimers of Dimers in Hyperpolarization-Activated Cyclic-Nucleotide-Modulated (HCN) Ion Channels.

The hyperpolarization-activated cyclic-nucleotide-modulated (HCN) ion channels control rhythmicity in neurons and cardiomyocytes. Cyclic AMP (cAMP) modulates HCN activity through the cAMP-dependent formation of a tetrameric gating ring spanning the intracellular region (IR) of HCN. In the absence of cAMP, the IR cAMP-binding domain (CBD) mainly samples its inactive conformation, resulting in steric clashes that destabilize the IR tetramer. Although these clashes with the inactive CBD are released through tetramer dissociation into monomers, functional mutagenesis suggests that the apo IR is not fully monomeric. To investigate the inhibitory nonmonomeric IR species, we performed molecular dynamics simulations starting from "hybrid" structures that are tetrameric but contain inactive apo-state CBD conformations. The ensemble of simulated trajectories reveals that full dissociation of the tetramer into monomers is not necessary to release the steric hindrance with the inactive CBD. Specifically, we found that partial dissociation of the tetramer into dimers is sufficient to accommodate four inactive CBDs, while reduction of the quaternary symmetry of the nondissociated tetramer from 4- to 2-fold permits accommodation of two inactive CBDs. Our findings not only rationalize available electrophysiological, fluorometry, and sedimentation equilibrium data, but also provide unprecedented structural insight into previously elusive nonmonomeric autoinhibitory HCN species.