Predicting Allosteric Effects from Orthosteric Binding in Hsp90-Ligand Interactions: Implications for Fragment-Based Drug Design

A key question in mapping dynamics of protein-ligand interactions is to distinguish changes at binding sites from those associated with long range conformational changes upon binding at distal sites. This assumes a greater challenge when considering the interactions of low affinity ligands (dissociation constants, K_D, in the μM range or lower). Amide hydrogen deuterium Exchange mass spectrometry (HDXMS) is a robust method that can provide both structural insights and dynamics information on both high affinity and transient protein-ligand interactions. In this study, an application of HDXMS for probing the dynamics of low affinity ligands to proteins is described using the N-terminal ATPase domain of Hsp90. Comparison of Hsp90 dynamics between high affinity natural inhibitors (K_D ~ nM) and fragment compounds reveal that HDXMS is highly sensitive in mapping the interactions of both high and low affinity ligands. HDXMS reports on changes that reflect both orthosteric effects and allosteric changes accompanying binding. Orthosteric sites can be identified by overlaying HDXMS onto structural information of protein-ligand complexes. Regions distal to orthosteric sites indicate long range conformational changes with implications for allostery. HDXMS, thus finds powerful utility as a high throughput method for compound library screening to identify binding sites and describe allostery with important implications for fragment-based ligand discovery (FBLD).     

Dynamics of phosphodiesterase-induced cAMP dissociation from protein kinase A: Capturing transient ternary complexes by HDXMS

cAMP signaling is a fundamental cellular process necessary for mediating responses to hormonal stimuli. In contrast to cAMP-dependent activation of protein kinase A (PKA), an important cellular target, far less is known on termination in cAMP signaling, specifically how phosphodiesterases (PDEs) facilitate dissociation and hydrolysis of bound cAMP. In this study, we have probed the dynamics of a ternary complex of PKA and a PDE–RegA with an excess of a PDE-nonhydrolyzable cAMP analog, Sp-cAMPS by amide hydrogen/deuterium exchange mass spectrometry (HDXMS). Our results highlight how HDXMS can be used to monitor reactions together with mapping conformational dynamics of transient signaling complexes. Our results confirm a two-state model for active RegA-mediated dissociation of bound cAMP. Further, our results reveal that Sp-cAMPS and RegA mediate mutually exclusive interactions with the same region of PKA and at specific concentrations of Sp-cAMPS, RegA is capable of blocking Sp-cAMPS reassociation to PKA. This provides a molecular basis for how PDEs modulate levels of intracellular cAMP so that PKA is better suited to responding to fluxes rather than constant levels of cAMP. This study underscores how HDXMS can be a powerful tool for monitoring reactions together with mapping conformational dynamics in signaling proteins. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

Insights into enzyme structure and dynamics elucidated by amide H/D exchange mass spectrometry

Conformational changes and protein dynamics play an important role in the catalytic efficiency of enzymes. Amide H/D exchange mass spectrometry (H/D exchange MS) is emerging as an efficient technique to study the local and global changes in protein structure and dynamics due to ligand binding, protein activation-inactivation by modification, and protein-protein interactions. By monitoring the selective exchange of hydrogen for deuterium along a peptide backbone, this sensitive technique probes protein motions and structural elements that may be relevant to allostery and function. In this report, several applications of H/D exchange MS are presented which demonstrate the unique capability of amide hydrogen/deuterium exchange mass spectrometry for examining dynamic and structural changes associated with enzyme catalysis.     

Conserved core of amyloid fibrils of wild type and A30P mutant α-synuclein

The major component of neural inclusions that are the pathological hallmark of Parkinson's disease are amyloid fibrils of the protein α-synuclein (aS). Here we investigated if the disease-related mutation A30P not only modulates the kinetics of aS aggregation, but also alters the structure of amyloid fibrils. To this end we optimized the method of quenched hydrogen/deuterium exchange coupled to NMR spectroscopy and performed two-dimensional proton-detected high-resolution magic angle spinning experiments. The combined data indicate that the A30P mutation does not cause changes in the number, location and overall arrangement of β-strands in amyloid fibrils of aS. At the same time, several residues within the fibrillar core retain nano-second dynamics. We conclude that the increased pathogenicity related to the familial A30P mutation is unlikely to be caused by a mutation-induced change in the conformation of aS aggregates.     

Allosteric signaling in the biotin repressor occurs via local folding coupled to global dampening of protein dynamics

