Phosphorylation releases constraints to domain motion in ERK2.

Protein motions control enzyme catalysis through mechanisms that are incompletely understood. Here NMR (13)C relaxation dispersion experiments were used to monitor changes in side-chain motions that occur in response to activation by phosphorylation of the MAP kinase ERK2. NMR data for the methyl side chains on Ile, Leu, and Val residues showed changes in conformational exchange dynamics in the microsecond-to-millisecond time regime between the different activity states of ERK2. In inactive, unphosphorylated ERK2, localized conformational exchange was observed among methyl side chains, with little evidence for coupling between residues. Upon dual phosphorylation by MAP kinase kinase 1, the dynamics of assigned methyls in ERK2 were altered throughout the conserved kinase core, including many residues in the catalytic pocket. The majority of residues in active ERK2 fit to a single conformational exchange process, with kex ≈ 300 s(-1) (kAB ≈ 240 s(-1)/kBA ≈ 60 s(-1)) and pA/pB ≈ 20%/80%, suggesting global domain motions involving interconversion between two states. A mutant of ERK2, engineered to enhance conformational mobility at the hinge region linking the N- and C-terminal domains, also induced two-state conformational exchange throughout the kinase core, with exchange properties of kex ≈ 500 s(-1) (kAB ≈ 15 s(-1)/kBA ≈ 485 s(-1)) and pA/pB ≈ 97%/3%. Thus, phosphorylation and activation of ERK2 lead to a dramatic shift in conformational exchange dynamics, likely through release of constraints at the hinge.

Dynamics of protein kinases: insights from nuclear magnetic resonance.

Protein kinases are ubiquitous enzymes with critical roles in cellular processes and pathology. As a result, researchers have studied their activity and regulatory mechanisms extensively. Thousands of X-ray structures give snapshots of the architectures of protein kinases in various states of activation and ligand binding. However, the extent of and manner by which protein motions and conformational dynamics underlie the function and regulation of these important enzymes is not well understood. Nuclear magnetic resonance (NMR) methods provide complementary information about protein conformation and dynamics in solution. However, until recently, the large size of these enzymes prevented researchers from using these methods with kinases. Developments in transverse relaxation-optimized spectroscopy (TROSY)-based techniques and more efficient isotope labeling strategies are now allowing researchers to carry out NMR studies on full-length protein kinases. In this Account, we describe recent insights into the role of dynamics in protein kinase regulation and catalysis that have been gained from NMR measurements of chemical shift changes and line broadening, residual dipolar couplings, and relaxation. These findings show strong associations between protein motion and events that control kinase activity. Dynamic and conformational changes occurring at ligand binding sites and other regulatory domains of these proteins propagate to conserved kinase core regions that mediate catalytic function. NMR measurements of slow time scale (microsecond to millisecond) motions also reveal that kinases carry out global exchange processes that synchronize multiple residues and allosteric interconversion between conformational states. Activating covalent modifications or ligand binding to form the Michaelis complex can induce these global processes. Inhibitors can also exploit the exchange properties of kinases by using conformational selection to form dynamically quenched states. These investigations have revealed that kinases are highly dynamic enzymes, whose regulation by interdomain interactions, ligand binding, and covalent modifications involve changes in motion and conformational equilibrium in a manner that can be correlated with function. Thus, NMR provides a unique window into the role of protein dynamics in kinase regulation and catalysis with important implications for drug design.

Slow inhibition and conformation selective properties of extracellular signal-regulated kinase 1 and 2 inhibitors.

The mitogen-activated protein (MAP) kinase pathway is a target for anticancer therapy, validated using inhibitors of B-Raf and MAP kinase kinase (MKK) 1 and 2. Clinical outcomes show a high frequency of acquired resistance in patient tumors, involving upregulation of activity of the MAP kinase, extracellular signal-regulated kinase (ERK) 1 and 2. Thus, inhibitors for ERK1/2 are potentially important for targeted therapeutics against cancer. The structures and potencies of different ERK inhibitors have been published, but their kinetic mechanisms have not been characterized. Here we perform enzyme kinetic studies on six representative ERK inhibitors, with potencies varying from 100 pM to 20 µM. Compounds with significant biological activity (IC50 < 100 nM) that inhibit in the subnanomolar range (Vertex-11e and SCH772984) display slow-onset inhibition and represent the first inhibitors of ERK2 known to demonstrate slow dissociation rate constants (values of 0.2 and 1.1 h(-1), respectively). Furthermore, we demonstrate using kinetic competition assays that Vertex-11e binds with differing affinities to ERK2 in its inactive, unphosphorylated and active, phosphorylated forms. Finally, two-dimensional heteronuclear multiple-quantum correlation nuclear magnetic resonance experiments reveal that distinct conformational states are formed in complexes of Vertex-11e with inactive and active ERK2. Importantly, two conformers interconvert in equilibrium in the active ERK2 apoenzyme, but Vertex-11e strongly shifts the equilibrium completely to one conformer. Thus, a high-affinity, slow dissociation inhibitor stabilizes different enzyme conformations depending on the activity state of ERK2 and reveals properties of conformational selection toward the active kinase.

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.

