Michael P. Latham

NMR Methods for Studying the Structure and Dynamics of RNA

Proper functioning of RNAs requires the formation of complex three‐dimensional structures combined with the ability to rapidly interconvert between multiple functional states. This review covers recent advances in isotope‐labeling strategies and NMR experimental approaches that have promise for facilitating solution structure determinations and dynamics studies of biologically active RNAs. Improved methods for the production of isotopically labeled RNAs combined with new multidimensional heteronuclear NMR experiments make it possible to dramatically reduce spectral crowding and simplify resonance assignments for RNAs. Several novel applications of experiments that directly detect hydrogen‐bonding interactions are discussed. These studies demonstrate how NMR spectroscopy can be used to distinguish between possible secondary structures and identify mechanisms of ligand binding in RNAs. A variety of recently developed methods for measuring base and sugar residual dipolar couplings are described. NMR residual dipolar coupling techniques provide valuable data for determining the long‐range structure and orientation of helical regions in RNAs. A number of studies are also presented where residual dipolar coupling constraints are used to determine the global structure and dynamics of RNAs. NMR relaxation data can be used to probe the dynamics of macromolecules in solution. The power dependence of transverse rotating‐frame relaxation rates was used here to study dynamics in the minimal hammerhead ribozyme. Improved methods for isotopically labeling RNAs combined with new types of structural data obtained from a growing repertoire of NMR experiments are facilitating structural and dynamic studies of larger RNAs.

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.

Phosphorylation releases constraints to domain motion in ERK2

Significance This paper uses NMR methods to compare the dynamics of a protein kinase in its active and inactive states. The results show that domain movements in the MAP kinase ERK2 are inherently constrained until the enzyme is activated by phosphorylation, with the constraint located at the hinge region. This represents an important mode for dynamical regulation in ERK2, not anticipated from previous X-ray structural analyses. Protein motions control enzyme catalysis through mechanisms that are incompletely understood. Here NMR 13C 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.

Is buffer a good proxy for a crowded cell-like environment? A comparative NMR study of calmodulin side-chain dynamics in buffer and E. coli lysate.

Biophysical studies of protein structure and dynamics are typically performed in a highly controlled manner involving only the protein(s) of interest. Comparatively fewer such studies have been carried out in the context of a cellular environment that typically involves many biomolecules, ions and metabolites. Recently, solution NMR spectroscopy, focusing primarily on backbone amide groups as reporters, has emerged as a powerful technique for investigating protein structure and dynamics in vivo and in crowded “cell-like” environments. Here we extend these studies through a comparative analysis of Ile, Leu, Val and Met methyl side-chain motions in apo, Ca2+-bound and Ca2+, peptide-bound calmodulin dissolved in aqueous buffer or in E. coli lysate. Deuterium spin relaxation experiments, sensitive to pico- to nano-second time-scale processes and Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, reporting on millisecond dynamics, have been recorded. Both similarities and differences in motional properties are noted for calmodulin dissolved in buffer or in lysate. These results emphasize that while significant insights can be obtained through detailed “test-tube” studies, experiments performed under conditions that are “cell-like” are critical for obtaining a comprehensive understanding of protein motion in vivo and therefore for elucidating the relation between motion and function.

Probing non-specific interactions of Ca2+-calmodulin in E. coli lysate

The biological environment in which a protein performs its function is a crowded milieu containing millions of molecules that can potentially lead to a great many transient, non-specific interactions. NMR spectroscopy is especially well suited to study these weak molecular contacts. Here, non-specific interactions between the Ca2+-bound form of calmodulin (CaM) and non-cognate proteins in Escherichia coli lysate are explored using Ile, Leu, Val and Met methyl probes. Changes in CaM methyl chemical shifts as a function of added E. coli lysate are measured to determine a minimum ‘average’ dissociation constant for interactions between Ca2+-CaM and E. coli lysate proteins. 2H R2 and 13C R1 spin relaxation rates report on the binding reaction as well. Our results further highlight the power of methyl containing side-chains for characterizing biomolecular interactions, even in complex in-cell like environments.

Understanding the mechanism of proteasome 20S core particle gating

