Dynamics of cAPK type IIbeta activation revealed by enhanced amide H/2H exchange mass spectrometry (DXMS)

cAMP-dependent protein kinase (cAPK) is a key component in numerous cell signaling pathways. The cAPK regulatory (R) subunit maintains the kinase in an inactive state until cAMP saturation of the R-subunit leads to activation of the enzyme. To delineate the conformational changes associated with cAPK activation, the amide hydrogen/deuterium exchange in the cAPK type IIbeta R-subunit was probed by electrospray mass spectrometry. Three states of the R-subunit, cAMP-bound, catalytic (C)-subunit bound, and apo, were incubated in deuterated water for various lengths of time and then, prior to mass spectrometry analysis, subjected to digestion by pepsin to localize the deuterium incorporation. High sequence coverage (>99%) by the pepsin-digested fragments enables us to monitor the dynamics of the whole protein. The effects of cAMP binding on RIIbeta amide hydrogen exchange are restricted to the cAMP-binding pockets, while the effects of C-subunit binding are evident across both cAMP-binding domains and the linker region. The decreased amide hydrogen exchange for residues 253-268 within cAMP binding domain A and for residues 102-115, which include the pseudosubstrate inhibitory site, support the prediction that these two regions represent the conserved primary and peripheral C-subunit binding sites. An increase in amide hydrogen exchange for a broad area within cAMP-binding domain B and a narrow area within cAMP-binding domain A (residues 222-232) suggest that C-subunit binding transmits long-distance conformational changes throughout the protein.     

Mapping intersubunit interactions of the regulatory subunit (RIalpha) in the type I holoenzyme of protein kinase A by amide hydrogen/deuterium exchange mass spectrometry (DXMS)

Protein kinase A (PKA), a central locus for cAMP signaling in the cell, is composed of regulatory (R) and catalytic (C) subunits. The C-subunits are maintained in an inactive state by binding to the R-subunit dimer in a tetrameric holoenzyme complex (R(2)C(2)). PKA is activated by cAMP binding to the R-subunits which induces a conformational change leading to release of the active C-subunit. Enzymatic activity of the C-subunit is thus regulated by cAMP via the R-subunit, which toggles between cAMP and C-subunit bound states. The R-subunit is composed of a dimerization/docking (D/D) domain connected to two cAMP-binding domains (cAMP:A and cAMP:B). While crystal structures of the free C-subunit and cAMP-bound states of a deletion mutant of the R-subunit are known, there is no structure of the holoenzyme complex or of the cAMP-free state of the R-subunit. An important step in understanding the cAMP-dependent activation of PKA is to map the R-C interface and characterize the mutually exclusive interactions of the R-subunit with cAMP and C-subunit. Amide hydrogen/deuterium exchange mass spectrometry is a suitable method that has provided insights into the different states of the R-subunit in solution, thereby allowing mapping of the effects of cAMP and C-subunit on different regions of the R-subunit. Our study has localized interactions with the C-subunit to a small contiguous surface on the cAMP:A domain and the linker region. In addition, C-subunit binding causes increased amide hydrogen exchange within both cAMP-domains, suggesting that these regions become more flexible in the holoenzyme and are primed to bind cAMP. Furthermore, the difference in the protection patterns between RIalpha and the previously studied RIIbeta upon cAMP-binding suggests isoform-specific differences in cAMP-dependent regulation of PKA activity.     

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.     

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 alpha1 and alpha2 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.     

A point mutation abolishes binding of cAMP to site A in the regulatory subunit of cAMP-dependent protein kinase

