Nuclear Magnetic Resonance Studies of Biopolymer Dynamics

Hydrodynamic modelling of protein conformation in solution: ELLIPS and HYDRO

The last three decades has seen some important advances in our ability to represent the conformation of proteins in solution on the basis of hydrodynamic measurements. Advances in theoretical modeling capabilities have been matched by commensurate advances in the precision of hydrodynamic measurements. We consider the advances in whole-body (simple ellipsoid-based) modeling—still useful for providing an overall idea of molecular shape, particularly for those systems where only a limited amount of data is available—and outline the ELLIPS suite of algorithms which facilitates the use of this approach. We then focus on bead modeling strategies, particularly the surface or shell–bead approaches and the HYDRO suite of algorithms. We demonstrate how these are providing great insights into complex issues such as the conformation of immunoglobulins and other multi-domain complexes.

An R(1ρ) expression for a spin in chemical exchange between two sites with unequal transverse relaxation rates

An analytical expression is derived for the rotating frame relaxation rate, R(1ρ), of a spin exchanging between two sites with different transverse relaxation times. A number of limiting cases are examined, with the equation reducing to formulae derived previously under the assumption of equivalent relaxation rates at each site. The measurement of a pair off-resonance R(1ρ) values, with the carrier displaced equally on either side of the observed correlation, forms the basis of one of the approaches for obtaining signs of chemical shift differences, Δω, of exchanging nuclei. The results presented here establish that this method is relatively insensitive to differential transverse relaxation rates between the exchaning states, greatly simplifying the calculation of optimal parameters in R(1ρ) based experiments that are used for measurement of signs of Δω.

Site-resolved measurement of microsecond-to-millisecond conformational-exchange processes in proteins by solid-state NMR spectroscopy

We demonstrate that conformational exchange processes in proteins on microsecond-to-millisecond time scales can be detected and quantified by solid-state NMR spectroscopy. We show two independent approaches that measure the effect of conformational exchange on transverse relaxation parameters, namely Carr-Purcell-Meiboom-Gill relaxation-dispersion experiments and measurement of differential multiple-quantum coherence decay. Long coherence lifetimes, as required for these experiments, are achieved by the use of highly deuterated samples and fast magic-angle spinning. The usefulness of the approaches is demonstrated by application to microcrystalline ubiquitin. We detect a conformational exchange process in a region of the protein for which dynamics have also been observed in solution. Interestingly, quantitative analysis of the data reveals that the exchange process is more than 1 order of magnitude slower than in solution, and this points to the impact of the crystalline environment on free energy barriers.

Protein function and allostery: a dynamic relationship

Allostery is a fundamental process by which distant sites within a protein system sense each other. Allosteric regulation is such an efficient mechanism that it is used to control protein activity in most biological processes, including signal transduction, metabolism, catalysis, and gene regulation. Over recent years, our view and understanding of the fundamental principles underlying allostery have been enriched and often utterly reshaped. This has been especially so for powerful techniques such as nuclear magnetic resonance spectroscopy, which offers an atomic view of the intrinsic motions of proteins. Here, I discuss recent results on the catabolite activator protein (CAP) that have drastically revised our view about how allosteric interactions are modulated. CAP has provided the first experimentally identified system showing that (i) allostery can be mediated through changes in protein motions, in the absence of changes in the mean structure of the protein, and (ii) favorable changes in protein motions may activate allosteric proteins that are otherwise structurally inactive.

NMR insights into protein allostery

Allosterism is one of nature's principal methods for regulating protein function. Allosterism utilizes ligand binding at one site to regulate the function of the protein by modulating the structure and dynamics of a distant binding site. In this review, we first survey solution NMR techniques and how they may be applied to the study of allostery. Subsequently, we describe several examples of application of NMR to protein allostery and highlight the unique insight provided by this experimental technique.

Progress in correlation spectroscopy at ultra-fast magic-angle spinning: basic building blocks and complex experiments for the study of protein structure and dynamics

Recent progress in multi-dimensional solid-state NMR correlation spectroscopy at high static magnetic fields and ultra-fast magic-angle spinning is discussed. A focus of the review is on applications to protein resonance assignment and structure determination as well as on the characterization of protein dynamics in the solid state. First, the consequences of ultra-fast spinning on sensitivity and sample heating are considered. Recoupling and decoupling techniques at ultra-fast MAS are then presented, as well as more complex experiments assembled from these basic building blocks. Furthermore, we discuss new avenues in biomolecular solid-state NMR spectroscopy that become feasible in the ultra-fast spinning regime, such as sensitivity enhancement based on paramagnetic doping, and the prospect of direct proton detection.

Examining the relationship between RNA function and motion using nuclear magnetic resonance

