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.

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.

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.

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 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.

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.

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.

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.

Role of Dynamics in the Autoinhibition and Activation of the 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 allosterically modulates HCN through the cAMP-dependent formation of a tetrameric gating ring spanning the intracellular region (IR) of HCN, to which cAMP binds. Although the apo versus holo conformational changes of the cAMP-binding domain (CBD) have been previously mapped, only limited information is currently available on the HCN IR dynamics, which have been hypothesized to play a critical role in the cAMP-dependent gating of HCN. Here, using molecular dynamics simulations validated and complemented by experimental NMR and CD data, we comparatively analyze HCN IR dynamics in the four states of the thermodynamic cycle arising from the coupling between cAMP binding and tetramerization equilibria. This extensive set of molecular dynamics trajectories captures the active-to-inactive transition that had remained elusive for other CBDs, and it provides unprecedented insight on the role of IR dynamics in HCN autoinhibition and its release by cAMP. Specifically, the IR tetramerization domain becomes more flexible in the monomeric states, removing steric clashes that the apo-CDB structure would otherwise impose. Furthermore, the simulations reveal that the active/inactive structural transition for the apo-monomeric CBD occurs through a manifold of pathways that are more divergent than previously anticipated. Upon cAMP binding, these pathways become disallowed, pre-confining the CBD conformational ensemble to a tetramer-compatible state. This conformational confinement primes the IR for tetramerization and thus provides a model of how cAMP controls HCN channel gating.

Mechanism of cAMP Partial Agonism in Protein Kinase G (PKG).

Protein kinase G (PKG) is a major receptor of cGMP and controls signaling pathways often distinct from those regulated by cAMP. Hence, the selective activation of PKG by cGMP versus cAMP is critical. However, the mechanism of cGMP-versus-cAMP selectivity is only limitedly understood. Although the C-terminal cyclic nucleotide-binding domain B of PKG binds cGMP with higher affinity than cAMP, the intracellular concentrations of cAMP are typically higher than those of cGMP, suggesting that the cGMP-versus-cAMP selectivity of PKG is not controlled uniquely through affinities. Here, we show that cAMP is a partial agonist for PKG, and we elucidate the mechanism for cAMP partial agonism through the comparative NMR analysis of the apo, cGMP-, and cAMP-bound forms of the PKG cyclic nucleotide-binding domain B. We show that although cGMP activation is adequately explained by a two-state conformational selection model, the partial agonism of cAMP arises from the sampling of a third, partially autoinhibited state.

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.

Allosteric linkers in cAMP signalling.

Weak interactions mediated by dynamic linkers are key determinants of allosteric regulation in multidomain signalling proteins. However, the mechanisms of linker-dependent control have remained largely elusive. In the present article, we review an allosteric model introduced recently to explain how signalling proteins effectively sense and respond to weak interactions, such as those elicited by flexible linkers flanking globular domains. Central to this model is the idea that near degeneracy within the free energy landscape of conformational selection maximally amplifies the response to weak (~2RT), but conformation-selective interactions. The model was tested as proof of principle using the prototypical regulatory subunit (R) of protein kinase A and led to the unanticipated finding that dynamic linkers control kinase activation and inhibition by tuning the inhibitory pre-equilibrium of a minimally populated intermediate (apo R). A practical implication of the proposed model is a new strategy to design kinase inhibitors with enhanced potency through frustration-relieving mutations.

Tapping the translation potential of cAMP signalling: molecular basis for selectivity in cAMP agonism and antagonism as revealed by NMR.

Eukaryotic CBDs (cAMP-binding domains) control multiple cellular functions (e.g. phosphorylation, guanine exchange and ion channel gating). Hence the manipulation of cAMP-dependent signalling pathways has a high translational potential. However, the ubiquity of eukaryotic CBDs also poses a challenge in terms of selectivity. Before the full translational potential of cAMP signalling can be tapped, it is critical to understand the structural basis for selective cAMP agonism and antagonism. Recent NMR investigations have shown that structurally homologous CBDs respond differently to several CBD ligands and that these unexpected differences arise at the level of either binding (i.e. affinity) or allostery (i.e. modulation of the autoinhibitory equilibria). In the present article, we specifically address how the highly conserved CBD fold binds cAMP with markedly different affinities in PKA (protein kinase A) relative to other eukaryotic cAMP receptors, such as Epac (exchange protein directly activated by cAMP) and HCN (hyperpolarization-activated cyclic-nucleotide-modulated channel). A major emerging determinant of cAMP affinity is hypothesized to be the position of the autoinhibitory equilibrium of the apo-CBD, which appears to vary significantly across different CBDs. These analyses may assist the development of selective CBD effectors that serve as potential drug leads for the treatment of cardiovascular diseases.

Mapping allostery through the covariance analysis of NMR chemical shifts.