Each regulatory subunit of cAMP-dependent protein kinase has two tandem cAMP-binding sites, A and B, at the carboxyl terminus. Based on sequence homologies with the cAMP-binding domain of the Escherichia coli catabolite gene activator protein, a model has been constructed for each cAMP-binding domain. Two of the conserved features of each cAMP-binding site are an arginine and a glutamic acid which interact with the negatively charged phosphate and with the 2'-OH on the ribose ring, respectively. In the type I regulatory subunit, this arginine in cAMP binding site A is Arg-209. Recombinant DNA techniques have been used to change this arginine to a lysine. The resulting protein binds cAMP with a high affinity and associates with the catalytic subunit to form holoenzyme. The mutant holoenzyme also is activated by cAMP. However, the mutant R-subunit binds only 1 mol of cAMP/R-monomer. Photoaffinity labeling confirmed that the mutant R-subunit has only one functional cAMP-binding site. In contrast to the native R-subunit which is labeled at Trp-260 and Tyr-371 by 8-N3cAMP, the mutant R-subunit is convalently modified at a single site, Tyr-371, which correlates with a functional cAMP-binding site B. The lack of functional cAMP-binding site A also was confirmed by activating the mutant holoenzyme with analogs of cAMP which have a high specificity for either site A or site B. 8-NH2-methyl cAMP which preferentially binds to site B was similar to cAMP in its ability to activate both mutant and wild type holoenzyme whereas N6-monobutyryl cAMP, a site A-specific analog, was a very poor activator of the mutant holoenzyme. The results support the conclusions that 1) Arg-209 is essential for cAMP binding to site A and 2) cAMP binding to domain A is not essential for dissociation of the mutant holoenzyme.     

Ligand-dependent perturbation of the conformational ensemble for the GPCR β2 adrenergic receptor revealed by HDX

Mechanism of G protein-coupled receptor (GPCR) activation and their modulation by functionally distinct ligands remains elusive. Using the technique of amide hydrogen/deuterium exchange coupled with mass spectrometry, we examined the ligand-induced changes in conformational states and stability within the beta-2-adrenergic receptor (β(2)AR). Differential HDX reveals ligand-specific alterations in the energy landscape of the receptor's conformational ensemble. The inverse agonists timolol and carazolol were found to be most stabilizing even compared with the antagonist alprenolol, notably in intracellular regions where G proteins are proposed to bind, while the agonist isoproterenol induced the largest degree of conformational mobility. The partial agonist clenbuterol displayed conformational effects found in both the inverse agonists and the agonist. This study highlights the?regional plasticity of the receptor and characterizes unique conformations spanning the entire receptor sequence stabilized by functionally selective ligands, all of which differ from the profile for the apo receptor.     

Solution X-ray scattering evidence for agonist- and antagonist-induced modulation of cleft closure in a glutamate receptor ligand-binding domain

Agonist-induced conformational changes in the ligand-binding domains (LBD) of glutamate receptor ion channels provide the driving force for molecular rearrangements that mediate channel opening and subsequent desensitization. The resulting regulated transmembrane ion fluxes form the basis for most excitatory neuronal signaling in the brain. Crystallographic analysis of the GluR2 LBD core has revealed a ligand-binding cleft located between two lobes. Channel antagonists stabilize an open cleft, whereas agonists stabilize a closed cleft. The crystal structure of the apo form is similar to the antagonist-bound, open state. To understand the conformational behavior of the LBD in the absence of crystal lattice constraints, and thus better to appreciate the thermodynamic constraints on ligand binding, we have undertaken a solution x-ray scattering study using two different constructs encoding either the core or an extended LBD. In agreement with the GluR2 crystal structures, the LBD is more compact in the presence of agonist than it is in the presence of antagonist. However, the time-averaged conformation of the ligand-free core in solution is intermediate between the open, antagonist-bound state and the closed, agonist-bound state, suggesting a conformational equilibrium. Addition of peptide moieties that connect the core domain to the other functional domains in each channel subunit appears to constrain the conformational equilibrium in favor of the open state.     

Conformational transitions in protein-protein association: binding of fasciculin-2 to acetylcholinesterase

The neurotoxin fasciculin-2 (FAS2) is a picomolar inhibitor of synaptic acetylcholinesterase (AChE). The dynamics of binding between FAS2 and AChE is influenced by conformational fluctuations both before and after protein encounter. Submicrosecond molecular dynamics trajectories of apo forms of fasciculin, corresponding to different conformational substates, are reported here with reference to the conformational changes of loop I of this three-fingered toxin. This highly flexible loop exhibits an ensemble of conformations within each substate corresponding to its functions. The high energy barrier found between the two major substates leads to transitions that are slow on the timescale of the diffusional encounter of noninteracting FAS2 and AChE. The more stable of the two apo substates may not be the one observed in the complex with AChE. It seems likely that the more stable apo form binds rapidly to AChE and conformational readjustments then occur in the resulting encounter complex.     

Mapping the Free Energy Landscape of PKA Inhibition and Activation: A Double-Conformational Selection Model for the Tandem cAMP-Binding Domains of PKA RIα