The biological function of proteins and nucleic acids relies on their complex structures, yet dynamics provides an additional layer of functional adaptability. Numerous studies have demonstrated that RNA is only able to perform the multitude of functions for which it is responsible by readily changing its conformation in response to binding of proteins or small molecules. Examination of RNA dynamics is therefore essential to understanding its biological function. Nuclear magnetic resonance (NMR) has emerged as a leading technique for the examination of RNA motion and conformational transitions. It can examine domain motions as well as motion with atomic level resolution over a wide range of time scales. This review examines how NMR spectroscopy can be applied to examine the relationship between function and dynamics in RNA.

Ergodicity and efficiency of cross-polarization in NMR of static solids

Cross-polarization transfer is employed in virtually every solid-state NMR experiment to enhance magnetization of low-gamma spins. Theory and experiment is used to assess the magnitude of the final quasistationary magnetization amplitude. The many-body density matrix equation is solved for relatively large (up to N=14) spin systems without the spin-temperature assumption for the final spin states. Simulations show that about 13% of the thermodynamic limit is still retained within the proton bath. To test this theoretical prediction, a combination of a reverse cross-polarization experiment and multiple contacts is employed to show that the thermodynamic limit of magnetization cannot be transferred from high- to low-gamma nuclei in a single contact. Multiple contacts, however, fully transfer the maximum magnetization. A simple diffusion on a cone model shows that slow dynamics can affect the build up profile for the transferred magnetization.

Quantification of protein backbone hydrogen-deuterium exchange rates by solid state NMR spectroscopy

We present the quantification of backbone amide hydrogen-deuterium exchange rates (HDX) for immobilized proteins. The experiments make use of the deuterium isotope effect on the amide nitrogen chemical shift, as well as on proton dilution by deuteration. We find that backbone amides in the microcrystalline α-spectrin SH3 domain exchange rather slowly with the solvent (with exchange rates negligible within the individual (15)N-T (1) timescales). We observed chemical exchange for 6 residues with HDX exchange rates in the range from 0.2 to 5 s(-1). Backbone amide (15)N longitudinal relaxation times that we determined previously are not significantly affected for most residues, yielding no systematic artifacts upon quantification of backbone dynamics (Chevelkov et al. 2008b). Significant exchange was observed for the backbone amides of R21, S36 and K60, as well as for the sidechain amides of N38, N35 and for W41ε. These residues could not be fit in our previous motional analysis, demonstrating that amide proton chemical exchange needs to be considered in the analysis of protein dynamics in the solid-state, in case D(2)O is employed as a solvent for sample preparation. Due to the intrinsically long (15)N relaxation times in the solid-state, the approach proposed here can expand the range of accessible HDX rates in the intermediate regime that is not accessible so far with exchange quench and MEXICO type experiments.

Nonexponential decay of internal rotational correlation functions of native proteins and self-similar structural fluctuations

Structural fluctuations of a protein are essential for the function of native proteins and for protein folding. To understand how the main chain in the native state of a protein fluctuates on different time scales, we examined the rotational correlation functions (RCFs), C(t), of the backbone N-H bonds and of the dihedral angles γ built on four consecutive C(α) atoms. Using molecular dynamics simulations of a model α/β protein (VA3) in its native state, we demonstrate that these RCFs decay as stretched exponentials, ln[C(t)] ≈ D(α)t(α) with a constant D(α) and an exponent α (0 < α < 0.35) varying with the free-energy profiles (FEPs) along the amino acid sequence. The probability distributions of the fluctuations of the main chain computed at short time scale (1 ps) were identical to those computed at large time scale (1 ns) if the time is rescaled by a factor depending on α < 1. This self-similar property and the nonexponential decays (α ≠ 1) of the RCFs are described by a rotational diffusion equation with a time-dependent diffusion coefficient D(t) = αD(α)t(α-1). The present findings agree with observations of subdiffusion (α < 1) of fluorescent probes within a protein molecule. The subdiffusion of (15)N-H bonds did not affect the value of the order parameter S(2) extracted from the NMR relaxation data by assuming normal diffusion (α = 1) of (15)N-H bonds on a nanosecond time scale. However, we found that the RCF does not converge to S(2) on the nanosecond time scale for residues with multiple-minima FEPs.

Interpretation of biomolecular NMR spin relaxation parameters

Biomolecular nuclear magnetic resonance (NMR) spin relaxation experiments provide exquisite information on the picosecond to nanosecond timescale motions of bond vectors. Spin-lattice (T1) and spin-spin (T2) relaxation times and the steady-state nuclear Overhauser effect (NOE) are the first set of parameters extracted from typical 15N or 13C NMR relaxation experiments. Therefore, verifying that T1, T2, and NOE are consistent with theoretical predictions is an important step before carrying out the more detailed model-free and reduced spectral density mapping analyses commonly employed. In this mini-review, we discuss the essential motional parameters used to describe biomolecular dynamics in the context of a variety of examples of folded and intrinsically disordered proteins and peptides in aqueous and membrane mimetic environments. Estimates of these parameters can be used as input for an online interface, introduced herein, allowing plotting of trends of T1, T2, and NOE with magnetic field strength. The plots may serve as a first-check to the spectroscopist preparing to embark on a detailed NMR relaxation analysis.