The biotin repressor is an allosterically regulated, site-specific DNA-binding protein. Binding of the small ligand bio-5′-AMP activates repressor dimerization, which is a prerequisite to DNA binding. Multiple disorder-to-order transitions, some of which are known to be important for the functional allosteric response, occur in the vicinity of the ligand-binding site concomitant with effector binding to the repressor monomer. In this work, the extent to which these local changes are coupled to additional changes in the structure/dynamics of the repressor was investigated using hydrogen/deuterium exchange coupled to mass spectrometry. Measurements were performed on the apo-protein and on complexes of the protein bound to four different effectors that elicit a range of thermodynamic responses in the repressor. Global exchange measurements indicate that binding of any effector to the intact protein is accompanied by protection from exchange. Mass spectrometric analysis of pepsin-cleavage products generated from the exchanged complexes reveals that the protection is distributed throughout the protein. Furthermore, the magnitude of the level of protection in each peptide from hydrogen/deuterium exchange correlates with the magnitude of the functional allosteric response elicited by a ligand. These results indicate that local structural changes in the binding site that occur concomitant with effector binding nucleate global dampening of dynamics. Moreover, the magnitude of dampening of repressor dynamics tracks with the magnitude of the functional response to effector binding.     

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.     

The Deuterator: software for the determination of backbone amide deuterium levels from H/D exchange MS data

Background  

The combination of mass spectrometry and solution phase amide hydrogen/deuterium exchange (H/D exchange) experiments is an effective method for characterizing protein dynamics, and protein-protein or protein-ligand interactions. Despite methodological advancements and improvements in instrumentation and automation, data analysis and display remains a tedious process. The factors that contribute to this bottleneck are the large number of data points produced in a typical experiment, each requiring manual curation and validation, and then calculation of the level of backbone amide exchange. Tools have become available that address some of these issues, but lack sufficient integration, functionality, and accessibility required to address the needs of the H/D exchange community. To date there is no software for the analysis of H/D exchange data that comprehensively addresses these issues.     

Hydrogen exchange shows peptide binding stabilizes motions in Hck SH2.

Src-homology-2 domains are small, 100 amino acid protein modules that are present in a number of signal transduction proteins. Previous NMR studies of SH2 domain dynamics indicate that peptide binding decreases protein motions in the pico- to nanosecond, and perhaps slower, time range. We suggest that amide hydrogen exchange and mass spectrometry may be useful for detecting changes in protein dynamics because hydrogen exchange rates are relatively insensitive to the time domains of the dynamics. In the present study, hydrogen exchange and mass spectrometry were used to probe hematopoietic cell kinase SH2 that was either free or bound to a 12-residue high-affinity peptide. Hydrogen exchange rates were determined by exposing free and bound SH2 to D(2)O, fragmenting the SH2 with pepsin, and determining the deuterium level in the peptic fragments. Binding generally decreased hydrogen exchange along much of the SH2 backbone, indicating a widespread reduction in dynamics. Alterations in the exchange of the most rapidly exchanging amide hydrogens, which was detected following acid quench and analysis by mass spectrometry, were used to locate differences in low-amplitude motion when SH2 was bound to the peptide. In addition, the results indicate that hydrogen exchange from the folded form of SH2 is an important process along the entire SH2 backbone.     

Dimer dissociation is a key energetic event in the fold-switch pathway of KaiB

Cyanobacteria possesses the simplest circadian clock, composed of three proteins that act as a phosphorylation oscillator: KaiA, KaiB, and KaiC. The timing of this oscillator is determined by the fold-switch of KaiB, a structural rearrangement of its C-terminal half that is accompanied by a change in the oligomerization state. During the day, KaiB forms a stable tetramer (gsKaiB), whereas it adopts a monomeric thioredoxin-like fold during the night (fsKaiB). Although the structures and functions of both native states are well studied, little is known about the sequence and structure determinants that control their structural interconversion. Here, we used confinement molecular dynamics (CCR-MD) and folding simulations using structure-based models to show that the dissociation of the gsKaiB dimer is a key energetic event for the fold-switch. Hydrogen-deuterium exchange mass spectrometry (HDXMS) recapitulates the local stability of protein regions reported by CCR-MD, with both approaches consistently indicating that the energy and backbone flexibility changes are solely associated with the region that fold-switches between gsKaiB and fsKaiB and that the localized regions that differentially stabilize gsKaiB also involve regions outside the dimer interface. Moreover, two mutants (R23C and R75C) previously reported to be relevant for altering the rhythmicity of the Kai clock were also studied by HDXMS. Particularly, R75C populates dimeric and monomeric states with a deuterium incorporation profile comparable to the one observed for fsKaiB, emphasizing the importance of the oligomerization state of KaiB for the fold-switch. These findings suggest that the information necessary to control the rhythmicity of the cyanobacterial biological clock is, to a great extent, encoded within the KaiB sequence.     

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.     