Protein Kinase A (PKA) is the major receptor for the cyclic adenosine monophosphate (cAMP) secondary messenger in eukaryotes. cAMP binds to two tandem cAMP-binding domains (CBD-A and -B) within the regulatory subunit of PKA (R), unleashing the activity of the catalytic subunit (C). While CBD-A in RIα is required for PKA inhibition and activation, CBD-B functions as a “gatekeeper” domain that modulates the control exerted by CBD-A. Preliminary evidence suggests that CBD-B dynamics are critical for its gatekeeper function. To test this hypothesis, here we investigate by Nuclear Magnetic Resonance (NMR) the two-domain construct RIα (91–379) in its apo, cAMP₂, and C-bound forms. Our comparative NMR analyses lead to a double conformational selection model in which each apo CBD dynamically samples both active and inactive states independently of the adjacent CBD within a nearly degenerate free energy landscape. Such degeneracy is critical to explain the sensitivity of CBD-B to weak interactions with C and its high affinity for cAMP. Binding of cAMP eliminates this degeneracy, as it selectively stabilizes the active conformation within each CBD and inter-CBD contacts, which require both cAMP and W260. The latter is contributed by CBD-B and mediates capping of the cAMP bound to CBD-A. The inter-CBD interface is dispensable for intra-CBD conformational selection, but is indispensable for full activation of PKA as it occludes C-subunit recognition sites within CBD-A. In addition, the two structurally homologous cAMP-bound CBDs exhibit marked differences in their residual dynamics profiles, supporting the notion that conservation of structure does not necessarily imply conservation of dynamics.     

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.     

Testing a flexible-receptor docking algorithm in a model binding site

Sampling receptor flexibility is challenging for database docking. We consider a method that treats multiple flexible regions of the binding site independently, recombining them to generate different discrete conformations. This algorithm scales linearly rather than exponentially with the receptor's degrees of freedom. The method was first evaluated for its ability to identify known ligands of a hydrophobic cavity mutant of T4 lysozyme (L99A). Some 200000 molecules of the Available Chemical Directory (ACD) were docked against an ensemble of cavity conformations. Surprisingly, the enrichment of known ligands from among a much larger number of decoys in the ACD was worse than simply docking to the apo conformation alone. Large decoys, accommodated in the larger cavity conformations sampled in the ensemble, were ranked better than known small ligands. The calculation was redone with an energy correction term that considered the cost of forming the larger cavity conformations. Enrichment improved, as did the balance between high-ranking large and small ligands. In a second retrospective test, the ACD was docked against a conformational ensemble of thymidylate synthase. Compared to docking against individual enzyme conformations, the flexible receptor docking approach improved enrichment of known ligands. Including a receptor conformational energy weighting term improved enrichment further. To test the method prospectively, the ACD database was docked against another cavity mutant of lysozyme (L99A/M102Q). A total of 18 new compounds predicted to bind this polar cavity and to change its conformation were tested experimentally; 14 were found to bind. The bound structures for seven ligands were determined by X-ray crystallography. The predicted geometries of these ligands all corresponded to the observed geometries to within 0.7A RMSD or better. Significant conformational changes of the cavity were observed in all seven complexes. In five structures, part of the observed accommodations were correctly predicted; in two structures, the receptor conformational changes were unanticipated and thus never sampled. These results suggest that although sampling receptor flexibility can lead to novel ligands that would have been missed when docking a rigid structure, it is also important to consider receptor conformational energy.     

Solution structure and fluctuation of the Mg(2+)-bound form of calmodulin C-terminal domain

Calmodulin (CaM) is a Ca(2+)-binding protein that functions as a ubiquitous Ca(2+)-signaling molecule, through conformational changes from the "closed" apo conformation to the "open" Ca(2+)-bound conformation. Mg(2+) also binds to CaM and stabilizes its folded structure, but the NMR signals are broadened by slow conformational fluctuations. Using the E104D/E140D mutant, designed to decrease the signal broadening in the presence of Mg(2+) with minimal perturbations of the overall structure, the solution structure of the Mg(2+)-bound form of the CaM C-terminal domain was determined by multidimensional NMR spectroscopy. The Mg(2+)-induced conformational change mainly occurred in EF hand IV, while EF-hand III retained the apo structure. The helix G and helix H sides of the binding sequence undergo conformational changes needed for the Mg(2+) coordination, and thus the helices tilt slightly. The aromatic rings on helix H move to form a new cluster of aromatic rings in the hydrophobic core. Although helix G tilts slightly to the open orientation, the closed conformation is maintained. The fact that the Mg(2+)-induced conformational changes in EF-hand IV and the hydrophobic core are also seen upon Ca(2+) binding suggests that the Ca(2+)-induced conformational changes can be divided into two categories, those specific to Ca(2+) and those common to Ca(2+) and Mg(2+).     

