Hydrogen-exchange mass spectrometry reveals activation-induced changes in the conformational mobility of p38alpha MAP kinase

Conformational Dynamics Analysis of MEK1 Using Hydrogen/Deuterium Exchange Mass Spectrometry

BACKGROUND

Activation of mitogen-activated protein kinases (MAPKs) is regulated by a phosphorylation cascade comprising three kinases, MAPK kinase kinase (MAP3K), MAPK kinase (MAP2K), and MAPK. MAP2K1 and MAPK2K2, also known as MEK1 and MEK2, activate ERK1 and ERK2. The structure of the MAPK signaling cascade has been studied, but high-resolution structural studies of MAP2Ks have often focused on kinase domains or docking sites, but not on full-length proteins.

OBJECTIVE

To understand the conformational dynamics of MEK1.

METHODS

Full-length MEK1 was purified from Escherichia coli (BL21), and its conformational dynamics were analyzed using hydrogen/deuterium exchange mass spectrometry (HDX-MS). The effects of ATP binding were examined by co-incubating MEK1 and adenylyl-imidodiphosphate (AMP- PNP), a non-hydrolysable ATP analog.

RESULTS

MEK1 exhibited mixed EX1/EX2 HDX kinetics within the N-terminal tail through β1, αI, and the C-terminal helix. AMP-PNP binding was found to reduce conformational dynamics within the glycine-rich loop and regions near the DFG motif, along with the activation lip.

CONCLUSION

We report for the first time that MEK1 has regions that slowly change its folded and unfolded states (mixed EX1/EX2 kinetics) and also report the conformational effects of ATP-binding to MEK1.

Conformational selection vs. induced fit: insights into the binding mechanisms of p38α MAP Kinase inhibitors

The conformational dynamics of a kinase's activation loop have been challenging to assess due to the activation loop's intrinsic flexibility. To directly probe the conformational equilibrium of the activation loop of mitogen-activated protein kinase p38α, we present an approach based on site-directed spin labeling, electron paramagnetic resonance (EPR) distance restraints, and multilateration. We demonstrate that the activation loop of apo p38α resides in a highly flexible equilibrium state and we reveal that binding of small molecules significantly alters this equilibrium and the populated sub-states.

Structural Insights into Pseudokinase Domains of Receptor Tyrosine Kinases

Despite their apparent lack of catalytic activity, pseudokinases are essential signaling molecules. Here, we describe the structural and dynamic properties of pseudokinase domains from the Wnt-binding receptor tyrosine kinases (PTK7, ROR1, ROR2, and RYK), which play important roles in development. We determined structures of all pseudokinase domains in this family and found that they share a conserved inactive conformation in their activation loop that resembles the autoinhibited insulin receptor kinase (IRK). They also have inaccessible ATP-binding pockets, occluded by aromatic residues that mimic a cofactor-bound state. Structural comparisons revealed significant domain plasticity and alternative interactions that substitute for absent conserved motifs. The pseudokinases also showed dynamic properties that were strikingly similar to those of IRK. Despite the inaccessible ATP site, screening identified ATP-competitive type-II inhibitors for ROR1. Our results set the stage for an emerging therapeutic modality of "conformational disruptors" to inhibit or modulate non-catalytic functions of pseudokinases deregulated in disease.

p38α Mitogen-Activated Protein Kinase Is a Druggable Target in Pancreatic Adenocarcinoma

p38 mitogen-activated protein kinases are signaling molecules with major involvement in cancer. A detailed mechanistic understanding of how p38 MAPK family members function is urgently warranted for cancer targeted therapy. The conformational dynamics of the most common member of p38 MAPK family, p38α, are crucial for its function but poorly understood. Here we found that, unlike in other cancer types, p38α is significantly activated in pancreatic adenocarcinoma samples, suggesting its potential for anti-pancreatic cancer therapy. Using a state of the art supercomputer, Anton, long-timescale (39 μs) unbiased molecular dynamics simulations of p38α show that apo p38α has high structural flexibility in six regions, and reveal potential catalysis mechanism involving a “butterfly” motion. Moreover, in vitro studies show the low-selectivity of the current p38α inhibitors in both human and mouse pancreatic cancer cell lines, while computational solvent mapping identified 17 novel pockets for drug design. Taken together, our study reveals the conformational dynamics and potentially druggable pockets of p38α, which may potentiate p38α-targeting drug development and benefit pancreatic cancer patients.

A temperature-dependent conformational shift in p38α MAPK substrate-binding region associated with changes in substrate phosphorylation profile

Febrile-range hyperthermia worsens and hypothermia mitigates lung injury, and temperature dependence of lung injury is blunted by inhibitors of p38 mitogen-activated protein kinase (MAPK). Of the two predominant p38 isoforms, p38α is proinflammatory and p38β is cytoprotective. Here, we analyzed the temperature dependence of p38 MAPK activation, substrate interaction, and tertiary structure. Incubating HeLa cells at 39.5 °C stimulated modest p38 activation, but did not alter tumor necrosis factor-α (TNFα)-induced p38 activation. In in vitro kinase assays containing activated p38α and MAPK-activated kinase-2 (MK2), MK2 phosphorylation was 14.5-fold greater at 39.5 °C than at 33 °C. By comparison, we observed only 3.1- and 1.9-fold differences for activating transcription factor-2 (ATF2) and signal transducer and activator of transcription-1α (STAT1α) and a 7.7-fold difference for p38β phosphorylation of MK2. The temperature dependence of p38α:substrate binding affinity, as measured by surface plasmon resonance, paralleled substrate phosphorylation. Hydrogen-deuterium exchange MS (HDX-MS) of p38α performed at 33, 37, and 39.5 °C indicated temperature-dependent conformational changes in an α helix near the common docking and glutamate:aspartate substrate-binding domains at the known binding site for MK2. In contrast, HDX-MS analysis of p38β did not detect significant temperature-dependent conformational changes in this region. We observed no conformational changes in the catalytic domain of either isoform and no corresponding temperature dependence in the C-terminal p38α-interacting region of MK2. Because MK2 participates in the pathogenesis of lung injury, the observed changes in the structure and function of proinflammatory p38α may contribute to the temperature dependence of acute lung injury.