Allostery is a fundamental mechanism of regulation in biology. The residues at the end points of long-range allosteric perturbations are commonly identified by the comparative analyses of structures and dynamics in apo and effector-bound states. However, the networks of interactions mediating the propagation of allosteric signals between the end points often remain elusive. Here we show that the covariance analysis of NMR chemical shift changes caused by a set of covalently modified analogs of the allosteric effector (i.e., agonists and antagonists) reveals extended networks of coupled residues. Unexpectedly, such networks reach not only sites subject to effector-dependent structural variations, but also regions that are controlled by dynamically driven allostery. In these regions the allosteric signal is propagated mainly by dynamic rather than structural modulations, which result in subtle but highly correlated chemical shift variations. The proposed chemical shift covariance analysis (CHESCA) identifies interresidue correlations based on the combination of agglomerative clustering (AC) and singular value decomposition (SVD). AC results in dendrograms that define functional clusters of coupled residues, while SVD generates score plots that provide a residue-specific dissection of the contributions to binding and allostery. The CHESCA approach was validated by applying it to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC) and the CHESCA results are in full agreement with independent mutational data on EPAC activation. Overall, CHESCA is a generally applicable method that utilizes a selected chemical library of effector analogs to quantitatively decode the binding and allosteric information content embedded in chemical shift changes.

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?

The projection analysis of NMR chemical shifts reveals extended EPAC autoinhibition determinants.

EPAC is a cAMP-dependent guanine nucleotide exchange factor that serves as a prototypical molecular switch for the regulation of essential cellular processes. Although EPAC activation by cAMP has been extensively investigated, the mechanism of EPAC autoinhibition is still not fully understood. The steric clash between the side chains of two conserved residues, L273 and F300 in EPAC1, has been previously shown to oppose the inactive-to-active conformational transition in the absence of cAMP. However, it has also been hypothesized that autoinhibition is assisted by entropic losses caused by quenching of dynamics that occurs if the inactive-to-active transition takes place in the absence of cAMP. Here, we test this hypothesis through the comparative NMR analysis of several EPAC1 mutants that target different allosteric sites of the cAMP-binding domain (CBD). Using what to our knowledge is a novel projection analysis of NMR chemical shifts to probe the effect of the mutations on the autoinhibition equilibrium of the CBD, we find that whenever the apo/active state is stabilized relative to the apo/inactive state, dynamics are consistently quenched in a conserved loop (ß2-ß3) and helix (α5) of the CBD. Overall, our results point to the presence of conserved and nondegenerate determinants of CBD autoinhibition that extends beyond the originally proposed L273/F300 residue pair, suggesting that complete activation necessitates the simultaneous suppression of multiple autoinhibitory mechanisms, which in turn confers added specificity for the cAMP allosteric effector.

Role of dynamics in the autoinhibition and activation of the exchange protein directly activated by cyclic AMP (EPAC).

The exchange protein directly activated by cAMP (EPAC) is a key receptor of cAMP in eukaryotes and controls critical signaling pathways. Currently, no residue resolution information is available on the full-length EPAC dynamics, which are known to be pivotal determinants of allostery. In addition, no information is presently available on the intermediates for the classical induced fit and conformational selection activation pathways. Here these questions are addressed through molecular dynamics simulations on five key states along the thermodynamic cycle for the cAMP-dependent activation of a fully functional construct of EPAC2, which includes the cAMP-binding domain and the integral catalytic region. The simulations are not only validated by the agreement with the experimental trends in cAMP-binding domain dynamics determined by NMR, but they also reveal unanticipated dynamic attributes, rationalizing previously unexplained aspects of EPAC activation and autoinhibition. Specifically, the simulations show that cAMP binding causes an extensive perturbation of dynamics in the distal catalytic region, assisting the recognition of the Rap1b substrate. In addition, analysis of the activation intermediates points to a possible hybrid mechanism of EPAC allostery incorporating elements of both the induced fit and conformational selection models. In this mechanism an entropy compensation strategy results in a low free-energy pathway of activation. Furthermore, the simulations indicate that the autoinhibitory interactions of EPAC are more dynamic than previously anticipated, leading to a revised model of autoinhibition in which dynamics fine tune the stability of the autoinhibited state, optimally sensitizing it to cAMP while avoiding constitutive activation.

cAMP-dependent allostery and dynamics in Epac: an NMR view.

Epac (exchange protein directly activated by cAMP) is a critical cAMP receptor, which senses cAMP and couples the cAMP signal to the catalysis of guanine exchange in the Rap substrate. In the present paper, we review the NMR studies that we have undertaken on the CBD (cyclic-nucleotide-binding domain) of Epac1. Our NMR investigations have shown that cAMP controls distal autoinhibitory interactions through long-range modulations in dynamics. Such dynamically mediated allosteric effects contribute not only to the cAMP-dependent activation of Epac, but also to the selectivity of Epac for cAMP in contrast with cGMP. In addition, we have mapped the interaction networks that couple the cAMP-binding site to the sites involved in the autoinhibitory interactions, using a method based on the covariance analysis of NMR chemical shifts. We anticipate that this approach is generally applicable to dissect allosteric networks in signalling domains.

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.