Significance The proteasome plays a critical role in regulating cellular homeostasis by degrading target protein substrates. Gating residues that are part of the barrel-like proteasome structure provide a barrier to prevent inadvertent proteolysis. With NMR spectroscopy, the mechanism by which the gates interconvert between open and closed states is studied. We show that collisions with water molecules provide a driving force for gate interconversion and that the process takes place via many small steps that involve protein segments smaller than 4 Å. The effect of cellular lysate on proteasome gating is explored, establishing that gating equilibria and kinetics are not perturbed relative to buffer solutions under conditions of similar viscosity. The 20S core particle proteasome is a molecular machine playing an important role in cellular function by degrading protein substrates that no longer are required or that have become damaged. Regulation of proteasome activity occurs, in part, through a gating mechanism controlling the sizes of pores at the top and bottom ends of the symmetric proteasome barrel and restricting access to catalytic sites sequestered in the lumen of the structure. Although atomic resolution models of both open and closed states of the proteasome have been elucidated, the mechanism by which gates exchange between these states remains to be understood. Here, this is investigated by using magnetization transfer NMR spectroscopy focusing on the 20S proteasome core particle from Thermoplasma acidophilum. We show from viscosity-dependent proteasome gating kinetics that frictional forces originating from random solvent motions are critical for driving the gating process. Notably, a small effective hydrodynamic radius (EHR; <4Å) is obtained, providing a picture in which gate exchange proceeds through many steps involving only very small segment sizes. A small EHR further suggests that the kinetics of gate interconversion will not be affected appreciably by large viscogens, such as macromolecules found in the cell, so long as they are inert. Indeed, measurements in cell lysate reveal that the gate interconversion rate decreases only slightly, demonstrating that controlled studies in vitro provide an excellent starting point for understanding regulation of 20S core particle function in complex, biologically relevant environments.

Viscosity-dependent kinetics of protein conformational exchange: microviscosity effects and the need for a small viscogen.

Conformational rearrangements are critical to a variety of biological processes including protein folding and misfolding, ligand binding, enzyme catalysis, and signal transduction. Viscosity-dependent kinetics measurements can provide crucial insights into the dynamics of protein conformational exchange by highlighting the relative importance of frictional forces derived from either solvent or from internal protein interactions in activating the exchange reaction. Here, we analyze the kinetics of interconversion between the native and intermediate states of the four helix bundle FF domain recorded in solutions containing the viscogens glycerol or bovine serum albumin (BSA), using the viscosity measured from the translational diffusion of probes of different sizes. In the large viscogen BSA, we demonstrate that vastly different internal friction values are obtained using the different viscosity measures, leading to conflicting interpretations of the role of solvent friction in the interconversion. We show that this can be a consequence of the small effective hydrodynamic radius of the protein conformational transition and differences between solution micro- and macroscopic viscosities that are germane in this case. In general, correct values of internal friction can only be obtained by carrying out measurements using small viscogens.

A similar in vitro and in cell lysate folding intermediate for the FF domain.

Understanding the mechanisms by which proteins fold into their three-dimensional structures, including a description of the intermediates that are formed during the folding process, remains a goal of protein science. Most studies are performed under carefully controlled conditions in which the folding reaction is monitored in a buffer solution that is far from the natural milieu of the cell. Here, we have used (13)C and (1)H relaxation dispersion NMR spectroscopy to study folding of the FF domain in both Escherichia coli and Saccharomyces cerevisiae cellular lysates. We find that a conformationally excited state is populated in both lysates, which is very similar in structure to a folding intermediate observed in previous studies in buffer, with the kinetics and thermodynamics of the interconversion between native and intermediate conformers somewhat changed. The results point to the importance of extending folding studies beyond the test tube yet emphasize that insights can be obtained through careful experiments recorded in controlled buffer solutions.

High-resolution pyrimidine- and ribose-specific 4D HCCH-COSY spectra of RNA using the filter diagonalization method

The NMR spectra of nucleic acids suffer from severe peak overlap, which complicates resonance assignments. 4D NMR experiments can overcome much of the degeneracy in 2D and 3D spectra; however, the linear increase in acquisition time with each new dimension makes it impractical to acquire high-resolution 4D spectra using standard Fourier transform (FT) techniques. The filter diagonalization method (FDM) is a numerically efficient algorithm that fits the entire multi-dimensional time-domain data to a set of multi-dimensional oscillators. Selective 4D constant-time HCCH-COSY experiments that correlate the H5–C5–C6–H6 base spin systems of pyrimidines or the H1′–C1′–C2′–H2′ spin systems of ribose sugars were acquired on the 13C-labeled iron responsive element (IRE) RNA. FDM-processing of these 4D experiments recorded with only 8 complex points in the indirect dimensions showed superior spectral resolution than FT-processed spectra. Practical aspects of obtaining optimal FDM-processed spectra are discussed. The results here demonstrate that FDM-processing can be used to obtain high-resolution 4D spectra on a medium sized RNA in a fraction of the acquisition time normally required for high-resolution, high-dimensional spectra.