Kinetic network model to explain gain-of-function mutations in ERK2 enzyme

ERK2 is a kinase protein that belongs to a Ras/Raf/MEK/ERK signaling pathway, which is activated in response to a range of extracellular signals. Malfunctioning of this cascade leads to a variety of serious diseases, including cancers. This is often caused by mutations in proteins belonging to the cascade, frequently leading to abnormally high activity of the cascade even in the absence of an external signal. One such "gain-of-function" mutation in the ERK2 protein, called a "sevenmaker" mutation (D319N), was discovered in 1994 in Drosophila. The mutation leads to disruption of interactions of other proteins with the D-site of ERK2 and results, contrary to expectations, in an increase of its activity in vivo. However, no molecular mechanism to explain this effect has been presented so far. The difficulty is that this mutation should equally negatively affect interactions of ERK2 with all substrates, activators, and deactivators. In this paper, we present a semiquantitative kinetic network model that gives a possible explanation of the increased activity of mutant ERK2 species. A simplified biochemical network for ERK2, viewed as a system of coupled Michaelis-Menten processes, is presented. Its dynamic properties are calculated explicitly using the method of first-passage processes. The effect of mutation is associated with changes in the strength of interaction energy between the enzyme and the substrates. It is found that the dependence of kinetic properties of the protein on the interaction energy is nonmonotonic, suggesting that some mutations might lead to more efficient catalytic properties, despite weakening intermolecular interactions. Our theoretical predictions agree with experimental observations for the sevenmaker mutation in ERK2. It is also argued that the effect of mutations might depend on the concentrations of substrates.

Substrate Locking Promotes Dimer-Dimer Docking of an Enzyme Antibiotic Target

Protein dynamics manifested through structural flexibility play a central role in the function of biological molecules. Here we explore the substrate-mediated change in protein flexibility of an antibiotic target enzyme, Clostridium botulinum dihydrodipicolinate synthase. We demonstrate that the substrate, pyruvate, stabilizes the more active dimer-of-dimers or tetrameric form. Surprisingly, there is little difference between the crystal structures of apo and substrate-bound enzyme, suggesting protein dynamics may be important. Neutron and small-angle X-ray scattering experiments were used to probe substrate-induced dynamics on the sub-second timescale, but no significant changes were observed. We therefore developed a simple technique, coined protein dynamics-mass spectrometry (ProD-MS), which enables measurement of time-dependent alkylation of cysteine residues. ProD-MS together with X-ray crystallography and analytical ultracentrifugation analyses indicates that pyruvate locks the conformation of the dimer that promotes docking to the more active tetrameric form, offering insight into ligand-mediated stabilization of multimeric enzymes.

Structural characterization of the catalytic γ and regulatory β subunits of phosphorylase kinase in the context of the hexadecameric enzyme complex

In the tightly regulated glycogenolysis cascade, the breakdown of glycogen to glucose-1-phosphate, phosphorylase kinase (PhK) plays a key role in regulating the activity of glycogen phosphorylase. PhK is a 1.3 MDa hexadecamer, with four copies each of four different subunits (α, β, γ and δ), making the study of its structure challenging. Using hydrogen-deuterium exchange, we have analyzed the regulatory β subunit and the catalytic γ subunit in the context of the intact non-activated PhK complex to study the structure of these subunits and identify regions of surface exposure. Our data suggest that within the non-activated complex the γ subunit assumes an activated conformation and are consistent with a previous docking model of the β subunit within the cryoelectron microscopy envelope of PhK.

The Hippo Pathway Maintains the Equatorial Division Plane in the Ciliate Tetrahymena

The mechanisms that govern pattern formation within the cell are poorly understood. Ciliates carry on their surface an elaborate pattern of cortical organelles that are arranged along the anteroposterior and circumferential axes by largely unknown mechanisms. Ciliates divide by tandem duplication: the cortex of the predivision cell is remodeled into two similarly sized and complete daughters. In the conditional cdaI-1 mutant of Tetrahymena thermophila, the division plane migrates from its initially correct equatorial position toward the cell's anterior, resulting in unequal cell division, and defects in nuclear divisions and cytokinesis. We used comparative whole genome sequencing to identify the cause of cdaI-1 as a mutation in a Hippo/Mst kinase. CdaI is a cortical protein with a cell cycle-dependent, highly polarized localization. Early in cell division, CdaI marks the anterior half of the cell, and later concentrates at the posterior end of the emerging anterior daughter. Despite the strong association of CdaI with the new posterior cell end, the cdaI-1 mutation does not affect the patterning of the new posterior cortical organelles. We conclude that, in Tetrahymena, the Hippo pathway maintains an equatorial position of the fission zone, and, by this activity, specifies the relative dimensions of the anterior and posterior daughter cell.