Resonance assignments are the first step in most NMR studies of protein structure, function, and dynamics. Standard protein assignment methods employ through-bond backbone experiments on uniformly (13)C/(15)N-labeled proteins. For larger proteins, this through-bond assignment procedure often breaks down due to rapid relaxation and spectral overlap. The challenges involved in studies of larger proteins led to efficient methods for (13)C labeling of side chain methyl groups, which have favorable relaxation properties and high signal-to-noise. These methyls are often still assigned by linking them to the previously assigned backbone, thus limiting the applications for larger proteins. Here, a structure-based procedure is described for assignment of (13)C(1)H3-labeled methyls by comparing distance information obtained from three-dimensional methyl-methyl nuclear Overhauser effect (NOE) spectroscopy with the X-ray structure. The Ile, Leu, or Val (ILV) methyl type is determined by through-bond experiments, and the methyl-methyl NOE data are analyzed in combination with the known structure. A hierarchical approach was employed that maps the largest observed "NOE-methyl cluster" onto the structure. The combination of identification of ILV methyl type with mapping of the NOE-methyl clusters greatly simplifies the assignment process. This method was applied to the inactive and active forms of the 42-kDa ILV (13)C(1)H3-methyl labeled extracellular signal-regulated kinase 2 (ERK2), leading to assignment of 60% of the methyls, including 90% of Ile residues. A series of ILV to Ala mutants were analyzed, which helped confirm the assignments. These assignments were used to probe the local and long-range effects of ligand binding to inactive and active ERK2.

Optimized reverse micelle surfactant system for high-resolution NMR spectroscopy of encapsulated proteins and nucleic acids dissolved in low viscosity fluids.

An optimized reverse micelle surfactant system has been developed for solution nuclear magnetic resonance studies of encapsulated proteins and nucleic acids dissolved in low viscosity fluids. Comprising the nonionic 1-decanoyl-rac-glycerol and the zwitterionic lauryldimethylamine-N-oxide (10MAG/LDAO), this mixture is shown to efficiently encapsulate a diverse set of proteins and nucleic acids. Chemical shift analyses of these systems show that high structural fidelity is achieved upon encapsulation. The 10MAG/LDAO surfactant system reduces the molecular reorientation time for encapsulated macromolecules larger than ~20 kDa leading to improved overall NMR performance. The 10MAG/LDAO system can also be used for solution NMR studies of lipid-modified proteins. New and efficient strategies for optimization of encapsulation conditions are described. 10MAG/LDAO performs well in both the low viscosity pentane and ultralow viscosity liquid ethane and therefore will serve as a general surfactant system for initiating solution NMR studies of proteins and nucleic acids.

Efficient detection of hydrogen bonds in dynamic regions of RNA by sensitivity-optimized NMR pulse sequences.

Improved Sensitivity: Efficient NMR experiments are presented for determining the secondary structure in large and dynamic RNAs using J-couplings across hydrogen bonds. The experiments provide up to eight-fold improved sensitivity and thus enable detection of base pairs in dynamic regions even in large RNAs.

The cavity-chaperone Skp protects its substrate from aggregation but allows independent folding of substrate domains.