HDX‐MS reveals orthosteric and allosteric changes in apolipoprotein‐D structural dynamics upon binding of progesterone

Apolipoprotein‐D is a glycosylated tetrameric lipocalin that binds and transports small hydrophobic molecules such as progesterone and arachidonic acid. Like other lipocalins, apolipoprotein‐D adopts an eight‐stranded β‐barrel fold stabilized by two intramolecular disulphide bonds, with an adjacent α‐helix. Crystallography studies of recombinant apolipoprotein‐D demonstrated no major conformational changes upon progesterone binding. Amide hydrogen‐deuterium exchange mass spectrometry (HDX‐MS) reports structural changes of proteins in solution by monitoring exchange of amide hydrogens in the protein backbone with deuterium. HDX‐MS detects changes in conformation and structural dynamics in response to protein function such as ligand binding that may go undetected in X‐ray crystallography, making HDX‐MS an invaluable orthogonal technique. Here, we report an HDX‐MS protocol for apolipoprotein‐D that solved challenges of high protein rigidity and low pepsin cleavage using rigorous quenching conditions and longer deuteration times, yielding 85% sequence coverage and 50% deuterium exchange. The relative fractional deuterium exchange of ligand‐free apolipoprotein‐D revealed apolipoprotein‐D to be a highly structured protein. Progesterone binding was detected by significant reduction in deuterium exchange in eight peptides. Stabilization of apolipoprotein‐D dynamics can be interpreted as a combined orthosteric effect in the ligand binding pocket and allosteric effect at the N‐terminus and C‐terminus. Together, our experiments provide insight into apolipoprotein‐D structural dynamics and map the effects of progesterone binding that are relayed to distal parts of the protein. The observed stabilization of apolipoprotein‐D dynamics upon progesterone binding demonstrates a common behaviour in the lipocalin family and may have implications for interactions of apolipoprotein‐D with receptors or lipoprotein particles. Statement: We reveal for the first time how apolipoprotein‐D, which is protective in Alzheimer's disease, becomes more ordered when bound to a molecule of steroid hormone. These results significantly extend the understanding of apolipoprotein‐D structure from X‐ray crystallography studies by incorporating information on how protein motion changes over time. To achieve these results an improved protocol was developed, suitable for proteins similar to apolipoprotein‐D, to elucidate how proteins change flexibility when binding to small molecules.     

Solution structure and dynamics of a heat shock protein assembly probed by hydrogen exchange and mass spectrometry

The solution conformation and dynamics of the 16.9 kDa small heat shock protein from wheat have been studied using a combination of hydrogen/deuterium exchange, proteolytic digestion, and mass spectrometry. At room temperature, HSP16.9 exists as a dodecameric assembly. Regions of HSP16.9 that form extensive and essential intersubunit contacts in the assembly, including residues 1-40 and 131-151, show little or no protection against hydrogen/deuterium exchange after incubation in D(2)O for 5 s. The high levels of hydrogen/deuterium exchange indicate that these regions have experienced large conformational fluctuations in solution, breaking intersubunit contacts and exposing buried amide hydrogens to solvent. When HSP16.9 is pulse labeled for 10 ms, residues 1-40 and 131-151 are substantially more protected than they are after 5 s. Thus, the breaking of intersubunit contacts occurs on a time scale between 10 milliseconds and 5 s. At 42 degrees C, HSP16.9 exists in a suboligomeric form. When the intrinsic temperature dependence of hydrogen/deuterium exchange is taken into account, exchange patterns at 25 and 42 degrees C are identical within experimental error, suggesting that the conformation of individual HSP16.9 subunits is the same in both the dodecameric and subdodecameric forms. Significant protection is seen in regions that form the dimeric interface, suggesting that the stable suboligomeric form is a dimer. Taken together, these results suggest that heat activation of HSP16.9 occurs by shifting the dodecamer <--> dimer equilibrium in favor of free dimers. The conformation of the dimers themselves does not appear to be altered with an increase in temperature.     

Relationship between Protein Stabilization and Protein Rigidification Induced by Mannosylglycerate