Changes in the free-energy landscape of p38α MAP kinase through its canonical activation and binding events as studied by enhanced molecular dynamics simulations

p38α is a Ser/Thr protein kinase involved in a variety of cellular processes and pathological conditions, which makes it a promising pharmacological target. Although the activity of the enzyme is highly regulated, its molecular mechanism of activation remains largely unexplained, even after decades of research. By using state-of-the-art molecular dynamics simulations, we decipher the key elements of the complex molecular mechanism refined by evolution to allow for a fine tuning of p38α kinase activity. Our study describes for the first time the molecular effects of different regulators of the enzymatic activity, and provides an integrative picture of the activation mechanism that explains the seemingly contradictory X-ray and NMR data.

DOI:http://dx.doi.org/10.7554/eLife.22175.001

Conformations of JNK3α splice variants analyzed by hydrogen/deuterium exchange mass spectrometry

c-Jun N-terminal kinases (JNKs) are members of the mitogen-activated protein kinase (MAPK) family that regulate apoptosis, inflammation, cytokine production, and metabolism. MAPKs undergo various splicing within their kinase domains. Unlike other MAPKs, JNKs have alternative splicing at the C-terminus, resulting in long and short variants. Functional or conformational effects due to the elongated C-terminal tail in the long splice variants have not been investigated nor has the conformation of the C-terminal tail been analyzed. Here, we analyzed the conformation of the elongated C-terminal tail and investigated conformational differences between long and short splice variants of JNKs using JNK3α2 and JNK3α1 as models. We adopted hydrogen/deuterium exchange mass spectrometry (HDX-MS) to analyze the conformation. HDX-MS revealed that the C-terminal tail is mostly intrinsically disordered, and that the conformation of the kinase domain of JNK3α2 is more dynamic than that of JNK3α1. The different conformation dynamics between long and short splice variants of JNK3α might affect the cellular functions of JNK3.

Theoretical Investigation of the Mechanisms of ERK2 Enzymatic Catalysis

ERK2 are protein kinases that during the enzymatic catalysis, in contrast to traditional enzymes, utilize additional interactions with substrates outside of the active sites. It is widely believed that these docking interactions outside of the enzymatic pockets enhance the specificity of these proteins. However, the molecular mechanisms of how the docking interactions affect the catalysis remain not well understood. Here, we develop a simple theoretical approach to analyze the enzymatic catalysis in ERK2 proteins. Our method is based on first-passage process analysis, and it provides explicit expressions for all dynamic properties of the system. It is found that there are specific binding energies for substrates in docking and catalytic domains that lead to maximal enzymatic reaction rates. Thus, we propose that the role of the docking interactions is not only to increase the enzymatic specificity but also to optimize the dynamics of the catalytic process. Our theoretical results are utilized to describe experimental observations on ERK2 enzymatic activities.

p38β Mitogen-Activated Protein Kinase Modulates Its Own Basal Activity by Autophosphorylation of the Activating Residue Thr180 and the Inhibitory Residues Thr241 and Ser261

Many enzymes are self-regulated and can either inhibit or enhance their own catalytic activity. Enzymes that do both are extremely rare. Many protein kinases autoactivate by autophosphorylating specific sites at their activation loop and are inactivated by phosphatases. Although mitogen-activated protein kinases (MAPKs) are usually activated by dual phosphorylation catalyzed by MAPK kinases (MAPKKs), the MAPK p38β is exceptional and is capable of self-activation by cis autophosphorylation of its activation loop residue T180. We discovered that p38β also autophosphorylates in trans two previously unknown sites residing within a MAPK-specific structural element known as the MAPK insert: T241 and S261. Whereas phosphorylation of T180 evokes catalytic activity, phosphorylation of S261 reduces the activity of T180-phosphorylated p38β, and phosphorylation of T241 reduces its autophosphorylation in trans Both phosphorylations do not affect the activity of dually phosphorylated p38β. T241 of p38β is found phosphorylated in vivo in bone and muscle tissues. In myogenic cell lines, phosphorylation of p38β residue T241 is correlated with differentiation to myotubes. T241 and S261 are also autophosphorylated in intrinsically active variants of p38α, but in this protein, they probably play a different role. We conclude that p38β is an unusual enzyme that automodulates its basal, MAPKK-independent activity by several autophosphorylation events, which enhance and suppress its catalytic activity.

Tighter αC-helix-αL16-helix interactions seem to make p38α less prone to activation by autophosphorylation than Hog1

A structural element termed ‘hydrophobic core’ is a suppressor of spontaneous autophosphorylation in Hog1 and p38s. Practically any mutation in this core of Hog1, but not of p38, evokes spontaneous autophosphorylation. This inherent autophosphorylation suppressor is tighter in mammalian's p38s.

