Structure-Based Assignment of Ile, Leu, and Val Methyl Groups in the Active and Inactive Forms of the Mitogen-Activated Protein Kinase Extracellular Signal-Regulated Kinase 2

Conformation selection by ATP-competitive inhibitors and allosteric communication in ERK2

Activation of the extracellular signal-regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named 'L' and 'R,' where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here, we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

Conformation Selection by ATP-competitive Inhibitors and Allosteric Communication in ERK2

Activation of the extracellular signal regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named "L" and "R", where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

Assignment of Ala, Ile, LeuproS, Met, and ValproS methyl groups of the protruding domain of murine norovirus capsid protein VP1 using methyl-methyl NOEs, site directed mutagenesis, and pseudocontact shifts

The protruding domain (P-domain) of the murine norovirus (MNV) capsid protein VP1 is essential for infection. It mediates receptor binding and attachment of neutralizing antibodies. Protein NMR studies into interactions of the P-domain with ligands will yield insights not easily available from other biophysical techniques and will extend our understanding of MNV attachment to host cells. Such studies require at least partial NMR assignments. Here, we describe the assignment of about 70% of the Ala, Ile, Leu^proS, Met, and Val^proS methyl groups. An unfavorable distribution of methyl group resonance signals prevents complete assignment based exclusively on 4D HMQC-NOESY-HMQC experiments, yielding assignment of only 55 out of 100 methyl groups. Therefore, we created point mutants and measured pseudo contact shifts, extending and validating assignments based on methyl-methyl NOEs. Of note, the P-domains are present in two different forms caused by an approximate equal distribution of trans- and cis-configured proline residues in position 361.

Backbone resonance assignments of the catalytic and regulatory domains of Ca2+/calmodulin-dependent protein kinase 1D

The CaMK subfamily of Ser/Thr kinases are regulated by calmodulin interactions with their C-terminal regions. They are exemplified by Ca²⁺/calmodulin dependent protein kinase 1δ which is known as CaMK1D, CaMKIδ or CKLiK. CaMK1D mediates intracellular signalling downstream of Ca²⁺ influx and thereby exhibits amplifications of Ca²⁺signals and polymorphisms that have been implicated in breast cancer and diabetes. Here we report the backbone ¹H, ¹³C, ¹⁵N assignments of the 38 kDa human CaMK1D protein in its free state, including both the canonical bi-lobed kinase fold as well as the autoinhibitory and calmodulin binding domains.

Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2

The activation loop segment in protein kinases is a common site for regulatory phosphorylation. In extracellular signal-regulated kinase 2 (ERK2), dual phosphorylation and conformational rearrangement of the activation loop accompany enzyme activation. X-ray structures show the active conformation to be stabilized by multiple ion pair interactions between phosphorylated threonine and tyrosine residues in the loop and six arginine residues in the kinase core. Despite the extensive salt bridge network, nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiments show that the phosphorylated activation loop is conformationally mobile on a microsecond to millisecond time scale. The dynamics of the loop match those of previously reported global exchange within the kinase core region and surrounding the catalytic site that have been found to facilitate productive nucleotide binding. Mutations in the core region that alter these global motions also alter the dynamics of the activation loop. Conversely, mutations in the activation loop perturb the global exchange within the kinase core. Together, these findings provide evidence for coupling between motions in the activation loop and those surrounding the catalytic site in the active state of the kinase. Thus, the activation loop segment in dual-phosphorylated ERK2 is not held statically in the active X-ray conformation but instead undergoes exchange between conformers separated by a small energetic barrier, serving as part of a dynamic allosteric network controlling nucleotide binding and catalytic function.

Automated assignment of methyl NMR spectra from large proteins

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.

Complete assignment of Ala, Ile, LeuProS, Met and ValProS methyl groups of the protruding domain from human norovirus GII.4 Saga

Attachment of human noroviruses to histo blood group antigens (HBGAs) is thought to be essential for infection, although how this binding event promotes infection is unknown. Recent studies have shown that 60% of all GII.4 epidemic strains may undergo a spontaneous post-translational modification (PTM) in an amino acid located adjacent to the binding pocket for HBGAs. This transformation proceeds with an estimated half-life of 1-2 days under physiological conditions, dramatically affecting HBGA recognition. The surface-exposed position of this PTM and its sequence conservation suggests a relevant role in immune escape and host-cell recognition. As a first step towards the understanding of the biological implications of this PTM at atomic resolution, we report the complete assignment of methyl resonances of a MIL^ProSV^ProSA methyl-labeled sample of a 72 kDa protruding domain from a GII.4 Saga human norovirus strain. Assignments were obtained from methyl-methyl NOESY experiments combined with site-directed mutagenesis and automated assignment. This data provides the basis for a detailed characterization of the PTM-driven modulation of immune recognition in human norovirus on a molecular level.