Analyzing protein folding cooperativity by differential scanning calorimetry and NMR spectroscopy

Some marginally stable proteins undergo microsecond time scale folding reactions that involve significant populations of partly ordered forms, making it difficult to discern individual steps in their folding pathways. It has been suggested that many of these proteins fold non-cooperatively, with no significant barriers to separate the energy landscape into distinct thermodynamic states. Here we present an approach for studying the cooperativity of rapid protein folding with a combination of differential scanning calorimetry (DSC), nuclear magnetic resonance (NMR) relaxation dispersion experiments, and an analysis of the temperature dependence of amide (1)H and (15)N chemical shifts. We applied this method to the PBX homeodomain (PBX-HD), which folds on the microsecond time scale and produces a broad DSC thermogram with an elevated and steeply sloping native-state heat capacity baseline, making it a candidate for barrierless folding. However, by globally fitting the NMR thermal melt and DSC data, and by comparing these results to those obtained from the NMR relaxation dispersion experiments, we show that the native form of the protein undergoes two-state exchange with a small population of the thermally denatured form, well below the melting temperature. This result directly demonstrates the coexistence of distinct folded and unfolded forms and firmly establishes that folding of PBX-HD is cooperative. Further, we see evidence of large-scale structural and dynamical changes within the native state by NMR, which helps to explain the broad and shallow DSC profile. This study illustrates the potential of combining calorimetry with NMR dynamics experiments to dissect mechanisms of protein folding.

Measurement of signs of chemical shift differences between ground and excited protein states: a comparison between H(S/M)QC and R1rho methods

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR spectroscopy has emerged as a powerful tool for quantifying the kinetics and thermodynamics of millisecond exchange processes between a major, populated ground state and one or more minor, low populated and often invisible 'excited' conformers. Analysis of CPMG data-sets also provides the magnitudes of the chemical shift difference(s) between exchanging states (|Deltavarpi|), that inform on the structural properties of the excited state(s). The sign of Deltavarpi is, however, not available from CPMG data. Here we present one-dimensional NMR experiments for measuring the signs of (1)H(N) and (13)C(alpha) Deltavarpi values using weak off-resonance R (1rho ) relaxation measurements, extending the spin-lock approach beyond previous applications focusing on the signs of (15)N and (1)H(alpha) shift differences. The accuracy of the method is established by using an exchanging system where the invisible, excited state can be converted to the visible, ground state by altering conditions so that the signs of Deltavarpi values obtained from the spin-lock approach can be validated with those measured directly. Further, the spin-lock experiments are compared with the established H(S/M)QC approach for measuring the signs of chemical shift differences. For the Abp1p and Fyn SH3 domains considered here it is found that while H(S/M)QC measurements provide signs for more residues than the spin-lock data, the two different methodologies are complementary, so that combining both approaches frequently produces signs for more residues than when the H(S/M)QC method is used alone.

Conformational flexibility of Y145Stop human prion protein amyloid fibrils probed by solid-state nuclear magnetic resonance spectroscopy

Amyloid aggregates of a C-truncated Y145Stop mutant of human prion protein, huPrP23-144, associated with a heritable amyloid angiopathy, have previously been shown to contain a compact, relatively rigid, and beta-sheet-rich approximately 30-residue amyloid core near the C-terminus under physiologically relevant conditions. In contrast, the remaining huPrP23-144 residues display considerable conformational dynamics, as evidenced by the absence of corresponding signals in cross-polarization (CP)-based solid-state NMR (SSNMR) spectra under ambient conditions and their emergence in analogous spectra recorded at low temperature on frozen fibril samples. Here, we present the direct observation of residues comprising the flexible N-terminal domain of huPrP23-144 amyloid by using 2D J-coupling-based magic-angle spinning (MAS) SSNMR techniques. Chemical shifts for these residues indicate that the N-terminal domain is effectively an ensemble of protein chains with random-coil-like conformations. Interestingly, a detailed analysis of signal intensities in CP-based 3D SSNMR spectra suggests that non-negligible molecular motions may also be occurring on the NMR time scale within the relatively rigid core of huPrP23-144 amyloid. To further investigate this hypothesis, quantitative measurements of backbone dipolar order parameters and transverse spin relaxation rates were performed for the core residues. The observed order parameters indicate that, on the submicrosecond time scale, these residues are effectively rigid and experience only highly restricted and relatively uniform motions similar to those characteristic for well-structured regions of microcrystalline proteins. On the other hand, significant variations in magnitude of transverse spin relaxation rates were noted for residues present at different locations within the core region and correlated with observed differences in spectral intensities. While interpreted only qualitatively at the present time, the extent of the observed variations in transverse relaxation rates is consistent with the presence of relatively slow, microsecond-millisecond time scale chemical exchange type phenomena within the huPrP23-144 amyloid core.