Outer membrane proteins (OMPs) of gram-negative bacteria are synthesized in the cytosol and must cross the periplasm before insertion into the outer membrane. The 17-kDa protein (Skp) is a periplasmic chaperone that assists the folding and insertion of many OMPs, including OmpA, a model OMP with a membrane embedded beta-barrel domain and a periplasmic alphabeta domain. Structurally, Skp belongs to a family of cavity-containing chaperones that bind their substrates in the cavity, protecting them from aggregation. However, some substrates, such as OmpA, exceed the capacity of the chaperone cavity, posing a mechanistic challenge. Here, we provide direct NMR evidence that, while bound to Skp, the beta-barrel domain of OmpA is maintained in an unfolded state, whereas the periplasmic domain is folded in its native conformation. Complementary cross-linking and NMR relaxation experiments show that the OmpA beta-barrel is bound deep within the Skp cavity, whereas the folded periplasmic domain protrudes outside of the cavity where it tumbles independently from the rest of the complex. This domain-based chaperoning mechanism allows the transport of beta-barrels across the periplasm in an unfolded state, which may be important for efficient insertion into the outer membrane.

Ultrafast dynamics show that the theophylline and 3-methylxanthine aptamers employ a conformational capture mechanism for binding their ligands.

RNAs often exhibit a high degree of conformational dynamics and heterogeneity, leading to a rugged energy landscape. However, the roles of conformational heterogeneity and rapid dynamics in molecular recognition or RNA function have not been extensively elucidated. Ultrafast time-resolved fluorescence spectroscopic experiments were used here to probe picosecond dynamics of the theophylline-binding RNA aptamer. These studies showed that multiple conformations are populated in the free RNA, indicating that this aptamer employs a conformational capture mechanism for ligand binding. The base on residue 27 in an internal loop exists in at least three conformational states in the free RNA, including binding competent and incompetent states that have distinct fluorescence decay signatures indicating different base stacking interactions. Picosecond dynamics were also detected by anisotropy experiments, where these motions indicate additional dynamics for base 27. The picosecond data show that theophylline binding shifts the equilibrium for conformations of base 27 from primarily stacked in the free RNA to mostly unstacked in the RNA-theophylline complex, as observed in the previous NMR structure. In contrast, base 10 in a second internal loop is mostly preorganized in the free RNA, consistent with it being stacked between G11 and G25, as is observed in the bound state. Picosecond dynamics were also measured on a modified aptamer that binds with higher affinity to 3-methylxanthine than theophylline. The modified aptamer shows less heterogeneity in the aptamer-3-methylxanthine complex than what is observed in the theophylline aptamer-theophylline complex.

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme.

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca(2+), Mg(2+), Mn(2+), and Sr(2+) have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.

NMR chemical exchange as a probe for ligand-binding kinetics in a theophylline-binding RNA aptamer.

The apparent on and off rate constants for binding of theophylline to its RNA aptamer in the absence of Mg(2+) were determined here by 2D (1)H-(1)H ZZ-exchange NMR spectroscopy. Analysis of the buildup rate of the exchange cross peaks for several base-paired imino protons in the RNA yielded an apparent k(on) of 600 M(-1) s(-1). This small apparent k(on) results because the free RNA exist as a dynamic equilibrium of inactive states rapidly interconverting with a low population of active species. The data found here indicate that the RNA aptamer employs a conformational selection mechanism for binding theophylline in the absence of Mg(2+). The kinetic data found here also explain a very unusual property of this RNA-theophylline system: slow exchange on the NMR chemical shift time scale for a weakly binding complex. To our knowledge, it is unprecedented to have such a weakly binding complex (K(d) approximately 3.0 mM at 15 degrees C) show slow exchange on the NMR chemical shift time scale, but the results clearly demonstrate that slow exchange and weak binding are readily rationalized by a small k(on). Comparisons with other ligand-receptor interactions are presented.

Chemical shift anisotropy of imino 15N nuclei in Watson-Crick base pairs from magic angle spinning liquid crystal NMR and nuclear spin relaxation.

Knowledge of (15)N chemical shift anisotropy is prerequisite both for quantitative interpretation of nuclear spin relaxation rates in terms of local dynamics and for the use of residual chemical shift anisotropy (RCSA) as a constraint in structure determination. Accurate measurement of the very small RCSA from the difference in (15)N chemical shift under isotropic and weakly aligning liquid crystalline conditions is very sensitive to minute differences in sample conditions, such as pH or ionic strength. For this reason, chemical shifts were measured for the same solution, under static liquid crystalline alignment, and under magic angle spinning conditions where alignment relative to the magnetic field is removed. Measurements were made for 14 well-resolved G-N(1) and 6 U-N(3) (15)N nuclei in a sample of tRNA(Val). Fitting these RCSA data together with (15)N-(1)H dipole-CSA cross-correlated relaxation measurements to the recently refined structural model of tRNA(Val) yields the magnitude, asymmetry, and orientation of the (15)N CSA tensors.

