Rajeevan Selvaratnam

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.

Signaling through dynamic linkers as revealed by PKA

Protein kinase A (PKA) is a prototype of multidomain signaling proteins functioning as allosteric conformational switches. Allosteric transitions have been the subject of extensive structural and dynamic investigations focusing mainly on folded domains. However, the current understanding of the allosteric role of partially unstructured linkers flanking globular domains is limited. Here, we show that a dynamic linker in the regulatory subunit (R) of PKA serves not only as a passive covalent thread, but also as an active allosteric element that controls activation of the kinase subunit (C) by tuning the inhibitory preequilibrium of a minimally populated intermediate (apo R). Apo R samples both C-binding competent (inactive) and incompetent (active) conformations within a nearly degenerate free-energy landscape and such degeneracy maximally amplifies the response to weak (∼2RT), but conformation-selective interactions elicited by the linker. Specifically, the R linker that in the R:C complex docks in the active site of C in apo R preferentially interacts with the C-binding incompetent state of the adjacent cAMP-binding domain (CBD). These unanticipated findings imply that the formation of the intermolecular R:C inhibitory interface occurs at the expense of destabilizing the intramolecular linker/CBD interactions in R. A direct implication of this model, which was not predictable solely based on protein structure, is that the disruption of a linker/CBD salt bridge in the R:C complex unexpectedly leads to increased affinity of R for C. The linker includes therefore sites of R:C complex frustration and frustration-relieving mutations enhance the kinase inhibitory potency of R without compromising its specificity.

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.

Dynamically driven ligand selectivity in cyclic nucleotide binding domains

One of the mechanisms that minimize the aberrant cross-talk between cAMP- and cGMP-dependent signaling pathways relies on the selectivity of cAMP binding domains (CBDs). For instance, the CBDs of two critical eukaryotic cAMP receptors, i.e. protein kinase A (PKA) and the exchange protein activated by cAMP (EPAC), are both selectively activated by cAMP. However, the mechanisms underlying their cAMP versus cGMP selectivity are quite distinct. In PKA this selectivity is controlled mainly at the level of ligand affinity, whereas in EPAC it is mostly determined at the level of allostery. Currently, the molecular basis for these different selectivity mechanisms is not fully understood. We have therefore comparatively analyzed by NMR the cGMP-bound states of the essential CBDs of PKA and EPAC, revealing key differences between them. Specifically, cGMP binds PKA preserving the same syn base orientation as cAMP at the price of local steric clashes, which lead to a reduced affinity for cGMP. Unlike PKA, cGMP is recognized by EPAC in an anti conformation and generates several short and long range perturbations. Although these effects do not alter significantly the structure of the EPAC CBD investigated, remarkable differences in dynamics between the cAMP- and cGMP-bound states are detected for the ionic latch region. These observations suggest that one of the determinants of cGMP antagonism in EPAC is the modulation of the entropic control of inhibitory interactions and illustrate the pivotal role of allostery in determining signaling selectivity as a function of dynamic changes, even in the absence of significant affinity variations.

Entropy-driven cAMP-dependent Allosteric Control of Inhibitory Interactions in Exchange Proteins Directly Activated by cAMP

Exchange proteins directly activated by cAMP (EPACs) are guanine nucleotide-exchange factors for the small GTPases Rap1 and Rap2 and represent a key receptor for the ubiquitous cAMP second messenger in eukaryotes. The cAMP-dependent activation of apoEPAC is typically rationalized in terms of a preexisting equilibrium between inactive and active states. Structural and mutagenesis analyses have shown that one of the critical determinants of the EPAC activation equilibrium is a cluster of salt bridges formed between the catalytic core and helices α1 and α2 at the N terminus of the cAMP binding domain and commonly referred to as ionic latch (IL). The IL stabilizes the inactive states in a closed topology in which access to the catalytic domain is sterically occluded by the regulatory moiety. However, it is currently not fully understood how the IL is allosterically controlled by cAMP. Chemical shift mapping studies consistently indicate that cAMP does not significantly perturb the structure of the IL spanning sites within the regulatory region, pointing to cAMP-dependent dynamic modulations as a key allosteric carrier of the cAMP-signal to the IL sites. Here, we have therefore investigated the dynamic profiles of the EPAC1 cAMP binding domain in its apo, cAMP-bound, and Rp-cAMPS phosphorothioate antagonist-bound forms using several 15N relaxation experiments. Based on the comparative analysis of dynamics in these three states, we have proposed a model of EPAC activation that incorporates the dynamic features allosterically modulated by cAMP and shows that cAMP binding weakens the IL by increasing its entropic penalty due to dynamic enhancements.

The Auto-Inhibitory Role of the EPAC Hinge Helix as Mapped by NMR

The cyclic-AMP binding domain (CBD) is the central regulatory unit of exchange proteins activated by cAMP (EPAC). The CBD maintains EPAC in a state of auto-inhibition in the absence of the allosteric effector, cAMP. When cAMP binds to the CBD such auto-inhibition is released, leading to EPAC activation. It has been shown that a key feature of such cAMP-dependent activation process is the partial destabilization of a structurally conserved hinge helix at the C-terminus of the CBD. However, the role of this helix in auto-inhibition is currently not fully understood. Here we utilize a series of progressive deletion mutants that mimic the hinge helix destabilization caused by cAMP to show that such helix is also a pivotal auto-inhibitory element of apo-EPAC. The effect of the deletion mutations on the auto-inhibitory apo/inactive vs. apo/active equilibrium was evaluated using recently developed NMR chemical shift projection and covariance analysis methods. Our results show that, even in the absence of cAMP, the C-terminal region of the hinge helix is tightly coupled to other conserved allosteric structural elements of the CBD and perturbations that destabilize the hinge helix shift the auto-inhibitory equilibrium toward the apo/active conformations. These findings explain the apparently counterintuitive observation that cAMP binds more tightly to shorter than longer EPAC constructs. These results are relevant for CBDs in general and rationalize why substrates sensitize CBD-containing systems to cAMP. Furthermore, the NMR analyses presented here are expected to be generally useful to quantitatively evaluate how mutations affect conformational equilibria.

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

Background: Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels control electrical activity through tetramerization of an intracellular linker. Results: NMR shows that the apo-cAMP-binding domain (CBD) of HCN4 destabilizes the tetramer through steric clashes. Conclusion: The apo-HCN4 CBD structure is compatible with monomeric and dimeric but not with tetrameric HCN4. Significance: The proposed mechanism explains HCN auto-inhibition and its relaxation 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.

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.

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.