Backbone dynamics of the antifungal Psd1 pea defensin and its correlation with membrane interaction by NMR spectroscopy

Plant defensins are cysteine-rich cationic peptides, components of the innate immune system. The antifungal sensitivity of certain exemplars was correlated to the level of complex glycosphingolipids in the membrane of fungi strains. Psd1 is a 46 amino acid residue defensin isolated from pea seeds which exhibit antifungal activity. Its structure is characterized by the so-called cysteine-stabilized alpha/beta motif linked by three loops as determined by two-dimensional NMR. In the present work we explored the measurement of heteronuclear Nuclear Overhauser Effects, R1 and R2 (15)N relaxation ratios, and chemical shift to probe the backbone dynamics of Psd1 and its interaction with membrane mimetic systems with phosphatidylcholine (PC) or dodecylphosphocholine (DPC) with glucosylceramide (CMH) isolated from Fusarium solani. The calculated R2 values predicted a slow motion around the highly conserved among Gly12 residue and also in the region of the Turn3 His36-Trp38. The results showed that Psd1 interacts with vesicles of PC or PC:CMH in slightly different forms. The interaction was monitored by chemical shift perturbation and relaxation properties. Using this approach we could map the loops as the binding site of Psd1 with the membrane. The major binding epitope showed conformation exchange properties in the mus-ms timescale supporting the conformation selection as the binding mechanism. Moreover, the peptide corresponding to part of Loop1 (pepLoop1: Gly12 to Ser19) is also able to interact with DPC micelles acquiring a stable structure and in the presence of DPC:CMH the peptide changes to an extended conformation, exhibiting NOE mainly with the carbohydrate and ceramide parts of CMH.

Intermediate rate atomic trajectories of RNA by solid-state NMR spectroscopy

Many RNAs undergo large conformational changes in response to the binding of proteins and small molecules. However, when RNA functional dynamics occur in the nanosecond-microsecond time scale, they become invisible to traditional solution NMR relaxation methods. Residual dipolar coupling methods have revealed the presence of extensive nanosecond-microsecond domain motions in HIV-1 TAR RNA, but this technique lacks information on the rates of motions. We have used solid-state deuterium NMR to quantitatively describe trajectories of key residues in TAR by exploiting the sensitivity of this technique to motions that occur in the nanosecond-microsecond regime. Deuterium line shape and relaxation data were used to model motions of residues within the TAR binding interface. The resulting motional models indicate two functionally essential bases within the single-stranded bulge sample both the free and Tat-bound conformations on the microsecond time scale in the complete absence of the protein. Thus, our results strongly support a conformational capture mechanism for recognition: the protein does not induce a new RNA structure, but instead captures an already-populated conformation.

NMR spectroscopy brings invisible protein states into focus

Molecular dynamics are essential for protein function. In some cases these dynamics involve the interconversion between ground state, highly populated conformers and less populated higher energy structures ('excited states') that play critical roles in biochemical processes. Here we describe recent advances in NMR spectroscopy methods that enable studies of these otherwise invisible excited states at an atomic level and that help elucidate their important relation to function. We discuss a range of examples from molecular recognition, ligand binding, enzyme catalysis and protein folding that illustrate the role that motion plays in 'funneling' conformers along preferred pathways that facilitate their biological function.

Dynamics of reassembled thioredoxin studied by magic angle spinning NMR: snapshots from different time scales

Solid-state NMR spectroscopy can be used to probe internal protein dynamics in the absence of the overall molecular tumbling. In this study, we report (15)N backbone dynamics in differentially enriched 1-73(U-(13)C,(15)N)/74-108(U-(15)N) reassembled thioredoxin on multiple time scales using a series of 2D and 3D MAS NMR experiments probing the backbone amide (15)N longitudinal relaxation, (1)H-(15)N dipolar order parameters, (15)N chemical shift anisotropy (CSA), and signal intensities in the temperature-dependent and (1)H T(2)'-filtered NCA experiments. The spin-lattice relaxation rates R(1) (R(1) = 1/T(1)) were observed in the range from 0.012 to 0.64 s(-1), indicating large site-to-site variations in dynamics on pico- to nanosecond time scales. The (1)H-(15)N dipolar order parameters, <S>, and (15)N CSA anisotropies, delta(sigma), reveal the backbone mobilities in reassembled thioredoxin, as reflected in the average <S> = 0.89 +/- 0.06 and delta(sigma) = 92.3 +/- 5.2 ppm, respectively. From the aggregate of experimental data from different dynamics methods, some degree of correlation between the motions on the different time scales has been suggested. Analysis of the dynamics parameters derived from these solid-state NMR experiments indicates higher mobilities for the residues constituting irregular secondary structure elements than for those located in the alpha-helices and beta-sheets, with no apparent systematic differences in dynamics between the alpha-helical and beta-sheet residues. Remarkably, the dipolar order parameters derived from the solid-state NMR measurements and the corresponding solution NMR generalized order parameters display similar qualitative trends as a function of the residue number. The comparison of the solid-state dynamics parameters to the crystallographic B-factors has identified the contribution of static disorder to the B-factors. The combination of longitudinal relaxation, dipolar order parameter, and CSA line shape analyses employed in this study provides snapshots of dynamics and a new insight on the correlation of these motions on multiple time scales.

