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.     

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.     

Normal-mode-analysis-guided investigation of crucial intersubunit contacts in the cAMP-dependent gating in HCN channels

Protein structures define a complex network of atomic interactions in three dimensions. Direct visualization of the structure and analysis of the interaction potential energy are not straightforward approaches to pinpoint the atomic contacts that are crucial for protein function. We used the tetrameric hyperpolarization-activated cAMP-regulated (HCN) channel as a model system to study the intersubunit contacts in cAMP-dependent gating. To obtain a systematic survey of the contacts between each pair of residues, we used normal-mode analysis, a computational approach for studying protein dynamics, and constructed the covariance matrix for C-α atoms. The significant contacts revealed by covariance analysis were further investigated by means of mutagenesis and functional assays. Among the mutant channels that show phenotypes different from those of the wild-type, we focused on two mutant channels that express opposite changes in cAMP-dependent gating. Subsequent biochemical assays on isolated C-terminal fragments, including the cAMP binding domain, revealed only minimal effects on cAMP binding, suggesting the necessity of interpreting the cAMP-dependent allosteric regulation at the whole-channel level. For this purpose, we applied the patch-clamp fluorometry technique and observed correlated changes in the dynamic, state-dependent cAMP binding in the mutant channels. This study not only provides further understanding of the intersubunit contacts in allosteric coupling in the HCN channel, it also illustrates an effective strategy for delineating important atomic contacts within a structure.     

cAMP Modulation of the cytoplasmic domain in the HCN2 channel investigated by molecular simulations

The hyperpolarization-activated cyclic nucleotide-modulated (HCN) cation channels are opened by membrane hyperpolarization, while their activation is modulated by the binding of cyclic adenosine monophosphate (cAMP) in the cytoplasm. Here we investigate the molecular basis of cAMP channel modulation by performing molecular dynamics simulations of a segment comprising the C-linker and the cyclic nucleotide binding domain (CNBD) in the presence and absence of cAMP, based on the available crystal structure of HCN2 from mouse. In presence of cAMP, the protein undergoes an oscillation of the quaternary structure on the order of 10 ns, not observed in the apoprotein. In contrast, the absence of ligand causes conformational rearrangements within the CNBDs, driving these domains to a more flexible state, similar to that described in CNBDs of other proteins. This increased flexibility causes a rather disordered movement of the CNBDs, resulting in an inhibitory effect on the channel. We propose that the cAMP-triggered large-scale oscillation plays an important role for the channel's function, being coupled to a motion of the C-linker which, in turn, modulates the gating of the channel.     

Structural Basis for the cAMP-dependent Gating in the Human HCN4 Channel

Hyperpolarization-activated cAMP-regulated (HCN) channels play important physiological roles in both cardiovascular and central nervous systems. Among the four HCN isoforms, HCN2 and HCN4 show high expression levels in the human heart, with HCN4 being the major cardiac isoform. The previously published crystal structure of the mouse HCN2 (mHCN2) C-terminal fragment, including the C-linker and the cyclic-nucleotide binding domain (CNBD), has provided many insights into cAMP-dependent gating in HCN channels. However, structures of other mammalian HCN channel isoforms have been lacking. Here we used a combination of approaches including structural biology, biochemistry, and electrophysiology to study cAMP-dependent gating in HCN4 channel. First we solved the crystal structure of the C-terminal fragment of human HCN4 (hHCN4) channel at 2.4 Å. Overall we observed a high similarity between mHCN2 and hHCN4 crystal structures. Functional comparison between two isoforms revealed that compared with mHCN2, the hHCN4 protein exhibited marked different contributions to channel function, such as a ∼3-fold reduction in the response to cAMP. Guided by structural differences in the loop region between β4 and β5 strands, we identified residues that could partially account for the differences in response to cAMP between mHCN2 and hHCN4 proteins. Moreover, upon cAMP binding, the hHCN4 C-terminal protein exerts a much prolonged effect in channel deactivation that could have significant physiological contributions.     

State-dependent cAMP binding to functioning HCN channels studied by patch-clamp fluorometry

One major goal of ion channel research is to delineate the molecular events from the detection of the stimuli to the movement of channel gates. For ligand-gated channels, it is challenging to separate ligand binding from channel gating. Here we studied the cyclic adenosine monophosphate (cAMP)-dependent gating in hyperpolarization-activated cAMP-regulated (HCN) channel by simultaneously recording channel opening and ligand binding, using the patch-clamp fluorometry technique with a unique fluorescent cAMP analog that fluoresces strongly in the hydrophobic binding pocket and exerts regulatory effects on HCN channels similar to those imposed by cAMP. Corresponding to voltage-dependent channel activation, we observed a robust, close-to-threefold increase in ligand binding, which was more pronounced at subsaturating ligand concentrations than higher concentrations. This observation supported the cyclic allosteric models and indicated that protein allostery can be implemented through differentiating ligand binding affinities between resting and active states. The kinetics of ligand binding largely matched channel activation. However, during channel deactivation, ligand unbinding was slower than channel closing, suggesting a delayed response to membrane potential by the ligand binding machinery. Our results provide what we believe to be new insights into the cAMP-dependent gating in HCN channel and the interpretation of protein allostery for general ligand-gated channels and receptors.     

