SARA: a software environment for the analysis of relaxation data acquired with accordion spectroscopy

We present SARA (Software for Accordion Relaxation Analysis), an interactive and user-friendly MATLAB software environment designed for analyzing relaxation data obtained with accordion spectroscopy. Accordion spectroscopy can be used to measure nuclear magnetic resonance (NMR) relaxation rates in a fraction of the time required by traditional methods, yet data analysis can be intimidating and no unified software packages are available to assist investigators. Hence, the technique has not achieved widespread use within the NMR community. SARA offers users a selection of analysis protocols spanning those presented in the literature thus far, with modifications permitting a more general application to crowded spectra such as those of proteins. We discuss the advantages and limitations of each fitting method and suggest a protocol combining the strengths of each procedure to achieve optimal results. In the end, SARA provides an environment for facile extraction of relaxation rates and should promote routine application of accordion relaxation spectroscopy.     

Using nuclear magnetic resonance spectroscopy to study molten globule states of proteins

Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for the study of the structure, dynamics, and folding of proteins in solution. It is particularly powerful when applied to dynamic or flexible systems, such as partially folded molten globule states of proteins, which are not usually amenable to X-ray crystallography. In this article, NMR methods suitable for the detailed characterisation of molten globule states are described. The specific method used to study the molten globule is determined by the quality of the NMR spectrum obtained. Molten globules are characterised by significant levels of secondary structure. Site-specific hydrogen-deuterium exchange experiments can be used to identify residues located in regions of secondary structure in the molten globule. If spectra characterised by sharp peaks are observed for the molten globule then information about secondary structure can be obtained by analysis of (1)H(alpha), (13)C(alpha), (13)C(beta), and (13)CO chemical shifts; this can be supplemented by (15)N relaxation studies. For molten globules characterised by extremely broad peaks (15)N-edited NMR experiments carried out in increasing concentrations of denaturants can be used to study the relative stabilities of different regions of structure. Examples of the application of these methods to the study of the low pH molten globule states of alpha-lactalbumin and apomyoglobin are presented.     

Spectroscopic characterization of a newly isolated cytochrome P450 from Rhodococcus rhodochrous.

Cytochrome P450 (P450) from Rhodococcus rhodochrous have been characterized through circular dichroism and nuclear magnetic resonance (NMR) spectroscopy, both in the substrate-free and substrate-bound forms. The data are compared with those of P450cam and indicate a close similarity of the structure of the active site in the two proteins. The substrate-free species contains low-spin iron(III), while the 2-ethoxyphenol bound species contains high-spin iron(III). The substrate is in slow exchange on the NMR time scale. The binding of CN- has been investigated and the final adduct characterized through NMR spectra. Nuclear relaxation times of the isotropically shifted signals turn out to be shorter than in other heme proteins, both in the high- and in the low-spin species. This is the result of longer electron relaxation times in P450s than in peroxidases and metmyoglobin. This property, as well as the electron paramagnetic resonance (EPR) spectrum of the substrate-free form, are discussed in terms of the presence of the cysteine as the fifth ligand of the iron ion instead of a histidine as it occurs in peroxidases and myoglobin.     

Biological membrane structure by solid-state NMR

Nuclear magnetic resonance (NMR) spectroscopy, and particularly solid-state NMR spectroscopy, is a method of choice to study the structure and dynamics of both the lipid and the protein components of model and biological membranes. Different approaches have been developed to study these systems in which the restricted molecular motions result in broad NMR spectra. This contribution will first present an overview of the different techniques used to study lipid bilayers, namely 31p, 2H and 13C solid-state NMR spectroscopy. On the other hand, the study of the structure of membrane peptides and proteins is a rapidly growing field and several methods developed in the last two decades will be presented. These methods allow the investigation of protein systems for which structural information is often difficult to obtain by techniques such as X-ray diffraction and multidimensional solution NMR.     

Multidimensional NMR methods for protein structure determination