A solid-state deuterium NMR and sum-frequency generation study of the side-chain dynamics of peptides adsorbed onto surfaces

The artificial amphiphilic peptide LKalpha14 adopts a helical structure at interfaces, with opposite orientation of its leucine (L, hydrophobic) and lysine (K, hydrophilic) side chains. When peptides are adsorbed onto surfaces, different residue side chains necessarily have different proximities to the surface, depending on both their position in the helix and the composition of the surface itself. Deuterating the individual leucine residues (isopropyl-d(7)) permits the use of solid-state deuterium NMR spectroscopy as a site-specific probe of side-chain dynamics. In conjunction with sum-frequency generation as a probe of the peptide-binding face, we demonstrate that the mobility of specific leucine side chains at the interface is quantifiable in terms of their surface proximity.

Structure and interdomain dynamics of apoptosis-associated speck-like protein containing a CARD (ASC)

The human protein ASC is a key mediator in apoptosis and inflammation. Through its two death domains (pyrin and CARD) ASC interacts with cell death executioners, acts as an essential adapter for inflammasome integrity, and oligomerizes into functional supramolecular assemblies. However, these functions are not understood at the structural-dynamic level. This study reports the solution structure and interdomain dynamics of full-length ASC. The pyrin and CARD domains are structurally independent six-helix bundle motifs connected by a 23-residue linker. The CARD structure reveals two distinctive characteristics; helix 1 is not fragmented as in all other known CARDs, and its electrostatic surface shows a uniform distribution of positive and negative charges, whereas these are commonly separated into two areas in other death domains. The linker adopts residual structure resulting in a back-to-back orientation of the domains, which avoids steric interference of each domain with the binding site of the other. NMR relaxation experiments show that the linker is flexible despite the residual structure. This flexibility could help expand the relative volume occupied by each domain, thus increasing the capture radius for effectors. Based on the ASC structure, a tentative model is proposed to illustrate how ASC oligomerizes via CARD and pyrin homophilic interactions. Moreover, ASC oligomers have been analyzed by atomic force microscopy, showing a predominant species of disk-like particles of approximately 12-nm diameter and approximately 1-nm height. Taken together, these results provide structural insight into the behavior of ASC as an adapter molecule.

Effects of amantadine on the dynamics of membrane-bound influenza A M2 transmembrane peptide studied by NMR relaxation

The molecular motions of membrane proteins in liquid-crystalline lipid bilayers lie at the interface between motions in isotropic liquids and in solids. Specifically, membrane proteins can undergo whole-body uniaxial diffusion on the microsecond time scale. In this work, we investigate the (1)H rotating-frame spin-lattice relaxation (T (1rho)) caused by the uniaxial diffusion of the influenza A M2 transmembrane peptide (M2TMP), which forms a tetrameric proton channel in lipid bilayers. This uniaxial diffusion was proved before by (2)H, (15)N and (13)C NMR lineshapes of M2TMP in DLPC bilayers. When bound to an inhibitor, amantadine, the protein exhibits significantly narrower linewidths at physiological temperature. We now investigate the origin of this line narrowing through temperature-dependent (1)H T (1rho) relaxation times in the absence and presence of amantadine. Analysis of the temperature dependence indicates that amantadine decreases the correlation time of motion from 2.8 +/- 0.9 mus for the apo peptide to 0.89 +/- 0.41 micros for the bound peptide at 313 K. Thus the line narrowing of the bound peptide is due to better avoidance of the NMR time scale and suppression of intermediate time scale broadening. The faster diffusion of the bound peptide is due to the higher attempt rate of motion, suggesting that amantadine creates better-packed and more cohesive helical bundles. Analysis of the temperature dependence of ln(T_1rho(-1)) indicates that the activation energy of motion increased from 14.0 +/- 4.0 kJ/mol for the apo peptide to 23.3 +/- 6.2 kJ/mol for the bound peptide. This higher activation energy indicates that excess amantadine outside the protein channel in the lipid bilayer increases the membrane viscosity. Thus, the protein-bound amantadine speeds up the diffusion of the helical bundles while the excess amantadine in the bilayer increases the membrane viscosity.