Intra-residue methyl-methyl correlations for valine and leucine residues in large proteins from a 3D-HMBC-HMQC experiment

Methyl-TROSY based NMR experiments have over the last two decades become one of the most important means to characterise dynamics and functional mechanisms of large proteins and macromolecular machines in solution. The chemical shift assignment of methyl groups in large proteins is, however, still not trivial and it is typically performed using backbone-dependent experiments in a 'divide and conquer' approach, mutations, structure-based assignments or a combination of these. Structure-based assignment of methyl groups is an emerging strategy, which reduces the time and cost required as well as providing a method that is independent of a backbone assignment. One crucial step in available structure-based assignment protocols is linking the two prochiral methyl groups of leucine and valine residues. This has previously been achieved by recording NOESY spectra with short mixing times or by comparing NOESY spectra. Herein, we present a method based on through-bond scalar coupling transfers, a 3D-HMBC-HMQC experiment, to link the intra-residue methyl groups of leucine and valine. It is shown that the HMBC-HMQC method has several advantages over solely using NOESY spectra since a unique intra-residue cross-peak is observed. Moreover, overlap in the methyl-TROSY HMQC spectrum can easily be identified with the HMBC-HMQC experiment, thereby removing possible ambiguities in the assignment.

Modulating multi-functional ERK complexes by covalent targeting of a recruitment site in vivo

Recently, the targeting of ERK with ATP-competitive inhibitors has emerged as a potential clinical strategy to overcome acquired resistance to BRAF and MEK inhibitor combination therapies. In this study, we investigate an alternative strategy of targeting the D-recruitment site (DRS) of ERK. The DRS is a conserved region that lies distal to the active site and mediates ERK-protein interactions. We demonstrate that the small molecule BI-78D3 binds to the DRS of ERK2 and forms a covalent adduct with a conserved cysteine residue (C159) within the pocket and disrupts signaling in vivo. BI-78D3 does not covalently modify p38MAPK, JNK or ERK5. BI-78D3 promotes apoptosis in BRAF inhibitor-naive and resistant melanoma cells containing a BRAF V600E mutation. These studies provide the basis for designing modulators of protein-protein interactions involving ERK, with the potential to impact ERK signaling dynamics and to induce cell cycle arrest and apoptosis in ERK-dependent cancers.

Automatic structure-based NMR methyl resonance assignment in large proteins

Isotopically labeled methyl groups provide NMR probes in large, otherwise deuterated proteins. However, the resonance assignment constitutes a bottleneck for broader applicability of methyl-based NMR. Here, we present the automated MethylFLYA method for the assignment of methyl groups that is based on methyl-methyl nuclear Overhauser effect spectroscopy (NOESY) peak lists. MethylFLYA is applied to five proteins (28–358 kDa) comprising a total of 708 isotope-labeled methyl groups, of which 612 contribute NOESY cross peaks. MethylFLYA confidently assigns 488 methyl groups, i.e. 80% of those with NOESY data. Of these, 459 agree with the reference, 6 were different, and 23 were without reference assignment. MethylFLYA assigns significantly more methyl groups than alternative algorithms, has an average error rate of 1%, modest runtimes of 0.4–1.2 h, and can handle arbitrary isotope labeling patterns and data from other types of NMR spectra.

The structures and dynamics of large proteins can be studied with methyl-based NMR but peak assignment is still challenging. Here the authors present MethylFLYA that allows automated assignment of methyl groups and apply it to five proteins with molecular weights in the range from 28 to 358 kDa.

Complete assignment of Ala, Ile, Leu, Met and Val methyl groups of human blood group A and B glycosyltransferases using lanthanide-induced pseudocontact shifts and methyl-methyl NOESY