Tetramerization dynamics of C-terminal domain underlies isoform-specific cAMP gating in hyperpolarization-activated cyclic nucleotide-gated channels

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are dually activated by hyperpolarization and binding of cAMP to their cyclic nucleotide binding domain (CNBD). HCN isoforms respond differently to cAMP; binding of cAMP shifts activation of HCN2 and HCN4 by 17 mV but shifts that of HCN1 by only 2-4 mV. To explain the peculiarity of HCN1, we solved the crystal structures and performed a biochemical-biophysical characterization of the C-terminal domain (C-linker plus CNBD) of the three isoforms. Our main finding is that tetramerization of the C-terminal domain of HCN1 occurs at basal cAMP concentrations, whereas those of HCN2 and HCN4 require cAMP saturating levels. Therefore, HCN1 responds less markedly than HCN2 and HCN4 to cAMP increase because its CNBD is already partly tetrameric. This is confirmed by voltage clamp experiments showing that the right-shifted position of V(½) in HCN1 is correlated with its propensity to tetramerize in vitro. These data underscore that ligand-induced CNBD tetramerization removes tonic inhibition from the pore of HCN channels.     

A novel dimerization interface of cyclic nucleotide binding domain, which is disrupted in presence of cAMP: implications for CNG channels gating

Cyclic nucleotide binding domain (CNBD) is a ubiquitous domain of effector proteins involved in signalling cascades of prokaryota and eukaryota. CNBD activation by cyclic nucleotide monophosphate (cNMP) is studied well in the case of several proteins. However, this knowledge is hardly applicable to cNMP-modulated cation channels. Despite the availability of CNBD crystal structures of bacterial cyclic nucleotide-gated (CNG) and mammalian hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels in presence and absence of the cNMP, the full understanding of CNBD conformational changes during activation is lacking. Here, we describe a novel CNBD dimerization interface found in crystal structures of bacterial CNG channel MlotiK1 and mammalian cAMP-activated guanine nucleotide-exchange factor Epac2. Molecular dynamics simulations show that the found interface is stable on the studied timescale of 100?ns, in contrast to the dimerization interface, reported previously. Comparisons with cN-bound structures of CNBD show that the dimerization is incompatible with cAMP binding. Thus, the cAMP-dependent monomerization of CNBD may be an alternative mechanism of the cAMP sensing. Based on these findings, we propose a model of the bacterial CNG channel modulation by cAMP.     

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.     

Regulation of hyperpolarization-activated HCN channel gating and cAMP modulation due to interactions of COOH terminus and core transmembrane regions

Members of the hyperpolarization-activated cation (HCN) channel family generate HCN currents (I(h)) that are directly regulated by cAMP and contribute to pacemaking activity in heart and brain. The four different HCN isoforms show distinct biophysical properties. In cell-free patches from Xenopus oocytes, the steady-state activation curve of HCN2 channels is 20 mV more hyperpolarized compared with HCN1. Whereas the binding of cAMP to a COOH-terminal cyclic nucleotide binding domain (CNBD) markedly shifts the activation curve of HCN2 by 17 mV to more positive potentials, the response of HCN1 is much less pronounced (4 mV shift). A previous deletion mutant study suggested that the CNBD inhibits hyperpolarization-gating in the absence of cAMP; the binding of cAMP shifts gating to more positive voltages by relieving this inhibition. The differences in basal gating and cAMP responsiveness between HCN1 and HCN2 were proposed to result from a greater inhibitory effect of the CNBD in HCN2 compared with HCN1. Here, we use a series of chimeras between HCN1 and HCN2, in which we exchange the NH(2) terminus, the transmembrane domain, or distinct domains of the COOH terminus, to investigate further the molecular bases for the modulatory action of cAMP and for the differences in the functional properties of the two channels. Differences in cAMP regulation between HCN1 and HCN2 are localized to sequence differences within the COOH terminus of the two channels. Surprisingly, exchange of the CNBDs between HCN1 and HCN2 has little effect on basal gating and has only a modest one on cAMP modulation. Rather, differences in cAMP modulation depend on the interaction between the CNBD and the C-linker, a conserved 80-amino acid region that connects the last (S6) transmembrane segment to the CNBD. Differences in basal gating depend on both the core transmembrane domain and the COOH terminus. These data, taken in the context of the previous data on deletion mutants, suggest that the inhibitory effect of the CNBD on basal gating depends on its interactions with both the C-linker and core transmembrane domain of the channel. The extent to which cAMP binding is able to relieve this inhibition is dependent on the interaction between the C-linker and the CNBD.     