Quantitative analysis of backbone motion in proteins using MAS solid-state NMR spectroscopy

We present a comprehensive analysis of protein dynamics for a micro-crystallin protein in the solid-state. Experimental data include (15)N T (1) relaxation times measured at two different magnetic fields as well as (1)H-(15)N dipole, (15)N CSA cross correlated relaxation rates which are sensitive to the spectral density function J(0) and are thus a measure of T (2) in the solid-state. In addition, global order parameters are included from a (1)H,(15)N dipolar recoupling experiment. The data are analyzed within the framework of the extended model-free Clore-Lipari-Szabo theory. We find slow motional correlation times in the range of 5 and 150 ns. Assuming a wobbling in a cone motion, the amplitude of motion of the respective amide moiety is on the order of 10 degrees for the half-opening angle of the cone in most of the cases. The experiments are demonstrated using a perdeuterated sample of the chicken alpha-spectrin SH3 domain.

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms-ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in (2)H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state (13)C relaxation measurements, we establish that ns-micros motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.

Domain features of the peripheral stalk subunit H of the methanogenic A1AO ATP synthase and the NMR solution structure of H(1-47)

A series of truncated forms of subunit H were generated to establish the domain features of that protein. Circular dichroism analysis demonstrated that H is divided at least into a C-terminal coiled-coil domain within residues 54-104, and an N-terminal domain formed by adjacent alpha-helices. With a cysteine at the C-terminus of each of the truncated proteins (H(1-47), H(1-54), H(1-59), H(1-61), H(1-67), H(1-69), H(1-71), H(1-78), H(1-80), H(1-91), and H(47-105)), the residues involved in formation of the coiled-coil interface were determined. Proteins H(1-54), H(1-61), H(1-69), and H(1-80) showed strong cross-link formation, which was weaker in H(1-47), H(1-59), H(1-71), and H(1-91). A shift in disulfide formation between cysteines at positions 71 and 80 reflected an interruption in the periodicity of hydrophobic residues in the region 71AEKILEETEKE81. To understand how the N-terminal domain of H is formed, we determined for the first time, to our knowledge, the solution NMR structure of H(1-47), which revealed an alpha-helix between residues 15-42 and a flexible N-terminal stretch. The alpha-helix includes a kink that would bring the two helices of the C-terminus into the coiled-coil arrangement. H(1-47) revealed a strip of alanines involved in dimerization, which were tested by exchange to single cysteines in subunit H mutants.

Accurate measurement of alpha proton chemical shifts of excited protein states by relaxation dispersion NMR spectroscopy

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy can provide detailed information about low populated, invisible states of protein molecules, including backbone chemical shifts of the invisible conformer and bond vector orientations that can be used as structural constraints. Notably, the measurement of 1Halpha chemical shifts in excited protein states has not been possible to date because, in the absence of suitable labeling, the homonuclear proton scalar coupling network in side chains of proteins leads to a significant degradation in the performance of proton-based relaxation dispersion experiments. Here we have overcome this problem through a labeling scheme in which proteins are prepared with U-2H glucose and 50% D2O/50% H2O that results in deuteration levels of between 50-88% at the Cbeta carbon. Effects from residual 1Halpha-1Hbeta scalar couplings can be suppressed through a new NMR experiment that is presented here. The utility of the methodology is demonstrated on a ligand binding exchanging system and it is shown that 1Halpha chemical shifts extracted from dispersion profiles are, on average, accurate to 0.03 ppm, an order of magnitude better than they can be predicted from structure using a database approach. The ability to measure 1Halpha chemical shifts of invisible conformers is particularly important because such shifts are sensitive to both secondary and tertiary structure. Thus, the methodology presented is a valuable addition to a growing list of experiments for characterizing excited protein states that are difficult to study using the traditional techniques of structural biology.

How binding of small molecule and peptide ligands to HIV-1 TAR alters the RNA motional landscape

The HIV-1 TAR RNA represents a well-known paradigm to study the role of dynamics and conformational change in RNA function. This regulatory RNA changes conformation in response to binding of Tat protein and of a variety of peptidic and small molecule ligands, indicating that its conformational flexibility and intrinsic dynamics play important roles in molecular recognition. We have used ¹³C NMR relaxation experiments to examine changes in the motional landscape of HIV-1 TAR in the presence of three ligands of different affinity and specificity. The ligands are argininamide, a linear peptide mimic of the Tat basic domain and a cyclic peptide that potently inhibits Tat-dependent activation of transcription. All three molecules induce the same motional characteristics within the three nucleotides bulge that represents the Tat-binding site. However, the cyclic peptide has a unique motional signature in the apical loop, which represents a binding site for the essential host co-factor cyclin T1. These results suggest that all peptidic mimics of Tat induce the same dynamics in TAR within this protein binding site. However, the new cyclic peptide mimic of Tat represents a new class of ligands with a unique effect on the dynamics and the structure of the apical loop.

