Protein dynamics from NMR

This review surveys recent investigations of conformational fluctuations of proteins in solution using NMR techniques. Advances in experimental methods have provided more accurate means of characterizing fast and slow internal motions as well as overall diffusion. The information obtained from NMR dynamics experiments provides insights into specific structural changes or configurational energetics associated with function. A variety of applications illustrate that studies of protein dynamics provide insights into protein-protein interactions, target recognition, ligand binding, and enzyme function.     

Dynamic coupling and allosteric behavior in a nonallosteric protein

Long-range intraprotein interactions give rise to many important protein behaviors. Understanding how energy is transduced through protein structures to either transmit a signal or elicit conformational changes is therefore a current challenge in structural biology. In an effort to understand such linkages, multiple V --> A mutations were made in the small globular protein eglin c. The physical responses, as mapped by NMR spin relaxation, residual dipolar couplings (RDCs), and scalar couplings, illustrate that the interior of this nonallosteric protein forms a dynamic network and that local perturbations are transmitted as dynamic and structural changes to distal sites as far as 16 A away. Two basic types of propagation responses were observed: contiguous pathways of enhanced (attenuated) dynamics with no change in structure; and dispersed (noncontiguous) changes in methyl rotation rates that appear to result from subtle deformation of backbone structure. In addition, energy transmission is found to be unidirectional. In one mutant, an allosteric conformational change of a side chain is seen in the context of a pathway of propagated changes in picosecond to nanosecond dynamics. The observation of so many long-range interactions in a small, rigid system lends experimental weight to the idea that all well-folded proteins inherently possess allosteric features [Gunasekaran et al. (2004) Proteins 57, 433-443], and that dynamics are a rich source of information for mapping and gaining mechanistic insight into communication pathways in individual proteins.     

A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Guanine nucleotides differentially modulate backbone dynamics of the STAS domain of the SulP/SLC26 transport protein Rv1739c of Mycobacterium tuberculosis

Enzymatic catalysis and protein signaling are dynamic processes that involve local and/or global conformational changes occurring across a broad range of time scales. (1) H-(15) N relaxation NMR provides a comprehensive understanding of protein backbone dynamics both in the apo (unliganded) and ligand-bound conformations, enabling both fast and slow internal motions of individual amino acid residues to be observed. We recently reported the structure and nucleotide binding properties of the sulfate transporter and anti-sigma factor antagonist (STAS) domain of Rv1739c, a SulP anion transporter protein of Mycobacterium tuberculosis. In the present study, we report (1) H-(15) N NMR backbone dynamics measurements [longitudinal (T(1) ), transverse (T(2) ) and steady-state ({(1) H}-(15) N) heteronuclear NOE] of the Rv1739c STAS domain, in the absence and presence of saturating concentrations of GTP and GDP. Analysis of measured relaxation data and estimated dynamic parameters indicated distinct features differentiating the binding of GTP and GDP to Rv1739c STAS. The 9.55 ns overall rotational correlation time of Rv1739c STAS increased to 10.48 ns in the presence of GTP, and to 13.25 ns in the presence of GDP, indicating significant nucleotide-induced conformational changes. These conformational changes were accompanied by slow time scale (μs to ms) motions in discrete regions of the protein, as reflected by guanine nucleotide-induced changes in relaxation parameters. The observed nucleotide-specific alterations in the relaxation properties of individual STAS residues reflect an increased molecular anisotropy and/or the emergence of conformational equilibria governing functional properties of the STAS domain.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Protein conformations explored by difference high-angle solution X-ray scattering: oxidation state and temperature dependent changes in cytochrome C

We demonstrate the use of high-angle X-ray scattering to explore protein conformational states in solution by resolving oxidation state- and temperature-dependent changes in the conformation of horse heart cytochrome c. Several detailed models exist for oxidation-dependent changes in mitochondrial class I c cytochromes determined by X-ray crystallography and solution NMR techniques. These models differ in the magnitude and locations of structural change. Our scattering measurements show that high-angle X-ray scattering can discriminate between these models, and that the experimental scattering data for horse cytochrome c can be best reconciled with selected NMR models for the same protein. These results demonstrate the ability to use high-angle X-ray scattering to resolve conformational states of proteins in solution, and to relate these measurements to detailed structural models. Furthermore, temperature-dependent changes are found in the high angle scattering patterns for horse cytochrome c, illustrating the sensitivity of these measurements to dynamic aspects of protein structure. These results demonstrate the ability to use difference high angle scattering as a quantitative monitor of reaction-linked changes in protein conformation and structural dynamics. Synchrotron-based high-angle scattering holds promise as a widely applicable, high throughput technique for exploring conformational states linked to physiological protein function, for resolving configurational differences between protein structures in solution and crystalline states, and for bridging the gap between solution NMR and crystallographic structure techniques.     