Many eukaryotic protein kinases (EPKs) are autoactivated through autophosphorylation of their activation loop. Mitogen-activated protein (MAP) kinases do not autophosphorylate spontaneously; relying instead upon mitogen-activated protein kinase (MAPK) kinases (MKKs) for their activation loop phosphorylation. Yet, in previous studies we identified mutations in the yeast MAPK high osmolarity glycerol (Hog1) that render it capable of spontaneous autophosphorylation and consequently intrinsically active (MKK-independent). Four of the mutations occurred in hydrophobic residues, residing in the αC-helix, which is conserved in all EPKs, and in the αL16-helix that is unique to MAPKs. These four residues interact together forming a structural element termed ‘hydrophobic core’. A similar element exists in the Hog1’s mammalian orthologues p38s. Here we show that the ‘hydrophobic core’ is a loose suppressor of Hog1’s autophosphorylation. We inserted 18 point mutations into this core, 17 of which were able to render Hog1 MKK-independent. In p38s, however, only a very few mutations in the equivalent residues rendered these proteins intrinsically active. Structural analysis revealed that a salt bridge between the αC-helix and the αL16-helix that exists in p38α may not exist in Hog1. This bond further stabilizes the ‘hydrophobic core’ of p38, making p38 less prone to de-repressing its concealed autophosphorylation. Mutating equivalent hydrophobic residues in Jnk1 and Erk2 has no effect on their autophosphorylation. We propose that specific structural elements developed in the course of evolution to suppress spontaneous autophosphorylation of Hog1/p38. The suppressors were kept wobbly, probably to allow activation by induced autophosphorylation, but became stricter in mammalian p38s than in the yeast Hog1.

Phosphorylation in the catalytic cleft stabilizes and attracts domains of a phosphohexomutase

Phosphorylation can modulate the activities of enzymes. The phosphoryl donor in the catalytic cleft of α-D-phosphohexomutases is transiently dephosphorylated while the reaction intermediate completes a 180° reorientation within the cleft. The phosphorylated form of 52 kDa bacterial phosphomannomutase/phosphoglucomutase is less accessible to dye or protease, more stable to chemical denaturation, and widely stabilized against NMR-detected hydrogen exchange across the core of domain 3 to juxtaposed domain 4 (each by ≥ 1.3 kcal/mol) and parts of domains 1 and 2. However, phosphorylation accelerates hydrogen exchange in specific regions of domains 1 and 2, including a metal-binding residue in the active site. Electrostatic field lines reveal attraction across the catalytic cleft between phosphorylated Ser-108 and domain 4, but repulsion when Ser-108 is dephosphorylated. Molecular dynamics (MD) simulated the dephosphorylated form to be expanded due to enhanced rotational freedom of domain 4. The contacts and fluctuations of the MD trajectories enabled correct simulation of more than 80% of sites that undergo either protection or deprotection from hydrogen exchange due to phosphorylation. Electrostatic attraction in the phosphorylated enzyme accounts for 1) domain 4 drawing closer to domains 1 and 3; 2) decreased accessibility; and 3) increased stability within these domains. The electrostriction due to phosphorylation may help capture substrate, whereas the opening of the cleft upon transient dephosphorylation allows rotation of the intermediate. The long-range effects of phosphorylation on hydrogen exchange parallel reports on protein kinases, suggesting a conceptual link among these multidomain, phosphoryl transfer enzymes.

Phosphorylation- and nucleotide-binding-induced changes to the stability and hydrogen exchange patterns of JNK1β1 provide insight into its mechanisms of activation

Many studies have characterized how changes to the stability and internal motions of a protein during activation can contribute to their catalytic function, even when structural changes cannot be observed. Here, unfolding studies and hydrogen-deuterium exchange (HX) mass spectrometry were used to investigate the changes to the stability and conformation/conformational dynamics of JNK1β1 induced by phosphorylative activation. Equivalent studies were also employed to determine the effects of nucleotide binding on both inactive and active JNK1β1 using the ATP analogue, 5'-adenylyl-imidodiphosphate (AMP-PNP). JNK1β1 phosphorylation alters HX in regions involved in catalysis and substrate binding, changes that can be ascribed to functional modifications in either structure and/or backbone flexibility. Increased HX in the hinge between the N- and C-terminal domains implied that it acquires enhanced flexibility upon phosphorylation that may be a prerequisite for interdomain closure. In combination with the finding that nucleotide binding destabilizes the kinase, the patterns of solvent protection by AMP-PNP were consistent with a novel mode of nucleotide binding to the C-terminal domain of a destabilized and open domain conformation of inactive JNK1β1. Solvent protection by AMP-PNP of both N- and C-terminal domains in active JNK1β1 revealed that the domains close around nucleotide upon phosphorylation, concomitantly stabilizing the kinase. This suggests that phosphorylation activates JNK1β1 in part by increasing hinge flexibility to facilitate interdomain closure and the creation of a functional active site. By uncovering the complex interplay that occurs between nucleotide binding and phosphorylation, we present new insight into the unique mechanisms by which JNK1β1 is regulated.

The p38β mitogen-activated protein kinase possesses an intrinsic autophosphorylation activity, generated by a short region composed of the α-G helix and MAPK insert

Protein kinases are regulated by a large number of mechanisms that vary from one kinase to another. However, a fundamental activation mechanism shared by all protein kinases is phosphorylation of a conserved activation loop threonine residue. This is achieved in many cases via autophosphorylation. The mechanism and structural basis for autophosphorylation are not clear and are in fact enigmatic because this phosphorylation occurs when the kinase is in its inactive conformation. Unlike most protein kinases, MAP kinases are not commonly activated by autophosphorylation but rather by MEK-dependent phosphorylation. Here we show that p38β, a p38 isoform that is almost identical to p38α, is exceptional and spontaneously autoactivates by autophosphorylation. We identified a 13-residue-long region composed of part of the αG-helix and the MAPK insert that triggers the intrinsic autophosphorylation activity of p38β. When inserted into p38α, this fragment renders it spontaneously active in vitro and in mammalian cells. We further found that an interaction between the N terminus and a particular region of the C-terminal extension suppresses the intrinsic autophosphorylation of p38β in mammalian cells. Thus, this study identified the structural motif responsible for the unique autophosphorylation capability of p38β and the motif inhibiting this activity in living cells. It shows that the MAPK insert and C-terminal extension, structural motifs that are unique to MAPKs, play a critical role in controlling autophosphorylation.