Monitoring molecular interactions by NMR

The ability of proteins to interact with small molecules or other proteins is essential in all aspects of biology. In many cases these interactions cause detectable changes in NMR chemical shifts, lineshapes, and relaxation rates and therefore provide a means by which to study these biologically important phenomena. Here we review the theory upon which this analysis is based, provide several illustrative examples, and highlight potential problems in the study of binding interactions by solution NMR.

Furanose dynamics in the HhaI methyltransferase target DNA studied by solution and solid-state NMR relaxation

Both solid-state and solution NMR relaxation measurements are routinely used to quantify the internal dynamics of biomolecules, but in very few cases have these two techniques been applied to the same system, and even fewer attempts have been made so far to describe the results obtained through these two methods through a common theoretical framework. We have previously collected both solution 13C and solid-state 2H relaxation measurements for multiple nuclei within the furanose rings of several nucleotides of the DNA sequence recognized by HhaI methyltransferase. The data demonstrated that the furanose rings within the GCGC recognition sequence are very flexible, with the furanose rings of the cytidine, which is the methylation target, experiencing the most extensive motions. To interpret these experimental results quantitatively, we have developed a dynamic model of furanose rings based on the analysis of solid-state 2H line shapes. The motions are modeled by treating bond reorientations as Brownian excursions within a restoring potential. By applying this model, we are able to reproduce the rates of 2H spin-lattice relaxation in the solid and 13C spin-lattice relaxation in solution using comparable restoring force constants and internal diffusion coefficients. As expected, the 13C relaxation rates in solution are less sensitive to motions that are slower than overall molecular tumbling than to the details of global molecular reorientation, but are somewhat more sensitive to motions in the immediate region of the Larmor frequency. Thus, we conclude that the local internal motions of this DNA oligomer in solution and in the hydrated solid state are virtually the same, and we validate an approach to the conjoint analysis of solution and solid-state NMR relaxation and line shapes data, with wide applicability to many biophysical problems.

Contact replacement for NMR resonance assignment

MOTIVATION

Complementing its traditional role in structural studies of proteins, nuclear magnetic resonance (NMR) spectroscopy is playing an increasingly important role in functional studies. NMR dynamics experiments characterize motions involved in target recognition, ligand binding, etc., while NMR chemical shift perturbation experiments identify and localize protein-protein and protein-ligand interactions. The key bottleneck in these studies is to determine the backbone resonance assignment, which allows spectral peaks to be mapped to specific atoms. This article develops a novel approach to address that bottleneck, exploiting an available X-ray structure or homology model to assign the entire backbone from a set of relatively fast and cheap NMR experiments.

RESULTS

We formulate contact replacement for resonance assignment as the problem of computing correspondences between a contact graph representing the structure and an NMR graph representing the data; the NMR graph is a significantly corrupted, ambiguous version of the contact graph. We first show that by combining connectivity and amino acid type information, and exploiting the random structure of the noise, one can provably determine unique correspondences in polynomial time with high probability, even in the presence of significant noise (a constant number of noisy edges per vertex). We then detail an efficient randomized algorithm and show that, over a variety of experimental and synthetic datasets, it is robust to typical levels of structural variation (1-2 AA), noise (250-600%) and missings (10-40%). Our algorithm achieves very good overall assignment accuracy, above 80% in alpha-helices, 70% in beta-sheets and 60% in loop regions.

AVAILABILITY

Our contact replacement algorithm is implemented in platform-independent Python code. The software can be freely obtained for academic use by request from the authors.

Conformational dynamics of an intact virus: order parameters for the coat protein of Pf1 bacteriophage

This study has examined the atomic-level dynamics of the protein in the capsid of filamentous phage Pf1. This capsid consists of approximately 7,300 small subunits of only 46 aa in a helical array around a highly extended, circular single-stranded DNA molecule of 7,349 nt. Measurements were made of site-specific, solid-state NMR order parameters, S, the values which are dimensionless quantities between 0 (mobile) and 1 (static) that characterize the amplitudes of molecular bond angular motions that are faster than microseconds. It was found that the protein subunit backbone is very static, and of particular interest, it appears to be static at residues glycine 15 and glutamine 16 where it had been previously thought to be mobile. In contrast to the backbone, several side chains display large-amplitude angular motions. Side chains on the virion exterior that interact with solvent are highly mobile, but surprisingly, the side chains of residues arginine 44 and lysine 45 near the DNA deep in the interior of the virion are also highly mobile. The large-amplitude dynamic motion of these positively charged side chains in their interactions with the DNA were not previously expected. The results reveal a highly dynamic aspect of a DNA-protein interface within a virus.