Structural studies of proteins are critical for understanding biological processes at the molecular level. Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for obtaining structural and dynamic information on proteins and protein-ligand complexes. In the present review, methodologies for NMR structure determination of proteins and macromolecular complexes are described. In addition, a number of recent advances that reduce the molecular weight limitations previously imposed on NMR studies of biomolecules are discussed, highlighting applications of these technologies to protein systems studied in our laboratories.     

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, “conservation of experimental data in analysis of Fourier” (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b ₆ f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.     

Practical considerations for investigation of protein conformational dynamics by <Superscript>15</Superscript>N <Emphasis Type="Italic">R</Emphasis><Subscript>1ρ</Subscript> relaxation dispersion

It is becoming increasingly apparent that proteins are not static entities and that their function often critically depends on accurate sampling of multiple conformational states in aqueous solution. Accordingly, the development of methods to study conformational states in proteins beyond their ground-state structure (“excited states”) has crucial biophysical importance. Here we investigate experimental schemes for optimally probing chemical exchange processes in proteins on the micro- to millisecond timescale by ¹⁵N R _1ρ relaxation dispersion. The schemes use selective Hartmann–Hahn cross-polarization (CP) transfer for excitation, and derive peak integrals from 1D NMR spectra (Korzhnev et al. in J Am Chem Soc 127:713–721, 2005; Hansen et al. in J Am Chem Soc 131:3818–3819, 2009). Simulation and experiment collectively show that in such CP-based schemes care has to be taken to achieve accurate suppression of undesired off-resonance coherences, when using weak spin-lock fields. This then (i) ensures that relaxation dispersion profiles in the absence of chemical exchange are flat, and (ii) facilitates extraction of relaxation dispersion profiles in crowded regions of the spectrum. Further improvement in the quality of the experimental data is achieved by recording the free-induction decays in an interleaved manner and including a heating-compensation element. The reported considerations will particularly benefit the use of CP-based R _1ρ relaxation dispersion to analyze conformational exchange processes in larger proteins, where resonance line overlap becomes the main limiting factor.     

Fast and accurate fitting of relaxation dispersion data using the flexible software package GLOVE

Relaxation dispersion spectroscopy is one of the most widely used techniques for the analysis of protein dynamics. To obtain a detailed understanding of the protein function from the view point of dynamics, it is essential to fit relaxation dispersion data accurately. The grid search method is commonly used for relaxation dispersion curve fits, but it does not always find the global minimum that provides the best-fit parameter set. Also, the fitting quality does not always improve with increase of the grid size although the computational time becomes longer. This is because relaxation dispersion curve fitting suffers from a local minimum problem, which is a general problem in non-linear least squares curve fitting. Therefore, in order to fit relaxation dispersion data rapidly and accurately, we developed a new fitting program called GLOVE that minimizes global and local parameters alternately, and incorporates a Monte-Carlo minimization method that enables fitting parameters to pass through local minima with low computational cost. GLOVE also implements a random search method, which sets up initial parameter values randomly within user-defined ranges. We demonstrate here that the combined use of the three methods can find the global minimum more rapidly and more accurately than grid search alone.     

Backbone dynamics of the N-terminal domain of the HIV-1 capsid protein and comparison with the G94D mutant conferring cyclosporin resistance/dependence.

Nuclear magnetic resonance (NMR) (15)N relaxation methods have been used to characterize the backbone dynamics of the N-terminal core domain of the HIV-1 capsid protein (CA(151)). The domain, which has an unusually flat, triangular shape, tumbles in solution at 28 degrees C with an effective rotational correlation time of 11.5 ns. Relaxation data for backbone amides in the domain's seven alpha-helices are indicative of fully anisotropic rotational diffusion. The principal axes of the rotational diffusion tensor calculated from the NMR data are aligned to within 12-23 degrees of the principal axes of the inertial tensor, with the axis of fastest rotational diffusion coincident with that of minimal inertia, and vice versa. Large variations in the (15)N-(1)H nuclear Overhauser effects for individual amino acids correlate with the degree of convergence in the previously calculated NMR structure. In particular, the partially disordered residues Val86-Arg97 that contain the human cyclophilin A (CypA) packaging signal have (15)N heteronuclear NOEs and transversal relaxation rates consistent with a high degree of dynamic conformational averaging. The N-terminal domain of a CA mutant (G94D) that confers both resistance to and dependence on cyclosporin A analogues was also analyzed. Our results indicate that this mutation does not influence the conformation or dynamics of CA(151), and therefore probably affects the function of the protein by modifying essential intermolecular CA-CA interactions.     

Thermodynamic interpretation of protein dynamics from NMR relaxation measurements

Protein dynamics and thermodynamics can be characterized through measurements of relaxation rates of side chain (2)H and (13)C, and backbone (15)N nuclei using NMR spectroscopy. The rates reflect protein motions on timescales from picoseconds to milliseconds. Backbone and methyl side chain NMR relaxation measurements for several proteins are beginning to reveal the role of protein dynamics in protein stability and ligand binding.     

Abp1p and Fyn SH3 domains fold through similar low-populated intermediate states

Src homology 3 (SH3) domains are small modules that are thought to fold via a two-state mechanism, without the accumulation of significant populations of intermediate states. Relaxation dispersion NMR studies of the folding of G48V and G48M mutants of the Fyn SH3 domain have established that, at least for these modules, folding proceeds through the formation of a transient on-pathway intermediate with an equilibrium population of 1-2% that can be readily detected [Korzhnev, D. M., et al. (2004) Nature 430, 586-590]. To investigate the generality of this result, we present an (15)N relaxation dispersion NMR study of a pair of additional SH3 domains, including a G48V mutant of a stabilized Abp1p SH3 domain that shares 36% sequence identity with the Fyn SH3 module, and a A39V/N53P/V55L mutant Fyn SH3 domain. A transient folding intermediate is detected for both of the proteins studied here, and the dispersion data are well fit to a folding model of the form F <--> I <--> U, where F, I, and U correspond to folded, intermediate, and unfolded states, respectively. The temperature dependencies of the folding/unfolding rate constants were obtained so that the thermodynamic properties of each of F, I, and U could be established. The detection of I states in folding pathways of all SH3 domains examined to date via relaxation dispersion NMR spectroscopy indicates that such intermediates may well be a conserved feature in the folding of such domains in general but that their transient nature along with their low population makes detection difficult using more well-established approaches to the study of folding.     

Fragment-based drug discovery using NMR spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy has evolved into a powerful tool for fragment-based drug discovery over the last two decades. While NMR has been traditionally used to elucidate the three-dimensional structures and dynamics of biomacromolecules and their interactions, it can also be a very valuable tool for the reliable identification of small molecules that bind to proteins and for hit-to-lead optimization. Here, we describe the use of NMR spectroscopy as a method for fragment-based drug discovery and how to most effectively utilize this approach for discovering novel therapeutics based on our experience.     

Can enzyme engineering benefit from the modulation of protein motions? Lessons learned from NMR relaxation dispersion experiments

Despite impressive progress in protein engineering and design, our ability to create new and efficient enzyme activities remains a laborious and time-consuming endeavor. In the past few years, intricate combinations of rational mutagenesis, directed evolution and computational methods have paved the way to exciting engineering examples and are now offering a new perspective on the structural requirements of enzyme activity. However, these structure-function analyses are usually guided by the time-averaged static models offered by enzyme crystal structures, which often fail to describe the functionally relevant 'invisible states' adopted by proteins in space and time. To alleviate such limitations, NMR relaxation dispersion experiments coupled to mutagenesis studies have recently been applied to the study of enzyme catalysis, effectively complementing 'structure-function' analyses with 'flexibility-function' investigations. In addition to offering quantitative, site-specific information to help characterize residue motion, these NMR methods are now being applied to enzyme engineering purposes, providing a powerful tool to help characterize the effects of controlling long-range networks of flexible residues affecting enzyme function. Recent advancements in this emerging field are presented here, with particular attention to mutagenesis reports highlighting the relevance of NMR relaxation dispersion tools in enzyme engineering.     

NMR Study of Chloroplast Water in Acer platanoides: Water Exchange at Membrane-Bound Sites and Across the Chloroplast Envelope Membrane

Abstract: Nuclear magnetic resonance (NMR) spin relaxation techniques were used to measure the rate of water exchange across the chloroplast envelope membrane and to detect exchange of water-hydrogens between the bound and fluid states in Acer platanoides shade leaves. Results indicate that individual water molecules resided within a chloroplast for 320 ± 65 ms, on average, before exchanging through the envelope, and that exchangable bound hydrogens are at least six times more abundant in the chloroplast than in the nonchloroplast compartments.     

Structural Insight into the Transmembrane Domain and the Juxtamembrane Region of the Erythropoietin Receptor in Micelles

Erythropoietin receptor (EpoR) dimerization is an important step in erythrocyte formation. Its transmembrane domain (TMD) and juxtamembrane (JM) region are essential for signal transduction across the membrane. A construct compassing residues S212–P259 and containing the TMD and JM region of the human EpoR was purified and reconstituted in detergent micelles. The solution structure of the construct was determined in dodecylphosphocholine (DPC) micelles by solution NMR spectroscopy. Structural and dynamic studies demonstrated that the TMD and JM region are an α-helix in DPC micelles, whereas residues S212–D224 at the N-terminus of the construct are not structured. The JM region is a helix that contains a hydrophobic patch formed by conserved hydrophobic residues (L253, I257, and W258). Nuclear Overhauser effect analysis, fluorescence spectroscopy, and paramagnetic relaxation enhancement experiments suggested that the JM region is exposed to the solvent. The structures of the TMD and JM region of the mouse EpoR were similar to those of the human EpoR.