Understanding protein stabilization by small organic compounds is a topic of great practical importance. The effect of mannosylglycerate, a charged compatible solute typical of thermophilic microorganisms, on a variant of staphylococcal nuclease was investigated using several NMR spectroscopy methods. No structural changes were apparent from the chemical shifts of amide protons. Measurements of ¹⁵N relaxation and model-free analysis, water-amide saturation transfer (phase-modulated CLEAN chemical exchange), and hydrogen/deuterium exchange rates provided a detailed picture of the effects of mannosylglycerate on the backbone dynamics and time-averaged structure of this protein. The widest movements of the protein backbone were significantly constrained in the presence of mannosylglycerate, as indicated by the average 5-fold decrease of the hydrogen/deuterium exchange rates, but the effect on the millisecond timescale was small. At high frequencies, internal motions of staphylococcal nuclease were progressively restricted with increasing concentrations of mannosylglycerate or reduced temperature, while the opposite effect was observed with urea (a destabilizing solute). The order parameters showed a strong correlation with the changes in the T_m values induced by different solutes, determined by differential scanning calorimetry. These data show that mannosylglycerate caused a generalised reduction of backbone motions and demonstrate a correlation between protein stabilization and protein rigidification.     

Protein flexibility directs DNA recognition by the papillomavirus E2 proteins

Although DNA flexibility is known to play an important role in DNA-protein interactions, the importance of protein flexibility is less well understood. Here, we show that protein dynamics are important in DNA recognition using the well-characterized human papillomavirus (HPV) type 6 E2 protein as a model system. We have compared the DNA binding properties of the HPV 6 E2 DNA binding domain (DBD) and a mutant lacking two C-terminal leucine residues that form part of the hydrophobic core of the protein. Deletion of these residues results in increased specific and non-specific DNA binding and an overall decrease in DNA binding specificity. Using (15)N NMR relaxation and hydrogen/deuterium exchange, we demonstrate that the mutation results in increased flexibility within the hydrophobic core and loop regions that orient the DNA binding helices. Stopped-flow kinetic studies indicate that increased flexibility alters DNA binding by increasing initial interactions with DNA but has little or no effect on the structural rearrangements that follow this step. Taken together these data demonstrate that subtle changes in protein dynamics have a major influence on protein-DNA interactions.     

Dissecting interdomain communication within cAPK regulatory subunit type IIbeta using enhanced amide hydrogen/deuterium exchange mass spectrometry (DXMS)

cAMP-dependent protein kinase (cAPK) is a heterotetramer containing a regulatory (R) subunit dimer bound to two catalytic (C) subunits and is involved in numerous cell signaling pathways. The C-subunit is activated allosterically when two cAMP molecules bind sequentially to the cAMP-binding domains, designated A and B (cAB-A and cAB-B, respectively). Each cAMP-binding domain contains a conserved Arg residue that is critical for high-affinity cAMP binding. Replacement of this Arg with Lys affects cAMP affinity, the structural integrity of the cAMP-binding domains, and cAPK activation. To better understand the local and long-range effects that the Arg-to-Lys mutation has on the dynamic properties of the R-subunit, the amide hydrogen/deuterium exchange in the RIIbeta subunit was probed by electrospray mass spectrometry. Mutant proteins containing the Arg-to-Lys substitution in either cAMP-binding domain were deuterated for various times and then, prior to mass spectrometry analysis, subjected to pepsin digestion to localize the deuterium incorporation. Mutation of this Arg in cAB-A (Arg230) causes an increase in amide hydrogen exchange throughout the mutated domain that is beyond the modest and localized effects of cAMP removal and is indicative of the importance of this Arg in domain organization. Mutation of Arg359 (cAB-B) leads to increased exchange in the adjacent cAB-A domain, particularly in the cAB-A domain C-helix that lies on top of the cAB-B domain and is believed to be functionally linked to the cAB-B domain. This interdomain communication appears to be a unidirectional pathway, as mutation of Arg230 in cAB-A does not effect dynamics of the cAB-B domain.     

Conformational analysis of HAMLET, the folding variant of human alpha-lactalbumin associated with apoptosis

A combination of hydrogen/deuterium (H/D) exchange and limited proteolysis experiments coupled to mass spectrometry analysis was used to depict the conformation in solution of HAMLET, the folding variant of human alpha-lactalbumin, complexed to oleic acid, that induces apoptosis in tumor and immature cells. Although near- and far-UV CD and fluorescence spectroscopy were not able to discriminate between HAMLET and apo-alpha-lactalbumin, H/D exchange experiments clearly showed that they correspond to two distinct conformational states, with HAMLET incorporating a greater number of deuterium atoms than the apo and holo forms. Complementary proteolysis experiments revealed that HAMLET and apo are both accessible to proteases in the beta-domain but showed substantial differences in accessibility to proteases at specific sites. The overall results indicated that the conformational changes associated with the release of Ca2+ are not sufficient to induce the HAMLET conformation. Metal depletion might represent the first event to produce a partial unfolding in the beta-domain of alpha-lactalbumin, but some more unfolding is needed to generate the active conformation HAMLET, very likely allowing the protein to bind the C18:1 fatty acid moiety. On the basis of these data, a putative binding site of the oleic acid, which stabilizes the HAMLET conformation, is proposed.     