Longitudinal-relaxation-enhanced NMR experiments for the study of nucleic acids in solution.

Atomic-resolution information on the structure and dynamics of nucleic acids is essential for a better understanding of the mechanistic basis of many cellular processes. NMR spectroscopy is a powerful method for studying the structure and dynamics of nucleic acids; however, solution NMR studies are currently limited to relatively small nucleic acids at high concentrations. Thus, technological and methodological improvements that increase the experimental sensitivity and spectral resolution of NMR spectroscopy are required for studies of larger nucleic acids or protein-nucleic acid complexes. Here we introduce a series of imino-proton-detected NMR experiments that yield an over 2-fold increase in sensitivity compared to conventional pulse schemes. These methods can be applied to the detection of base pair interactions, RNA-ligand titration experiments, measurement of residual dipolar (15)N-(1)H couplings, and direct measurements of conformational transitions. These NMR experiments employ longitudinal spin relaxation enhancement techniques that have proven useful in protein NMR spectroscopy. The performance of these new experiments is demonstrated for a 10 kDa TAR-TAR*(GA) RNA kissing complex and a 26 kDa tRNA.

Measurement of imino 1H-1H residual dipolar couplings in RNA.

Imino (1)H-(15)N residual dipolar couplings (RDCs) provide additional structural information that complements standard (1)H-(1)H NOEs leading to improvements in both the local and global structure of RNAs. Here, we report measurement of imino (1)H-(1)H RDCs for the Iron Responsive Element (IRE) RNA and native E. coli tRNA(Val) using a BEST-Jcomp-HMQC2 experiment. (1)H-(1)H RDCs are observed between the imino protons in G-U wobble base pairs and between imino protons on neighboring base pairs in both RNAs. These imino (1)H-(1)H RDCs complement standard (1)H-(15)N RDCs because the (1)H-(1)H vectors generally point along the helical axis, roughly perpendicular to (1)H-(15)N RDCs. The use of longitudinal relaxation enhancement increased the signal-to-noise of the spectra by ~3.5-fold over the standard experiment. The ability to measure imino (1)H-(1)H RDCs offers a new restraint, which can be used in NMR domain orientation and structural studies of RNAs.

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme.

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH3, HPO4(2-) or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2.G23 base pair was the most stable base pair in all the duplexes. The conserved U-1*G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2'-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1*G22, C-2*G23 and A2*U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.

Imino proton exchange rates imply an induced-fit binding mechanism for the VEGF165-targeting aptamer, Macugen.

The 2'-fluoro/2'-O-methyl modified RNA aptamer Macugen is a potent inhibitor of the angiogenic regulatory protein, VEGF165. Macugen binds with high affinity to the heparin-binding domain (HBD) of VEGF165. Hydrogen exchange rates of the imino protons were measured for free Macugen and Macugen bound to the HBD or full-length VEGF to better understand the mechanism for high affinity binding. The results here show that the internal loop and hairpin loop of Macugen are highly dynamic in the free state and are greatly stabilized and/or protected from solvent upon protein binding.

Flexibility in the Periplasmic Domain of BamA Is Important for Function.

The ß-barrel assembly machine (BAM) mediates the biogenesis of outer membrane proteins (OMPs) in Gram-negative bacteria. BamA, the central BAM subunit composed of a transmembrane ß-barrel domain linked to five polypeptide transport-associated (POTRA) periplasmic domains, is thought to bind nascent OMPs and undergo conformational cycling to catalyze OMP folding and insertion. One model is that conformational flexibility between POTRA domains is part of this conformational cycling. Nuclear magnetic resonance (NMR) spectroscopy was used here to study the flexibility of the POTRA domains 1-5 in solution. NMR relaxation studies defined effective rotational correlational times and together with residual dipolar coupling data showed that POTRA1-2 is flexibly linked to POTRA3-5. Mutants of BamA that restrict flexibility between POTRA2 and POTRA3 by disulfide crosslinking displayed impaired function in vivo. Together these data strongly support a model in which conformational cycling of hinge motions between POTRA2 and POTRA3 in BamA is required for biological function.

Refined solution structure of the iron-responsive element RNA using residual dipolar couplings.