Human blood group A and B glycosyltransferases (GTA, GTB) are highly homologous glycosyltransferases. A number of high-resolution crystal structures is available showing that these enzymes convert from an open conformation into a catalytically active closed conformation upon substrate binding. However, the mechanism of glycosyltransfer is still under debate, and the precise nature as well as the time scales of conformational transitions are unknown. NMR offers a variety of experiments to shine more light on these unresolved questions. Therefore, in a first step we have assigned all methyl resonance signals in MILVA labeled samples of GTA and GTB, still a challenging task for 70 kDa homodimeric proteins. Assignments were obtained from methyl-methyl NOESY experiments, and from measurements of lanthanide-induced pseudocontact shifts (PCS) using high resolution crystal structures as templates. PCSs and chemical shift perturbations, induced by substrate analogue binding, suggest that the fully closed state is not adopted in the presence of lanthanide ions.

Mapping Antibody Epitopes by Solution NMR Spectroscopy: Practical Considerations

Identifying an epitope, the region of the antigen in contact with an antibody, is useful in both basic and pharmaceutical research, as well as in vaccine design. Solution NMR spectroscopy is particularly well suited to the residue level characterization of intermolecular interfaces, including antibody-antigen interactions, and thus to epitope identification. Here, we describe the use of NMR for residue level characterization of protein epitopes, focusing on experimental protocols and practical considerations, highlighting advantages and drawbacks of the approach.

Sparse multidimensional iterative lineshape-enhanced (SMILE) reconstruction of both non-uniformly sampled and conventional NMR data

Implementation of a new algorithm, SMILE, is described for reconstruction of non-uniformly sampled two-, three- and four-dimensional NMR data, which takes advantage of the known phases of the NMR spectrum and the exponential decay of underlying time domain signals. The method is very robust with respect to the chosen sampling protocol and, in its default mode, also extends the truncated time domain signals by a modest amount of non-sampled zeros. SMILE can likewise be used to extend conventional uniformly sampled data, as an effective multidimensional alternative to linear prediction. The program is provided as a plug-in to the widely used NMRPipe software suite, and can be used with default parameters for mainstream application, or with user control over the iterative process to possibly further improve reconstruction quality and to lower the demand on computational resources. For large data sets, the method is robust and demonstrated for sparsities down to ca 1%, and final all-real spectral sizes as large as 300 Gb. Comparison between fully sampled, conventionally processed spectra and randomly selected NUS subsets of this data shows that the reconstruction quality approaches the theoretical limit in terms of peak position fidelity and intensity. SMILE essentially removes the noise-like appearance associated with the point-spread function of signals that are a default of five-fold above the noise level, but impacts the actual thermal noise in the NMR spectra only minimally. Therefore, the appearance and interpretation of SMILE-reconstructed spectra is very similar to that of fully sampled spectra generated by Fourier transformation.

NMR assignments of sparsely labeled proteins using a genetic algorithm

Sparse isotopic labeling of proteins for NMR studies using single types of amino acid (15N or 13C enriched) has several advantages. Resolution is enhanced by reducing numbers of resonances for large proteins, and isotopic labeling becomes economically feasible for glycoproteins that must be expressed in mammalian cells. However, without access to the traditional triple resonance strategies that require uniform isotopic labeling, NMR assignment of crosspeaks in heteronuclear single quantum coherence (HSQC) spectra is challenging. We present an alternative strategy which combines readily accessible NMR data with known protein domain structures. Based on the structures, chemical shifts are predicted, NOE cross-peak lists are generated, and residual dipolar couplings (RDCs) are calculated for each labeled site. Simulated data are then compared to measured values for a trial set of assignments and scored. A genetic algorithm uses the scores to search for an optimal pairing of HSQC crosspeaks with labeled sites. While none of the individual data types can give a definitive assignment for a particular site, their combination can in most cases. Four test proteins previously assigned using triple resonance methods and a sparsely labeled glycosylated protein, Robo1, previously assigned by manual analysis, are used to validate the method and develop a criterion for identifying sites assigned with high confidence.

Solution NMR Spectroscopy Provides an Avenue for the Study of Functionally Dynamic Molecular Machines: The Example of Protein Disaggregation

Solution-based NMR spectroscopy has been an important tool for studying the structure and dynamics of relatively small proteins and protein complexes with aggregate molecular masses under approximately 50 kDa. The development of new experiments and labeling schemes, coupled with continued improvements in hardware, has significantly reduced this size limitation, enabling atomic-resolution studies of molecular machines in the 1 MDa range. In this Perspective, some of the important advances are highlighted in the context of studies of molecular chaperones involved in protein disaggregation. New insights into the structural biology of disaggregation obtained from NMR studies are described, focusing on the unique capabilities of the methodology for obtaining atomic-resolution descriptions of dynamic systems.