The inner membrane histidine kinase EnvZ senses osmolality via helix-coil transitions in the cytoplasm

Two-component systems mediate bacterial signal transduction, employing a membrane sensor kinase and a cytoplasmic response regulator (RR). Environmental sensing is typically coupled to gene regulation. Understanding how input stimuli activate kinase autophosphorylation remains obscure. The EnvZ/OmpR system regulates expression of outer membrane proteins in response to osmotic stress. To identify EnvZ conformational changes associated with osmosensing, we used HDXMS to probe the effects of osmolytes (NaCl, sucrose) on the cytoplasmic domain of EnvZ (EnvZ(c)). Increasing osmolality decreased deuterium exchange localized to the four-helix bundle containing the autophosphorylation site (His(243)). EnvZ(c) exists as an ensemble of multiple conformations and osmolytes favoured increased helicity. High osmolality increased autophosphorylation of His(243), suggesting that these two events are linked. In-vivo analysis showed that the cytoplasmic domain of EnvZ was sufficient for osmosensing, transmembrane domains were not required. Our results challenge existing claims of robustness in EnvZ/OmpR and support a model where osmolytes promote intrahelical H-bonding enhancing helix stabilization, increasing autophosphorylation and downstream signalling. The model provides a conserved mechanism for signalling proteins that respond to diverse physical and mechanical stimuli.     

The use of H/D exchange for secondary structure characterization of supermetallized complexes of ubiquitin with cerium(III)

The approach of hydrogen/deuterium exchange combined with ultrahigh resolution mass spectrometry was applied for investigation of conformational changes of supermetallized ubiquitin ions with cerium(III) atoms. The dependencies of the hydrogen/deuterium exchange efficiency on the charge state of ubiquitin ion, the number of associated cerium atoms, as well as on the temperature were obtained. The reaction of hydrogen/deuterium exchange was performed directly in the ionization source according to previously described method. It was found that the number of exchanges is hardly altered under the addition of cerium atoms. This result indirectly suggests that the conformation of small protein supermetallized ions does not significantly change during electrospray ionization.     

Influence of modulated structural dynamics on the kinetics of alpha-chymotrypsin catalysis. Insights through chemical glycosylation, molecular dynamics and domain motion analysis

Although the chemical nature of the catalytic mechanism of the serine protease alpha-chymotrypsin (alpha-CT) is largely understood, the influence of the enzyme's structural dynamics on its catalysis remains uncertain. Here we investigate whether alpha-CT's structural dynamics directly influence the kinetics of enzyme catalysis. Chemical glycosylation [Solá RJ & Griebenow K (2006) FEBS Lett 580, 1685-1690] was used to generate a series of glycosylated alpha-CT conjugates with reduced structural dynamics, as determined from amide hydrogen/deuterium exchange kinetics (k(HX)). Determination of their catalytic behavior (K(S), k(2), and k(3)) for the hydrolysis of N-succinyl-Ala-Ala-Pro-Phe p-nitroanilide (Suc-Ala-Ala-Pro-Phe-pNA) revealed decreased kinetics for the catalytic steps (k(2) and k(3)) without affecting substrate binding (K(S)) at increasing glycosylation levels. Statistical correlation analysis between the catalytic (DeltaG( not equal)k(i)) and structurally dynamic (DeltaG(HX)) parameters determined revealed that the enzyme acylation and deacylation steps are directly influenced by the changes in protein structural dynamics. Molecular modelling of the alpha-CT glycoconjugates coupled with molecular dynamics simulations and domain motion analysis employing the Gaussian network model revealed structural insights into the relation between the protein's surface glycosylation, the resulting structural dynamic changes, and the influence of these on the enzyme's collective dynamics and catalytic residues. The experimental and theoretical results presented here not only provide fundamental insights concerning the influence of glycosylation on the protein biophysical properties but also support the hypothesis that for alpha-CT the global structural dynamics directly influence the kinetics of enzyme catalysis via mechanochemical coupling between domain motions and active site chemical groups.