Structural mobility of the monomeric C-terminal domain of the HIV-1 capsid protein

The capsid protein of HIV-1 (p24) (CA) forms the mature capsid of the human immunodeficiency virus. Capsid assembly involves hexamerization of the N-terminal domain and dimerization of the C-terminal domain of CA (CAC), and both domains constitute potential targets for anti-HIV therapy. CAC homodimerization occurs mainly through its second helix, and it is abolished when its sole tryptophan is mutated to alanine. This mutant, CACW40A, resembles a transient monomeric intermediate formed during dimerization. Its tertiary structure is similar to that of the subunits in the dimeric, non-mutated CAC, but the segment corresponding to the second helix samples different conformations. The present study comprises a comprehensive examination of the CACW40A internal dynamics. The results obtained, with movements sampling a wide time regime (from pico- to milliseconds), demonstrate the high flexibility of the whole monomeric protein. The conformational exchange phenomena on the micro-to-millisecond time scale suggest a role for internal motions in the monomer-monomer interactions and, thus, flexibility of the polypeptide chain is likely to contribute to the ability of the protein to adopt different conformational states, depending on the biological environment.

Measurement of bond vector orientations in invisible excited states of proteins

The focus of structural biology is on studies of the highly populated, ground states of biological molecules; states that are only sparsely and transiently populated are more difficult to probe because they are invisible to most structural methods. Yet, such states can play critical roles in biochemical processes such as ligand binding, enzyme catalysis, and protein folding. A description of these states in terms of structure and dynamics is, therefore, of great importance. Here, we present a method, based on relaxation dispersion NMR spectroscopy of weakly aligned molecules in a magnetic field, that can provide such a description by direct measurement of backbone amide bond vector orientations in transient, low populated states that are not observable directly. Such information, obtained through the measurement of residual dipolar couplings, has until now been restricted to proteins that produce observable spectra. The methodology is applied and validated in a study of the binding of a target peptide to an SH3 domain from the yeast protein Abp1p and subsequently used in an application to protein folding of a mutational variant of the Fyn SH3 domain where (1)H-(15)N dipolar couplings of the invisible unfolded state of the domain are obtained. The approach, which can be used to obtain orientational restraints at other sites in proteins as well, promises to significantly extend the available information necessary for providing a site-specific characterization of structural properties of transient, low populated states that have to this point remained recalcitrant to detailed analysis.

Solution NMR of supramolecular complexes: providing new insights into function

Solution NMR spectroscopy is an extremely powerful technology for the study of biomolecular dynamics and site-specific molecular interactions. An important limitation in the past has been molecule size, with molecular weights of targets seldom exceeding 50 kDa. New labeling technology and NMR experiments are changing this paradigm so that applications for investigating supramolecular complexes are starting to become feasible. Here we describe a strategy developed in our laboratory that involves the use of labeled methyl groups of isoleucine, leucine and valine residues in proteins as probes, along with experiments that significantly enhance the lifetimes of the resulting signals. We describe the application of these methods to a number of systems with molecular weights in the hundreds of kilodaltons.

Conformation-specific binding of alpha-synuclein to novel protein partners detected by phage display and NMR spectroscopy

Alpha-synuclein (AS) is an intrinsically unstructured protein in aqueous solution but is capable of forming beta-sheet-rich fibrils that accumulate as intracytoplasmic inclusions in Parkinson disease and certain other neurological disorders. However, AS binding to phospholipid membranes leads to a distinct change in protein conformation, stabilizing an extended amphipathic alpha-helical domain reminiscent of the exchangeable apolipoproteins. To better understand the significance of this conformational change, we devised a novel bacteriophage display screen to identify protein binding partners of helical AS and have identified 20 proteins with roles in diverse cellular processes related to membrane trafficking, ion channel modulation, redox metabolism, and gene regulation. To verify that the screen identifies proteins with specificity for helical AS, we further characterized one of these candidates, endosulfine alpha (ENSA), a small cAMP-regulated phosphoprotein implicated in the regulation of insulin secretion but also expressed abundantly in the brain. We used solution NMR to probe the interaction between ENSA and AS on the surface of SDS micelles. Chemical shift perturbation mapping experiments indicate that ENSA interacts specifically with residues in the N-terminal helical domain of AS in the presence of SDS but not in aqueous buffer lacking SDS. The ENSA-related protein ARPP-19 (cAMP-regulated phosphoprotein 19) also displays specific interactions with helical AS. These results confirm that the helical N terminus of AS can mediate specific interactions with other proteins and suggest that membrane binding may regulate the physiological activity of AS in vivo.

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.