The iron-responsive element (IRE) is a 30nt RNA motif located in the non-coding regions of mRNAs of proteins involved in iron regulation. In humans, the IRE plays a direct role in the control of iron levels by post-transcriptional regulation of the ferritin and transferrin receptor proteins through highly specific recognition by IRE-binding proteins. The IRE fold is representative of many RNA motifs that contain helical domains separated by a bulge or internal loop. The global structures of such extended multi-domain RNAs are not well defined by conventional NMR-distance and torsion angle structural restraints. Residual dipolar couplings (RDCs) are employed here to better define the global structure of the IRE RNA in solution. RDCs contain valuable long-range structural information that compliments the short-range structural data derived from standard NOE-distance and torsion angle restraints. Several approaches for estimating alignment tensor parameters and incorporating RDCs into RNA structure determinations are compared. Both the local and global structure of the IRE are improved significantly by refinement with RDCs. These RDC refinements provide insight on the conformational dynamics of the IRE. These studies highlight some issues that need to be addressed when incorporating RDCs in solution structure determinations of nucleic acids. The approach used here should prove valuable for structure determinations of various multi-domain systems.

Structure of the BamC two-domain protein obtained by Rosetta with a limited NMR data set.

The CS-RDC-NOE Rosetta program was used to generate the solution structure of a 27-kDa fragment of the Escherichia coli BamC protein from a limited set of NMR data. The BamC protein is a component of the essential five-protein ß-barrel assembly machine in E. coli. The first 100 residues in BamC were disordered in solution. The Rosetta calculations showed that BamC101₋344 forms two well-defined domains connected by an ~18-residue linker, where the relative orientation of the domains was not defined. Both domains adopt a helix-grip fold previously observed in the Bet v 1 superfamily. ¹5N relaxation data indicated a high degree of conformational flexibility for the linker connecting the N-terminal domain and the C-terminal domain in BamC. The results here show that CS-RDC-NOE Rosetta is robust and has a high tolerance for misassigned nuclear Overhauser effect restraints, greatly simplifying NMR structure determinations.

Structure and flexibility of the complete periplasmic domain of BamA: the protein insertion machine of the outer membrane.

Folding and insertion of ß-barrel outer membrane proteins (OMPs) is essential for Gram-negative bacteria. This process is mediated by the multiprotein complex BAM, composed of the essential ß-barrel OMP BamA and four lipoproteins (BamBCDE). The periplasmic domain of BamA is key for its function and contains five "polypeptide transport-associated" (POTRA) repeats. Here, we report the crystal structure of the POTRA4-5 tandem, containing the essential for BAM complex formation and cell viability POTRA5. The domain orientation observed in the crystal is validated by solution NMR and SAXS. Using previously determined structures of BamA POTRA1-4, we present a spliced model of the entire BamA periplasmic domain validated by SAXS. Solution scattering shows that conformational flexibility between POTRA2 and 3 gives rise to compact and extended conformations. The length of BamA in its extended conformation suggests that the protein may bridge the inner and outer membranes across the periplasmic space.

Solution structure of tRNAVal from refinement of homology model against residual dipolar coupling and SAXS data.

A procedure is presented for refinement of a homology model of E. coli tRNA(Val), originally based on the X-ray structure of yeast tRNA(Phe), using experimental residual dipolar coupling (RDC) and small angle X-ray scattering (SAXS) data. A spherical sampling algorithm is described for refinement against SAXS data that does not require a globbic approximation, which is particularly important for nucleic acids where such approximations are less appropriate. Substantially higher speed of the algorithm also makes its application favorable for proteins. In addition to the SAXS data, the structure refinement employed a sparse set of NMR data consisting of 24 imino N-H(N) RDCs measured with Pf1 phage alignment, and 20 imino N-H(N) RDCs obtained from magnetic field dependent alignment of tRNA(Val). The refinement strategy aims to largely retain the local geometry of the 58% identical tRNA(Phe) by ensuring that the atomic coordinates for short, overlapping segments of the ribose-phosphate backbone and the conserved base pairs remain close to those of the starting model. Local coordinate restraints are enforced using the non-crystallographic symmetry (NCS) term in the XPLOR-NIH or CNS software package, while still permitting modest movements of adjacent segments. The RDCs mainly drive the relative orientation of the helical arms, whereas the SAXS restraints ensure an overall molecular shape compatible with experimental scattering data. The resulting structure exhibits good cross-validation statistics (jack-knifed Q (free) = 14% for the Pf1 RDCs, compared to 25% for the starting model) and exhibits a larger angle between the two helical arms than observed in the X-ray structure of tRNA(Phe), in agreement with previous NMR-based tRNA(Val) models.