Investigating the disorder–order transition of calmodulin binding domain upon binding calmodulin using molecular dynamics simulation

We have studied the conformational transition of the calmodulin binding domains (CBD) in several calmodulin‐binding kinases, in which CBD changes from the disordered state to the ordered state when binding with calmodulin (CaM). Targeted molecular dynamics simulation was used to investigate the binding process of CaM and CBD of CaM‐dependent kinase I (CaMKI–CBD). The results show that CaMKI–CBD began to form an α‐helix and the interaction free energy between CaM and CaMKI–CBD increased once CaM fully encompassed CaMKI–CBD. Two series of CaM/CBD complex systems, including the complexes of CaM with the initially disordered and the final ordered CBD, were constructed to study the interaction using molecular dynamics simulations. Our analyses suggest that the VDW interaction plays a dominant role in CaM/CBD binding and is a key factor in the disorder–order transition of CBD. Additionally, the entropy effect is not in favor of the formation of the CaM/CBD complex, which is consistent with the experimental evidence. Based on the results, it appears that the CBD conformational change from a non‐compact extended structure to compact α‐helix is critical in gaining a favorable VDW interaction and interaction free energy. Copyright © 2009 John Wiley & Sons, Ltd.     

A conserved tripeptide in CNG and HCN channels regulates ligand gating by controlling C-terminal oligomerization

Cyclic nucleotides directly enhance the opening of the tetrameric CNG and HCN channels, although the mechanism remains unclear. We examined why HCN and certain CNG subunits form functional homomeric channels, whereas other CNG subunits only function in heteromeric channels. The "defect" in the CNGA4 subunit that prevents its homomeric expression was localized to its C-linker, which connects the transmembrane domain to the binding domain and contains a tripeptide that decreases the efficacy of ligand gating. Remarkably, replacement of the homologous HCN tripeptide with the CNGA4 sequence transformed cAMP into an inverse agonist that inhibits HCN channel opening. Using analytical ultracentrifugation, we identified the structural basis for this gating switch: whereas cAMP normally enhances the assembly of HCN C-terminal domains into a tetrameric gating ring, inclusion of the CNGA4 tripeptide reversed this action so that cAMP now causes gating ring disassembly. Thus, ligand gating depends on the dynamic oligomerization of C-terminal binding domains.     

Transmembrane signaling of chemotaxis receptor tar: insights from molecular dynamics simulation studies

Transmembrane signaling of chemotaxis receptors has long been studied, but how the conformational change induced by ligand binding is transmitted across the bilayer membrane is still elusive at the molecular level. To tackle this problem, we carried out a total of 600-ns comparative molecular dynamics simulations (including model-building simulations) of the chemotaxis aspartate receptor Tar (a part of the periplasmic domain/transmembrane domain/HAMP domain) in explicit lipid bilayers. These simulations reveal valuable insights into the mechanistic picture of Tar transmembrane signaling. The piston-like movement of a transmembrane helix induced by ligand binding on the periplasmic side is transformed into a combination of both longitudinal and transversal movements of the helix on the cytoplasmic side as a result of different protein-lipid interactions in the ligand-off and ligand-on states of the receptor. This conformational change alters the dynamics and conformation of the HAMP domain, which is presumably a mechanism to deliver the signal from the transmembrane domain to the cytoplasmic domain. The current results are consistent with the previously suggested dynamic bundle model in which the HAMP dynamics change is a key to the signaling. The simulations provide further insights into the conformational changes relevant to the HAMP dynamics changes in atomic detail.     

Trafficking and gating of hyperpolarization-activated cyclic nucleotide-gated channels are regulated by interaction with tetratricopeptide repeat-containing Rab8b-interacting protein (TRIP8b) and cyclic AMP at distinct sites