Effect of pressure on the tertiary structure and dynamics of folded basic pancreatic trypsin inhibitor

The on-line high-pressure cell NMR technique was used to study pressure-induced changes in the tertiary structure and dynamics of a globular protein, basic pancreatic trypsin inhibitor (BPTI). Practically all the proton signals of BPTI were observed with (1)H two-dimensional NMR spectroscopy at 750 MHz at variable pressure between 1 and 2000 bar. Chemical shifts, nuclear Overhauser effect (NOE), and line shapes were used to analyze conformational and dynamic changes of the protein as functions of pressure. Linear, reversible, but nonuniform pressure-induced chemical shift changes of practically all the C(alpha) protons and side chain protons showed that the entire secondary and tertiary structures are altered by pressure within the folded ensemble of BPTI. The high field shift tendency of most side chain proton signals and the increase in NOE intensities of some specific side chain protons indicated a site-specific compaction of the tertiary structure. Pressure dependence of ring flip rates was deduced from resonance line shapes of the slices of the two-dimensional NMR spectrum for ring proton signals of Tyr-35 and Phe-45. The rates of the flip-flop motions were considerably reduced at high pressure, from which activation volumes were determined to be 85 +/- 20 A(3) (or 51.2 ml/mol) and 46 +/- 9 A(3) (or 27.7 ml/mol) for Tyr-35 and Phe-45, respectively, at 57 degrees C. The present experiments confirm that pressure affects the entire secondary and tertiary structures of a globular protein with specific compaction of a core, leading to quite significant changes in slow internal dynamics of a globular protein.     

The dynamic duo: Combining NMR and small angle scattering in structural biology

Structural biology provides essential information for elucidating molecular mechanisms that underlie biological function. Advances in hardware, sample preparation, experimental methods, and computational approaches now enable structural analysis of protein complexes with increasing complexity that more closely represent biologically entities in the cellular environment. Integrated multidisciplinary approaches are required to overcome limitations of individual methods and take advantage of complementary aspects provided by different structural biology techniques. Although X‐ray crystallography remains the method of choice for structural analysis of large complexes, crystallization of flexible systems is often difficult and does typically not provide insights into conformational dynamics present in solution. Nuclear magnetic resonance spectroscopy (NMR) is well‐suited to study dynamics at picosecond to second time scales, and to map binding interfaces even of large systems at residue resolution but suffers from poor sensitivity with increasing molecular weight. Small angle scattering (SAS) methods provide low resolution information in solution and can characterize dynamics and conformational equilibria complementary to crystallography and NMR. The combination of NMR, crystallography, and SAS is, thus, very useful for analysis of the structure and conformational dynamics of (large) protein complexes in solution. In high molecular weight systems, where NMR data are often sparse, SAS provides additional structural information and can differentiate between NMR‐derived models. Scattering data can also validate the solution conformation of a crystal structure and indicate the presence of conformational equilibria. Here, we review current state‐of‐the‐art approaches for combining NMR, crystallography, and SAS data to characterize protein complexes in solution.     

Thermal Adaptation of Conformational Dynamics in Ribonuclease H

The relationship between inherent internal conformational processes and enzymatic activity or thermodynamic stability of proteins has proven difficult to characterize. The study of homologous proteins with differing thermostabilities offers an especially useful approach for understanding the functional aspects of conformational dynamics. In particular, ribonuclease HI (RNase H), an 18 kD globular protein that hydrolyzes the RNA strand of RNA:DNA hybrid substrates, has been extensively studied by NMR spectroscopy to characterize the differences in dynamics between homologs from the mesophilic organism E. coli and the thermophilic organism T. thermophilus. Herein, molecular dynamics simulations are reported for five homologous RNase H proteins of varying thermostabilities and enzymatic activities from organisms of markedly different preferred growth temperatures. For the E. coli and T. thermophilus proteins, strong agreement is obtained between simulated and experimental values for NMR order parameters and for dynamically averaged chemical shifts, suggesting that these simulations can be a productive platform for predicting the effects of individual amino acid residues on dynamic behavior. Analyses of the simulations reveal that a single residue differentiates between two different and otherwise conserved dynamic processes in a region of the protein known to form part of the substrate-binding interface. Additional key residues within these two categories are identified through the temperature-dependence of these conformational processes.     