Extracellular-regulated kinase 2 is activated by the enhancement of hinge flexibility

Protein motions underlie conformational and entropic contributions to enzyme catalysis; however, relatively little is known about the ways in which this occurs. Studies of the mitogen-activated protein kinase ERK2 (extracellular-regulated protein kinase 2) by hydrogen-exchange mass spectrometry suggest that activation enhances backbone flexibility at the linker between N- and C-terminal domains while altering nucleotide binding mode. Here, we address the hypothesis that enhanced backbone flexibility within the hinge region facilitates kinase activation. We show that hinge mutations enhancing flexibility promote changes in the nucleotide binding mode consistent with domain movement, without requiring phosphorylation. They also lead to the activation of monophosphorylated ERK2, a form that is normally inactive. The hinge mutations bypass the need for pTyr but not pThr, suggesting that Tyr phosphorylation controls hinge motions. In agreement, monophosphorylation of pTyr enhances both hinge flexibility and nucleotide binding mode, measured by hydrogen-exchange mass spectrometry. Our findings demonstrate that regulated protein motions underlie kinase activation. Our working model is that constraints to domain movement in ERK2 are overcome by phosphorylation at pTyr, which increases hinge dynamics to promote the active conformation of the catalytic site.

Promotion of enzyme flexibility by dephosphorylation and coupling to the catalytic mechanism of a phosphohexomutase

The enzyme phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa catalyzes an intramolecular phosphoryl transfer across its phosphosugar substrates, which are precursors in the synthesis of exoproducts involved in bacterial virulence. Previous structural studies of PMM/PGM have established a key role for conformational change in its multistep reaction, which requires a dramatic 180° reorientation of the intermediate within the active site. Here hydrogen-deuterium exchange by mass spectrometry and small angle x-ray scattering were used to probe the conformational flexibility of different forms of PMM/PGM in solution, including its active, phosphorylated state and the unphosphorylated state that occurs transiently during the catalytic cycle. In addition, the effects of ligand binding were assessed through use of a substrate analog. We found that both phosphorylation and binding of ligand produce significant effects on deuterium incorporation. Phosphorylation of the conserved catalytic serine has broad effects on residues in multiple domains and is supported by small angle x-ray scattering data showing that the unphosphorylated enzyme is less compact in solution. The effects of ligand binding are generally manifested near the active site cleft and at a domain interface that is a site of conformational change. These results suggest that dephosphorylation of the enzyme may play two critical functional roles: a direct role in the chemical step of phosphoryl transfer and secondly through propagation of structural flexibility. We propose a model whereby increased enzyme flexibility facilitates the reorientation of the reaction intermediate, coupling changes in structural dynamics with the unique catalytic mechanism of this enzyme.

Importance of protein dynamics during enzymatic C-H bond cleavage catalysis

Quantum tunneling and protein dynamics have emerged as important components of enzyme function. This review focuses on soybean lipoxygenase-1, to illustrate how the properties of enzymatic C-H bond activation link protein motions to the fundamental bond making-breaking processes.

Integration of phosphoproteomic, chemical, and biological strategies for the functional analysis of targeted protein phosphorylation

Reversible phosphorylation, tightly controlled by protein kinases and phosphatases, plays a central role in mediating biological processes, such as protein-protein interactions, subcellular translocation, and activation of cellular enzymes. MS-based phosphoproteomics has now allowed the detection and quantification of tens of thousands of phosphorylation sites from a typical biological sample in a single experiment, which has posed new challenges in functional analysis of each and every phosphorylation site on specific signaling phosphoproteins of interest. In this article, we review recent advances in the functional analysis of targeted phosphorylation carried out by various chemical and biological approaches in combination with the MS-based phosphoproteomics. This review focuses on three types of strategies, including forward functional analysis, defined for the result-driven phosphoproteomics efforts in determining the substrates of a specific protein kinase; reverse functional analysis, defined for tracking the kinase(s) for specific phosphosite(s) derived from the discovery-driven phosphoproteomics efforts; and MS-based analysis on the structure-function relationship of phosphoproteins. It is expected that this review will provide a state-of-the-art overview of functional analysis of site-specific phosphorylation and explore new perspectives and outline future challenges.

Hydrogen tunneling links protein dynamics to enzyme catalysis

The relationship between protein dynamics and function is a subject of considerable contemporary interest. Although protein motions are frequently observed during ligand binding and release steps, the contribution of protein motions to the catalysis of bond making/breaking processes is more difficult to probe and verify. Here, we show how the quantum mechanical hydrogen tunneling associated with enzymatic C-H bond cleavage provides a unique window into the necessity of protein dynamics for achieving optimal catalysis. Experimental findings support a hierarchy of thermodynamically equilibrated motions that control the H-donor and -acceptor distance and active-site electrostatics, creating an ensemble of conformations suitable for H-tunneling. A possible extension of this view to methyl transfer and other catalyzed reactions is also presented. The impact of understanding these dynamics on the conceptual framework for enzyme activity, inhibitor/drug design, and biomimetic catalyst design is likely to be substantial.

Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.