Ion channel trafficking and gating are often influenced by interactions with auxiliary subunits. Tetratricopeptide repeat-containing Rab8b-interacting protein (TRIP8b) is an auxiliary subunit for neuronal hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. TRIP8b interacts directly with two distinct sites of HCN channel pore-forming subunits to control channel trafficking and gating. Here we use mutagenesis combined with electrophysiological studies to define and distinguish the functional importance of the HCN/TRIP8b interaction sites. Interaction with the last three amino acids of the HCN1 C terminus governed the effect of TRIP8b on channel trafficking, whereas TRIP8b interaction with the HCN1 cyclic nucleotide binding domain (CNBD) affected trafficking and gating. Biochemical studies revealed that direct interaction between TRIP8b and the HCN1 CNBD was disrupted by cAMP and that TRIP8b binding to the CNBD required an arginine residue also necessary for cAMP binding. In accord, increasing cAMP levels in cells antagonized the up-regulation of HCN1 channels mediated by a TRIP8b construct binding the CNBD exclusively. These data illustrate the distinct roles of the two TRIP8b-HCN interaction domains and suggest that TRIP8b and cAMP may directly compete for binding the HCN CNBD to control HCN channel gating, kinetics, and trafficking.     

Flavonoid Regulation of HCN2 Channels

The hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are pacemaker channels whose currents contribute to rhythmic activity in the heart and brain. HCN channels open in response to hyperpolarizing voltages, and the binding of cAMP to their cyclic nucleotide-binding domain (CNBD) facilitates channel opening. Here, we report that, like cAMP, the flavonoid fisetin potentiates HCN2 channel gating. Fisetin sped HCN2 activation and shifted the conductance-voltage relationship to more depolarizing potentials with a half-maximal effective concentration (EC₅₀) of 1.8 μm. When applied together, fisetin and cAMP regulated HCN2 gating in a nonadditive fashion. Fisetin did not potentiate HCN2 channels lacking their CNBD, and two independent fluorescence-based binding assays reported that fisetin bound to the purified CNBD. These data suggest that the CNBD mediates the fisetin potentiation of HCN2 channels. Moreover, binding assays suggest that fisetin and cAMP partially compete for binding to the CNBD. NMR experiments demonstrated that fisetin binds within the cAMP-binding pocket, interacting with some of the same residues as cAMP. Together, these data indicate that fisetin is a partial agonist for HCN2 channels.     

RIalpha subunit of PKA: a cAMP-free structure reveals a hydrophobic capping mechanism for docking cAMP into site B

In eukaryotes the primary target for cAMP, a ubiquitous second messenger, is cAMP-dependent protein kinase (PKA). Understanding how binding and release of cAMP changes the cAMP binding domains and then triggers long-range allosteric responses is an important challenge. This conformational switching requires structure solutions of cAMP binding domains in cAMP-bound and cAMP-free states. We describe for the first time a crystal structure of the cAMP binding domains of PKA type Ialpha regulatory subunit where site A is occupied by cGMP and site B is unoccupied. The structure reveals that the carboxyl terminus of domain B serves as a hydrophobic cap, locking the cyclic nucleotide via its adenine ring into the beta-barrel. In the absence of cAMP, the "cap" is released via an extension of the C-terminal helix. This simple hinge mechanism for binding and release of cAMP also provides a mechanism for allosteric communication between sites A and B.     

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.     

Integrated allosteric model of voltage gating of HCN channels

Hyperpolarization-activated (pacemaker) channels are dually gated by negative voltage and intracellular cAMP. Kinetics of native cardiac f-channels are not compatible with HH gating, and require closed/open multistate models. We verified that members of the HCN channel family (mHCN1, hHCN2, hHCN4) also have properties not complying with HH gating, such as sigmoidal activation and deactivation, activation deviating from fixed power of an exponential, removal of activation "delay" by preconditioning hyperpolarization. Previous work on native channels has indicated that the shifting action of cAMP on the open probability (Po) curve can be accounted for by an allosteric model, whereby cAMP binds more favorably to open than closed channels. We therefore asked whether not only cAMP-dependent, but also voltage-dependent gating of hyperpolarization-activated channels could be explained by an allosteric model. We hypothesized that HCN channels are tetramers and that each subunit comprises a voltage sensor moving between "reluctant" and "willing" states, whereas voltage sensors are independently gated by voltage, channel closed/open transitions occur allosterically. These hypotheses led to a multistate scheme comprising five open and five closed channel states. We estimated model rate constants by fitting first activation delay curves and single exponential time constant curves, and then individual activation/deactivation traces. By simply using different sets of rate constants, the model accounts for qualitative and quantitative aspects of voltage gating of all three HCN isoforms investigated, and allows an interpretation of the different kinetic properties of different isoforms. For example, faster kinetics of HCN1 relative to HCN2/HCN4 are attributable to higher HCN1 voltage sensors' rates and looser voltage-independent interactions between subunits in closed/open transitions. It also accounts for experimental evidence that reduction of sensors' positive charge leads to negative voltage shifts of Po curve, with little change of curve slope. HCN voltage gating thus involves two processes: voltage sensor gating and allosteric opening/closing.     

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.     

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.