Observation of slow dynamic exchange processes in Ras protein crystals by 31P solid state NMR spectroscopy

The folding, structure and biological function of many proteins are inherently dynamic properties of the protein molecule. Often, the respective molecular processes are preserved upon protein crystallization, leading, in X-ray diffraction experiments, to a blurring of the electron density map and reducing the resolution of the derived structure. Nuclear magnetic resonance (NMR) is known to be an alternative method to study molecular structure and dynamics. We designed and built a probe for phosphorus solid state NMR that allows for the first time to study static properties as well as dynamic processes in single-crystals of a protein by NMR spectroscopy. The sensitivity achieved is sufficient to detect the NMR signal from individual phosphorus sites in a 0.3mm(3) size single-crystal of GTPase Ras bound to the nucleotide GppNHp, that is, the signal from approximately 10(15) phosphorus nuclei. The NMR spectra obtained are discussed in terms of the conformational variability of the active center of the Ras-nucleotide complex. We conclude that, in the crystal, the protein complex exists in three different conformations. Magic angle spinning (MAS) NMR spectra of a powder sample of Ras-GppNHp show a splitting of one of the phosphate resonances and thus confirm this conclusion. The MAS spectra provide, furthermore, evidence of a slow, temperature-dependent dynamic exchange process in the Ras protein crystal.     

Conformational Dynamics of Light-Harvesting Complex II in a Native Membrane Environment

Photosynthetic light-harvesting complexes (LHCs) of higher plants, moss, and green algae can undergo dynamic conformational transitions, which have been correlated to their ability to adapt to fluctuations in the light environment. Herein, we demonstrate the application of solid-state NMR spectroscopy on native, heterogeneous thylakoid membranes of Chlamydomonas reinhardtii (Cr) and on Cr light-harvesting complex II (LHCII) in thylakoid lipid bilayers to detect LHCII conformational dynamics in its native membrane environment. We show that membrane-reconstituted LHCII contains selective sites that undergo fast, large-amplitude motions, including the phytol tails of two chlorophylls. Protein plasticity is also observed in the N-terminal stromal loop and in protein fragments facing the lumen, involving sites that stabilize the xanthophyll-cycle carotenoid violaxanthin and the two luteins. The results report on the intrinsic flexibility of LHCII pigment-protein complexes in a membrane environment, revealing putative sites for conformational switching. In thylakoid membranes, fast dynamics of protein and pigment sites is significantly reduced, which suggests that in their native organelle membranes, LHCII complexes are locked in specific conformational states.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Conformational ensemble of a multidomain protein explored by Gd3+ electron paramagnetic resonance

Despite their importance in function, the conformational state of proteins and its changes are often poorly understood, mainly because of the lack of an efficient tool. MurD, a 47-kDa protein enzyme responsible for peptidoglycan biosynthesis, is one of those proteins whose conformational states and changes during their catalytic cycle are not well understood. Although it has been considered that MurD takes a single conformational state in solution as shown by a crystal structure, the solution nuclear magnetic resonance (NMR) study suggested the existence of multiple conformational state of apo MurD in solution. However, the conformational distribution has not been evaluated. In this work, we investigate the conformational states of MurD by the use of electron paramagnetic resonance (EPR), especially intergadolinium distance measurement using double electron-electron resonance (DEER) measurement. The gadolinium ions are fixed on specific positions on MurD via a rigid double-arm paramagnetic lanthanide tag that has been originally developed for paramagnetic NMR. The combined use of NMR and EPR enables accurate interpretation of the DEER distance information to the structural information of MurD. The DEER distance measurement for apo MurD shows a broad distance distribution, whereas the presence of the inhibitor narrows the distance distribution. The results suggest that MurD exists in a wide variety of conformational states in the absence of ligands, whereas binding of the inhibitor eliminates variation in conformational states. The multiple conformational states of MurD were previously implied by NMR experiments, but our DEER data provided structural characterization of the conformational variety of MurD.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Localized conformational changes trigger the pH-induced fibrillogenesis of an amyloidogenic λ light chain protein