Structure of yeast regulatory subunit: a glimpse into the evolution of PKA signaling

The major cAMP receptors in eukaryotes are the regulatory (R) subunits of PKA, an allosteric enzyme conserved in fungi through mammals. While mammals have four R-subunit genes, Saccharomyces cerevisiae has only one, Bcy1. To achieve a molecular understanding of PKA activation in yeast and to explore the evolution of cyclic-nucleotide binding (CNB) domains, we solved the structure of cAMP-bound Bcy1(168-416). Surprisingly, the relative orientation of the two CNB domains in Bcy1 is very different from mammalian R-subunits. This quaternary structure is defined primarily by a fungi-specific sequence in the hinge between the αB/αC helices of the CNB-A domain. The unique interface between the two CNB domains in Bcy1 defines the allosteric mechanism for cooperative activation of PKA by cAMP. Some interface motifs are isoform-specific while others, although conserved, play surprisingly different roles in each R-subunit. Phylogenetic analysis shows that structural differences in Bcy1 are shared by fungi of the subphylum Saccharomycotina.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Influence of orthosteric ligand binding on the conformational dynamics of the β-2-adrenergic receptor via essential dynamics sampling simulation

The influence of the binding of orthosteric ligands on the conformational dynamics of the β-2-adrenoreceptor was identified using the molecular dynamics method. It was found that there was a small fraction of active states of the receptor in its apo (ligand free) ensemble. An analysis of the MD trajectories indicated that this spontaneous activation of the receptor was accompanied by the motion of its VI helix. Thus, the receptor’s constitutive activity is a direct result of its conformational dynamics. On the other hand, the binding of the full agonist resulted in a significant shift in the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized receptor in its inactive state. It is likely that the binding of the inverse agonists might be a universal method of the constitutive activity inhibition. Our results indicate that ligand binding redistributes preexisting conformational degrees of freedom (in accordance to the Monod-Wyman-Changeux Model), rather than causes an induced fit. Therefore, the ensemble of the biologically relevant receptors conformations has been encoded in its spatial structure and individual conformations from that ensemble and might be used by the cell according to the physiological behavior.     

Opening and closing motions in the periplasmic vitamin B12 binding protein BtuF

BtuF is the periplasmic binding protein (PBP) in the vitamin B(12) uptake system in Escherichia coli where it is associated with the ABC transporter BtuCD. When the ligand binds, PBPs generally display large conformational changes, commonly termed the Venus flytrap mechanism. BtuF belongs to a group of PBPs that, on the basis of crystal structures, does not appear to display such behavior. Using 480 ns multicopy molecular dynamics simulations of apo and holo forms of the protein, we investigate the dynamics of BtuF. We find BtuF to be more flexible than previously assumed, displaying clear opening and closing motions which are more pronounced in the apo form. The protein behavior is compatible with a PBP functional model that postulates a closed conformation for the ligand-bound state, whereas the empty form fluctuates between open and closed conformations. Elastic network normal-mode analysis suggests that all BtuF-like PBPs are capable of similar opening and closing motions. It also makes the typical Venus flytrap domain motions a likely common means of how PBP-ABC transporter interaction could occur.     

Structural Analyses of a Constitutively Active Mutant of Exchange Protein Directly Activated by cAMP

Exchange proteins directly activated by cAMP (EPACs) are important allosteric regulators of cAMP-mediated signal transduction pathways. To understand the molecular mechanism of EPAC activation, we have combined site-directed mutagenesis, X-ray crystallography, and peptide amide hydrogen/deuterium exchange mass spectrometry (DXMS) to probe the structural and conformational dynamics of EPAC2-F435G, a constitutively active EPAC2 mutant. Our study demonstrates that conformational dynamics plays a critical role in cAMP-induced EPAC activation. A glycine mutation at 435 position shifts the equilibrium of conformational dynamics towards the extended active conformation.