Conformational dynamics of a membrane transport protein probed by H/D exchange and covalent labeling: the glycerol facilitator

Glycerol facilitator (GF) is a tetrameric membrane protein responsible for the selective permeation of glycerol and water. Each of the four GF subunits forms a transmembrane channel. Every subunit consists of six helices that completely span the lipid bilayer, as well as two half-helices (TM7 and TM3). X-ray crystallography has revealed that the selectivity of GF is due to its unique amphipathic channel interior. To explore the structural dynamics of GF, we employ hydrogen/deuterium exchange (HDX) and oxidative labeling with mass spectrometry (MS). HDX-MS reveals that transmembrane helices are generally more protected than extramembrane segments, consistent with data previously obtained for other membrane proteins. Interestingly, TM7 does not follow this trend. Instead, this half-helix undergoes rapid deuteration, indicative of a highly dynamic local structure. The oxidative labeling behavior of most GF residues is consistent with the static crystal structure. However, the side chains of C134 and M237 undergo labeling although they should be inaccessible according to the X-ray structure. In agreement with our HDX-MS data, this observation attests to the fact that TM7 is only marginally stable. We propose that the highly mobile nature of TM7 aids in the efficient diffusion of guest molecules through the channel ("molecular lubrication"). In the absence of such dynamics, host-guest molecular recognition would favor semipermanent binding of molecules inside the channel, thereby impeding transport. The current work highlights the complementary nature of HDX, covalent labeling, and X-ray crystallography for the characterization of membrane proteins.

Comparative hydrogen-deuterium exchange for a mesophilic vs thermophilic dihydrofolate reductase at 25 °C: identification of a single active site region with enhanced flexibility in the mesophilic protein

The technique of hydrogen–deuterium exchange coupled to mass spectrometry (HDX-MS) has been applied to a mesophilic (E. coli) dihydrofolate reductase under conditions that allow direct comparison to a thermophilic (B. stearothermophilus) ortholog, Ec-DHFR and Bs-DHFR, respectively. The analysis of hydrogen–deuterium exchange patterns within proteolytically derived peptides allows spatial resolution, while requiring a series of controls to compare orthologous proteins with only ca. 40% sequence identity. These controls include the determination of primary structure effects on intrinsic rate constants for HDX as well as the use of existing 3-dimensional structures to evaluate the distance of each backbone amide hydrogen to the protein surface. Only a single peptide from the Ec-DHFR is found to be substantially more flexible than the Bs-DHFR at 25 °C in a region located within the protein interior at the intersection of the cofactor and substrate-binding sites. The surrounding regions of the enzyme are either unchanged or more flexible in the thermophilic DHFR from B. stearothermophilus. The region with increased flexibility in Ec-DHFR corresponds to one of two regions previously proposed to control the enthalpic barrier for hydride transfer in Bs-DHFR [Oyeyemi et al. (2010) Proc. Natl. Acad. Sci. U.S.A.50, 1007420534574].

Distinct patterns of activation-dependent changes in conformational mobility between ERK1 and ERK2

Hydrogen/deuterium exchange measurements by mass spectrometry (HX-MS) can be used to report localized conformational mobility within folded proteins, where exchange predominantly occurs through low energy fluctuations in structure, allowing transient solvent exposure. Changes in conformational mobility may impact protein function, even in cases where structural changes are unobservable. Previous studies of the MAP kinase, ERK2, revealed increases in HX upon activation occured at the hinge between conserved N- and C-terminal domains, which could be ascribed to enhanced backbone flexibility. This implied that kinase activation modulates interdomain closure, and was supported by evidence for two modes of nucleotide binding that were consistent with closed vs open conformations in active vs inactive forms of ERK2, respectively. Thus, phosphorylation of ERK2 releases constraints to interdomain closure, by modulating hinge flexibility. In this study, we examined ERK1, which shares 90% sequence identity with ERK2. HX-MS measurements of ERK1 showed similarities with ERK2 in overall deuteration, consistent with their similar tertiary structures. However, the patterns of HX that were altered upon activation of ERK1 differed from those in ERK2. In particular, alterations in HX at the hinge region upon activation of ERK2 did not occur in ERK1, suggesting that the two enzymes differ with respect to their regulation of hinge mobility and interdomain closure. In agreement, HX-MS measurements of nucleotide binding suggested revealed domain closure in both inactive and active forms of ERK1. We conclude that although ERK1 and ERK2 are closely related with respect to primary sequence and tertiary structure, they utilize distinct mechanisms for controlling enzyme function through interdomain interactions.

Carboxyl-group footprinting maps the dimerization interface and phosphorylation-induced conformational changes of a membrane-associated tyrosine kinase