Solvent conditions modulate the expression of the amyloidogenic potential of proteins. In this work the effect of pH on the fibrillogenic behavior and the conformational properties of 6aJL2, a model protein of the highly amyloidogenic variable light chain λ6a gene segment, was examined. Ordered aggregates showing the ultrastructural and spectroscopic properties observed in amyloid fibrils were formed in the 2.0–8.0?pH range. At pH <3.0 a drastic decrease in lag time and an increase in fibril formation rate were found. In the 4.0–8.0?pH range there was no spectroscopic evidence for significant conformational changes in the native state. Likewise, heat capacity measurements showed no evidence for residual structure in the unfolded state. However, at pH <3.0 stability is severely decreased and the protein suffers conformational changes as detected by circular dichroism, tryptophan and ANS fluorescence, as well as by NMR spectroscopy. Molecular dynamics simulations indicate that acid-induced conformational changes involve the exposure of the loop connecting strands E and F. These results are compatible with pH-induced changes in the NMR spectra. Overall, the results indicate that the mechanism involved in the acid-induced increase in the fibrillogenic potential of 6aJL2 is profoundly different to that observed in κ light chains, and is promoted by localized conformational changes in a region of the protein that was previously not known to be involved in acid-induced light chain fibril formation. The identification of this region opens the potential for the design of specific inhibitors.     

Characterization of dynamic behavior of side groups of globular proteins by spin-label, NMR and luminescence polarization methods

In this work we suggest a quantitative estimation of a complicated motion of side groups of globular proteins. In the general case, three basic parameters determine the motion: (a) rotational correlation time of a side unit under study, or covalently bound spin label, or dye, (b) parameter S that reflects sterical restrictions for re-orientation of the given unit (these two parameters depending on the side-chain structure and its conformational change within the immediate dynamic protein surrounding whereas correlation times of side units on microviscosity in addition), (c) rotational correlation time of protein globule. These parameters can be measured by spin-label, NMR and fluorescence polarization techniques. An attempt to describe a complicated dynamic behaviour of side units of protein macromolecules with a single dynamics parameter--rotational correlation time--not only leads to a loss of part of information about the local structural dynamics of macromolecules but also can diminish the tau value.     

Conformational fluctuations of proteins revealed by variable pressure NMR

With the high-resolution variable-pressure NMR spectroscopy, one can study conformational fluctuations of proteins in a much wider conformational space than hitherto explored by NMR and other spectroscopic techniques. This is because a protein in solution generally exists as a dynamic mixture of conformers mutually differing in partial molar volume, and pressure can select the population of a conformer according to its relative volume. In this review, we describe how variable-pressure NMR can be used to probe conformational fluctuations of proteins in a wide conformational space from the folded to the fully unfolded structures, with actual examples. Furthermore, the newly emerging technique "NMR snapshots" expresses amply fluctuating protein structures as changes in atomic coordinates. Finally, the concept of conformational fluctuation is extended to include intermolecular association leading to amyloidosis.     

Reconciling Mediating and Slaving Roles of Water in Protein Conformational Dynamics

Proteins accomplish their physiological functions with remarkably organized dynamic transitions among a hierarchical network of conformational substates. Despite the essential contribution of water molecules in shaping functionally important protein dynamics, their exact role is still controversial. Water molecules were reported either as mediators that facilitate or as masters that slave protein dynamics. Since dynamic behaviour of a given protein is ultimately determined by the underlying energy landscape, we systematically analysed protein self energies and protein-water interaction energies obtained from extensive molecular dynamics simulation trajectories of barstar. We found that protein-water interaction energy plays the dominant role when compared with protein self energy, and these two energy terms on average have negative correlation that increases with increasingly longer time scales ranging from 10 femtoseconds to 100 nanoseconds. Water molecules effectively roughen potential energy surface of proteins in the majority part of observed conformational space and smooth in the remaining part. These findings support a scenario wherein water on average slave protein conformational dynamics but facilitate a fraction of transitions among different conformational substates, and reconcile the controversy on the facilitating and slaving roles of water molecules in protein conformational dynamics.