Her4 is a transmembrane receptor tyrosine kinase belonging to the ErbB-EGFR family. It plays a vital role in the cardiovascular and nervous systems, and mutations in Her4 have been found in melanoma and lung cancer. The kinase domain of Her4 forms a dimer complex, called the asymmetric dimer, which results in kinase activation. Although a crystal structure of the Her4 asymmetric dimer is known, the dimer affinity and the effect of the subsequent phosphorylation steps on kinase domain conformation are unknown. We report here the use of carboxyl-group footprinting MS on a recombinant expressed, Her4 kinase-domain construct to address these questions. Carboxyl-group footprinting uses a water-soluble carbodiimide, 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide, in the presence of glycine ethyl ester, to modify accessible carboxyl groups on glutamate and aspartate residues. Comparisons of Her4 kinase-domain monomers versus dimers and of unphosphorylated versus phosphorylated dimers were made to map the dimerization interface and to determine phosphorylation induced-conformational changes. We detected 37 glutamate and aspartate residues that were modified, and we quantified their extents of modification by liquid chromatography MS. Five residues showed changes in carboxyl-group modification. Three of these residues are at the predicted dimer interface, as shown by the crystal structure, and the remaining two residues are on loops that likely have altered conformation in the kinase dimer. Incubating the Her4 kinase dimers with ATP resulted in dramatic increase in Tyr-850 phosphorylation, located on the activation loop, and this resulted in a conformational change in this loop, as evidenced by reduction in carboxyl-group modification. The kinase monomer-dimer equilibrium was measured using a titration format in which the extent of carboxyl-group footprinting was mathematically modeled to give the dimer association constant (1.5-6.8 × 10(12) dm(2)/mol). This suggests that the kinase-domain makes a significant contribution to the overall dimerization affinity of the full-length Her4 protein.

p38γ activation triggers dynamical changes in allosteric docking sites

Mitogen-activated protein kinases (MAPKs) are serine-threonine kinases that participate in signal transduction pathways. p38 MAPKs have four isoforms (p38α, p38β, p38γ, and p38δ) which are involved in multiple cellular functions such as proliferation, differentiation, survival, and migration. MAPK kinases phosphorylate p38s in the dual-phosphorylation motif, Thr-Gly-Tyr, located in their activation loop, which induces a conformational change that increases ATP binding affinity and catalytic activity. Several works have proposed that MAPK dynamics is a key factor in determining their function. However, we still do not understand the dynamical changes that lead to MAPK activation. In this work we have used molecular dynamics techniques to study the dynamical changes associated with p38γ activation, the only fully active MAPK crystallized so far. We performed MD simulations of p38γ in three different states, fully active with ATP, active without ATP, and inactive. We found that the dynamical fluctuations of the docking sites, important for protein-protein interactions, are regulated allosterically by changes in the active site. Interestingly, in the phosphorylated and ATP-bound states the whole protein dynamics lead to concerted motions of whole protein domains in contrast to the inactive state. The binding/unbinding of ATP participates in the reorientation of the two domains and in the regulation of protein plasticity. Our study shows that beyond the conformational changes associated with MAPK activation their correlated dynamics are highly regulated by phosphorylation and ATP binding. This means that MAPK plasticity may have a role in their catalytic activity, specificity, and protein-protein interactions and, therefore, in the outcome of the signaling network.

Cyclic AMP- and (Rp)-cAMPS-induced conformational changes in a complex of the catalytic and regulatory (RI{alpha}) subunits of cyclic AMP-dependent protein kinase

We took a discovery approach to explore the actions of cAMP and two of its analogs, one a cAMP mimic ((S(p))-adenosine cyclic 3':5'-monophosphorothioate ((S(p))-cAMPS)) and the other a diastereoisomeric antagonist ((R(p))-cAMPS), on a model system of the type Iα cyclic AMP-dependent protein kinase holoenzyme, RIα(91-244)·C-subunit, by using fluorescence spectroscopy and amide H/(2)H exchange mass spectrometry. Specifically, for the fluorescence experiments, fluorescein maleimide was conjugated to three cysteine single residue substitution mutants, R92C, T104C, and R239C, of RIα(91-244), and the effects of cAMP, (S(p))-cAMPS, and (R(p))-cAMPS on the kinetics of R-C binding and the time-resolved anisotropy of the reporter group at each conjugation site were measured. For the amide exchange experiments, ESI-TOF mass spectrometry with pepsin proteolytic fragmentation was used to assess the effects of (R(p))-cAMPS on amide exchange of the RIα(91-244)·C-subunit complex. We found that cAMP and its mimic perturbed at least parts of the C-subunit interaction Sites 2 and 3 but probably not Site 1 via reduced interactions of the linker region and αC of RIα(91-244). Surprisingly, (R(p))-cAMPS not only increased the affinity of RIα(91-244) toward the C-subunit by 5-fold but also produced long range effects that propagated through both the C- and R-subunits to produce limited unfolding and/or enhanced conformational flexibility. This combination of effects is consistent with (R(p))-cAMPS acting by enhancing the internal entropy of the R·C complex. Finally, the (R(p))-cAMPS-induced increase in affinity of RIα(91-244) toward the C-subunit indicates that (R(p))-cAMPS is better described as an inverse agonist because it decreases the fractional dissociation of the cyclic AMP-dependent protein kinase holoenzyme and in turn its basal activity.

Analysis of MAP kinases by hydrogen exchange mass spectrometry

Hydrogen exchange mass spectrometry (HX-MS) is an experimental technique that can be used to -examine solvent accessibility and conformational mobility in biological macromolecules. This chapter summarizes studies using HX-MS to examine the regulation of conformation, protein mobility, and ligand binding to MAP kinases. We describe the planning and design of HX-MS experiments, strategies for data analysis and interpretation, and available software.

Global consequences of activation loop phosphorylation on protein kinase A

Phosphorylation of the activation loop is one of the most common mechanisms for regulating protein kinase activity. The catalytic subunit of cAMP-dependent protein kinase autophosphorylates Thr(197) in the activation loop when expressed in Escherichia coli. Although mutation of Arg(194) to Ala prevents autophosphorylation, phosphorylation of Thr(197) can still be achieved by a heterologous protein kinase, phosphoinositide-dependent protein kinase (PDK1), in vitro. In this study, we examined the structural and functional consequences of adding a single phosphate to the activation loop of cAMP-dependent protein kinase by comparing the wild type C-subunit to the R194A mutant either in the presence or the absence of activation loop phosphorylation. Phosphorylation of Thr(197) decreased the K(m) by approximately 15- and 7-fold for kemptide and ATP, respectively, increased the stability of the enzyme as measured by fluorescence and circular dichroism, and enhanced the binding between the C-subunit and IP20, a protein kinase inhibitor peptide. Additionally, deuterium exchange coupled to mass spectrometry was used to compare the structural dynamics of these proteins. All of the regions of the C-subunit analyzed underwent amide hydrogen exchange at a higher or equal rate in the unphosphorylated enzyme compared with the phosphorylated enzyme. The largest changes occurred at the C terminus of the activation segment in the p + 1 loop/APE regions and the alphaH-alphaI loop motifs and leads to the prediction of a coordinated phosphorylation-induced salt bridge between two conserved residues, Glu(208) and Arg(280).

An Evaluation of Explicit Receptor Flexibility in Molecular Docking Using Molecular Dynamics and Torsion Angle Molecular Dynamics

Incorporating receptor flexibility into molecular docking should improve results for flexible proteins. However, the incorporation of explicit all-atom flexibility with molecular dynamics for the entire protein chain may also introduce significant error and "noise" that could decrease docking accuracy and deteriorate the ability of a scoring function to rank native-like poses. We address this apparent paradox by comparing the success of several flexible receptor models in cross-docking and multiple receptor ensemble docking for p38α mitogen-activated protein (MAP) kinase. Explicit all-atom receptor flexibility has been incorporated into a CHARMM-based molecular docking method (CDOCKER) using both molecular dynamics (MD) and torsion angle molecular dynamics (TAMD) for the refinement of predicted protein-ligand binding geometries. These flexible receptor models have been evaluated, and the accuracy and efficiency of TAMD sampling is directly compared to MD sampling. Several flexible receptor models are compared, encompassing flexible side chains, flexible loops, multiple flexible backbone segments, and treatment of the entire chain as flexible. We find that although including side chain and some backbone flexibility is required for improved docking accuracy as expected, docking accuracy also diminishes as additional and unnecessary receptor flexibility is included into the conformational search space. Ensemble docking results demonstrate that including protein flexibility leads to to improved agreement with binding data for 227 active compounds. This comparison also demonstrates that a flexible receptor model enriches high affinity compound identification without significantly increasing the number of false positives from low affinity compounds.

T cell receptor-mediated activation of p38{alpha} by mono-phosphorylation of the activation loop results in altered substrate specificity

p38 MAPKs are typically activated by upstream MAPK kinases that phosphorylate a Thr-X-Tyr motif in the activation loop. An exception is the T cell antigen receptor signaling pathway, which bypasses the MAPK cascade and activates p38alpha and p38beta by phosphorylation of Tyr-323 and subsequent autophosphorylation of the activation loop. Here we show that, unlike the classic MAPK cascade, the alternative pathway results primarily in mono-phosphorylation of the activation loop residue Thr-180. Recombinant mono-phosphorylated and dual phosphorylated p38alpha differed widely with regard to activity and substrate preference. Altered substrate specificity was reproduced in T cells in which p38 was activated by the alternative or classical MAPK pathways. These findings suggest that T cells have evolved a mechanism to utilize p38 in a specialized manner independent of and distinct from the classical p38 MAPK signaling cascade.

Energy-based analysis and prediction of the orientation between light- and heavy-chain antibody variable domains

Diversity in antibody structure is crucial to the ability of the adaptive immune system to recognize the tremendously diverse set of potential antigens. The diversity in structure is most apparent in the six hypervariable loops of the complementarity-determining regions. However, given that these loops occur at the interface of the heavy- and light-chain variable domains and form the antigen-binding site, the relative orientation of the heavy- and light-chain variable domains can create another source of structural diversity leading to changes in antigen binding. Here, we first reexamine the diversity of V(L):V(H) orientations in existing antibody crystal structures using 153 nonredundant sequences, demonstrating that the variation in V(L):V(H) orientation is greater than that expected from effects of crystal packing, antigen binding, or the presence of antibody constant regions and increases, on average, as sequence similarity decreases for residues in the interface between the domains. We developed a tool for predicting the relative orientations of the heavy- and light-chain variable domains using side-chain rotamer sampling in the interface and molecular-mechanics-based energy calculations. When using variable domain backbones from the crystal structures, the predicted orientation is very close (<1 A RMSD) to the crystallographically observed orientation in most cases, confirming that the V(L):V(H) orientation is determined by the antibody sequence and suggesting an approach to predicting the relative orientation of the variable domains when building homology models of antibodies. When applied to antibody homology models generated from templates with 55-75% sequence identity, we predict the V(L):V(H) orientation of 20 antibodies with an average/median RMSD of 2.1/1.6 A to the crystal structures.