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.     

Intein-based biosynthetic incorporation of unlabeled protein tags into isotopically labeled proteins for NMR studies

Segmental isotopic labeling of proteins using protein ligation is a recently established in vitro method for incorporating isotopes into one domain or region of a protein to reduce the complexity of NMR spectra, thereby facilitating the NMR analysis of larger proteins. Here we demonstrate that segmental isotopic labeling of proteins can be conveniently achieved in Escherichia coli using intein-based protein ligation. Our method is based on a dual expression system that allows sequential expression of two precursor fragments in media enriched with different isotopes. Using this in vivo approach, unlabeled protein tags can be incorporated into isotopically labeled target proteins to improve protein stability and solubility for study by solution NMR spectroscopy.     

Resolution-by-proxy: a simple measure for assessing and comparing the overall quality of NMR protein structures

In protein X-ray crystallography, resolution is often used as a good indicator of structural quality. Diffraction resolution of protein crystals correlates well with the number of X-ray observables that are used in structure generation and, therefore, with protein coordinate errors. In protein NMR, there is no parameter identical to X-ray resolution. Instead, resolution is often used as a synonym of NMR model quality. Resolution of NMR structures is often deduced from ensemble precision, torsion angle normality and number of distance restraints per residue. The lack of common techniques to assess the resolution of X-ray and NMR structures complicates the comparison of structures solved by these two methods. This problem is sometimes approached by calculating "equivalent resolution" from structure quality metrics. However, existing protocols do not offer a comprehensive assessment of protein structure as they calculate equivalent resolution from a relatively small number (<5) of protein parameters. Here, we report a development of a protocol that calculates equivalent resolution from 25 measurable protein features. This new method offers better performance (correlation coefficient of 0.92, mean absolute error of 0.28 ?) than existing predictors of equivalent resolution. Because the method uses coordinate data as a proxy for X-ray diffraction data, we call this measure "Resolution-by-Proxy" or ResProx. We demonstrate that ResProx can be used to identify under-restrained, poorly refined or inaccurate NMR structures, and can discover structural defects that the other equivalent resolution methods cannot detect. The ResProx web server is available at http://www.resprox.ca.     

High-Resolution 1H NMR Spectroscopy of Fish Muscle,Eggs and Small Whole Fish via Hadamard-Encoded Intermolecular Multiple-Quantum Coherence

Background and Purpose

Nuclear magnetic resonance (NMR) spectroscopy has become an important technique for tissue studies. Since tissues are in semisolid-state, their high-resolution (HR) spectra cannot be obtained by conventional NMR spectroscopy. Because of this restriction, extraction and high-resolution magic angle spinning (HR MAS) are widely applied for HR NMR spectra of tissues. However, both of the methods are subject to limitations. In this study, the feasibility of HR ¹H NMR spectroscopy based on intermolecular multiple-quantum coherence (iMQC) technique is explored using fish muscle, fish eggs, and a whole fish as examples.

Materials and Methods

Intact salmon muscle tissues, intact eggs from shishamo smelt and a whole fish (Siamese algae eater) are studied by using conventional 1D one-pulse sequence, Hadamard-encoded iMQC sequence, and HR MAS.

Results

When we use the conventional 1D one-pulse sequence, hardly any useful spectral information can be obtained due to the severe field inhomogeneity. By contrast, HR NMR spectra can be obtained in a short period of time by using the Hadamard-encoded iMQC method without shimming. Most signals from fatty acids and small metabolites can be observed. Compared to HR MAS, the iMQC method is non-invasive, but the resolution and the sensitivity of resulting spectra are not as high as those of HR MAS spectra.

Conclusion

Due to the immunity to field inhomogeneity, the iMQC technique can be a proper supplement to HR MAS, and it provides an alternative for the investigation in cases with field distortions and with samples unsuitable for spinning. The acquisition time of the proposed method is greatly reduced by introduction of the Hadamard-encoded technique, in comparison with that of conventional iMQC method.     

Naturally-occurring pituitary growth hormone is phosphorylated

The results of this communication show that ovine growth hormone (oGH) contains organically-bound phosphorous. The phosphorous content of growth hormone, lot S-11, is 1:3 (mol/mol) and that of lot S-12 is 1:6 (mol/mol). Results of 31P NMR studies suggest that the phosphorous exists in two chemical forms: as a monophosphoryl ester and as a phosphodiester. Evidence is provided which demonstrates that growth hormone can be phosphorylated in vitro with the catalytic subunit of protein kinase.     

COCO: a simple tool to enrich the representation of conformational variability in NMR structures

NMR structures are typically deposited in databases such as the PDB in the form of an ensemble of structures. Generally, each of the models in such an ensemble satisfies the experimental data and is equally valid. No unique solution can be calculated because the experimental NMR data is insufficient, in part because it reflects the conformational variability and dynamical behavior of the molecule in solution. Even for relatively rigid molecules, the limited number of structures that are typically deposited cannot completely encompass the structural diversity allowed by the observed NMR data, but they can be chosen to try and maximize its representation. We describe here the adaptation and application of techniques more commonly used to examine large ensembles from molecular dynamics simulations, to the analysis of NMR ensembles. The approach, which is based on principal component analysis, we call COCO ("Complementary Coordinates"). The COCO approach analyses the distribution of an NMR ensemble in conformational space, and generates a new ensemble that fills "gaps" in the distribution. The method is very rapid, and analysis of a 25-member ensemble and generation of a new 25 member ensemble typically takes 1-2 min on a conventional workstation. Applied to the 545 structures in the RECOORD database, we find that COCO generates new ensembles that are as structurally diverse-both from each other and from the original ensemble-as are the structures within the original ensemble. The COCO approach does not explicitly take into account the NMR restraint data, yet in tests on selected structures from the RECOORD database, the COCO ensembles are frequently good matches to this data, and certainly are structures that can be rapidly refined against the restraints to yield high-quality, novel solutions. COCO should therefore be a useful aid in NMR structure refinement and in other situations where a richer representation of conformational variability is desired-for example in docking studies. COCO is freely accessible via the website www.ccpb.ac.uk/COCO.     

Isotope-edited multidimensional NMR of calcineurin B in the presence of the non-deuterated detergent CHAPS

Summary At the concentration needed for NMR, the calcium-saturated form of calcineurin B dissolved in water shows resonance line widths that indicate aggregation of this protein. Although the line width or aggregation state can be influenced to some degree by temperature, pH, and salt concentrations, in the absence of detergent no conditions could be found where the protein behaved as a monomeric unit. In the presence of a 10- to 20-fold molar excess of the zwitterionic detergent 3-[(3-cholamidopropyl)-dimethyl-ammonio]-1-propanesulfonate (CHAPS), resonance line widths were considerably narrower and were compatible with a protein of 25 kDa. The presence of the NMR signals of the non-deuterated CHAPS does not interfere with modern isotope-directed NMR studies as the signals from protons not attached to ¹⁵N or ¹³C are removed by isotope filtering and purge pulses.     

Characterization of the yeast peroxiredoxin Ahp1 in its reduced active and overoxidized inactive forms using NMR

Peroxiredoxins (Prx's) are a superfamily of thiol-specific antioxidant proteins present in all organisms and involved in the hydroperoxide detoxification of the cell. The catalytic cysteine of Prx's reduces hydroperoxides and is transformed into a transient sulfenic acid (Cys-SOH). At high hydroperoxide concentration, the sulfenic acid can be overoxidized into a sulfinate, or even a sulfonate. We present here the first peroxiredoxin characterization by solution NMR of the Saccharomyces cerevisiae alkylhydroperoxide reductase (Ahp1) in its reduced and in vitro overoxidized forms. NMR (15)N relaxation data and ultracentrifugation experiments indicate that the protein behaves principally as a homodimer (2 x 19 kDa) in solution, regardless of the redox state. In vitro treatment of Ahp1 by a large excess of tBuOOH leads to an inactive form, with the catalytic cysteine overoxidized into sulfonate, as demonstrated by (13)C NMR. Depending on the amino acid sequence of their active site, Prx's are classified into five different families. In this classification, Ahp1 is a member of the scarcely studied D-type Prx's. Ahp1 is unique among the D-type Prx's in its ability to form an intermolecular disulfide. The peptidic sequence of Ahp1 was analyzed and compared to other D-type Prx sequences.     

The NMR structure of dematin headpiece reveals a dynamic loop that is conformationally altered upon phosphorylation at a distal site

Dematin (band 4.9) is found in the junctional complex of the spectrin cytoskeleton that supports the erythrocyte cell membrane. Dematin is a member of the larger class of cytoskeleton-associated proteins that contain a modular "headpiece" domain at their extreme C termini. The dematin headpiece domain provides the second F-actin-binding site required for in vitro F-actin bundling. The dematin headpiece is found in two forms in the cell, one of 68 residues (DHP) and one containing a 22-amino acid insert near its N terminus (DHP+22). In addition, dematin contains the only headpiece domain that is phosphorylated, in vivo. The 22-amino acid insert in DHP+22 appeared unstructured in NMR spectra; therefore, we have determined the three-dimensional structure of DHP by multidimensional NMR methods. Although the overall three-dimensional structure of DHP is similar to that of the villin headpiece, there are two novel characteristics revealed by this structure. First, unlike villin headpiece that contains a single buried salt bridge, DHP contains a buried charged cluster comprising residues Glu(39), Arg(66), Lys(70), and the C-terminal carboxylate of Phe(76). Second, (15)N relaxation experiments indicate that the longer "variable loop" region near the N terminus of DHP (residues 20-29) is dynamic, undergoing significantly greater motions that the rest of the structure. Furthermore, NMR chemical shift changes indicate that the conformation of the dynamic variable loop is altered by phosphorylation of serine 74, which is far in the sequence from the variable loop region. Our results suggest that phosphorylation of the dematin headpiece acts as a conformational switch within this headpiece domain.     

Direct NMR analysis of cannabis water extracts and tinctures and semi-quantitative data on delta9-THC and delta9-THC-acid

Cannabis sativa L. is the source for a whole series of chemically diverse bioactive compounds that are currently under intensive pharmaceutical investigation. In this work, hot and cold water extracts as well as ethanol/water mixtures (tinctures) of cannabis were compared in order to better understand how these extracts differ in their overall composition. NMR analysis and in vitro cell assays of crude extracts and fractions were performed. Manufacturing procedures to produce natural remedies can strongly affect the final composition of the herbal medicines. Temperature and polarity of the solvents used for the extraction resulted to be two factors that affect the total amount of Delta(9)-THC in the extracts and its relative quantity with respect to Delta(9)-THC-acid and other metabolites. Diffusion-edited (1)H NMR (1D DOSY) and (1)H NMR with suppression of the ethanol and water signals were used. With this method it was possible, without any evaporation or separation step, to distinguish between tinctures from different cannabis cultivars. This approach is proposed as a direct analysis of plant tinctures.     

Rapid prediction of multi-dimensional NMR data sets

We present a computational environment for Fast Analysis of multidimensional NMR DAta Sets (FANDAS) that allows assembling multidimensional data sets from a variety of input parameters and facilitates comparing and modifying such ??in silico?? data sets during the various stages of the NMR data analysis. The input parameters can vary from (partial) NMR assignments directly obtained from experiments to values retrieved from in silico prediction programs. The resulting predicted data sets enable a rapid evaluation of sample labeling in light of spectral resolution and structural content, using standard NMR software such as Sparky. In addition, direct comparison to experimental data sets can be used to validate NMR assignments, distinguish different molecular components, refine structural models or other parameters derived from NMR data. The method is demonstrated in the context of solid-state NMR data obtained for the cyclic nucleotide binding domain of a bacterial cyclic nucleotide-gated channel and on membrane-embedded sensory rhodopsin II. FANDAS is freely available as web portal under WeNMR (http://www.wenmr.eu/services/FANDAS).     

Practical Model Fitting Approaches to the Direct Extraction of NMR Parameters Simultaneously from All Dimensions of Multidimensional NMR Spectra

A maximum likelihood (ML)-based approach has been established for the direct extraction of NMR parameters (e.g., frequency, amplitude, phase, and decay rate) simultaneously from all dimensions of a D-dimensional NMR spectrum. The approach, referred to here as HTFD-ML (hybrid time frequency domain maximum likelihood), constructs a time-domain model composed of a sum of exponentially-decaying sinusoidal signals. The apodized Fourier transform of this time-domain signal is a model spectrum that represents the best fit to the equivalent frequency-domain data spectrum. The desired amplitude and frequency parameters can be extracted directly from the signal model constructed by the HTFD-ML algorithm. The HTFD-ML approach presented here, as embodied in the software package CHIFIT, is designed to meet the challenges posed by model fitting of D-dimensional NMR data sets, where each consists of many data points (108 is not uncommon) encoding information about numerous signals (up to 105 for a protein of moderate size) that exhibit spectral overlap. The suitability of the approach is demonstrated by its application to the concerted analysis of a series of ten 2D 1H-15N HSQC experiments measuring 15N T1 relaxation. In addition to demonstrating the practicality of performing maximum likelihood analysis on large, multidimensional NMR spectra, the results demonstrate that this parametric model-fitting approach provides more accurate amplitude and frequency estimates than those obtained from conventional peak-based analysis of the FT spectrum. The improved performance of the model fitting approach derives from its ability to take into account the simultaneous contributions of all signals in a crowded spectral region (deconvolution) as well as to incorporate prior knowledge in constructing models to fit the data.     

NMR characterization of a pH-dependent equilibrium between two folded solution conformations of the pheromone-binding protein from Bombyx mori

NMR spectroscopic changes as a function of pH in solutions of the pheromone-binding protein of Bombyx mori (BmPBP) show that BmPBP undergoes a conformational transition between pH 4.9 and 6.0. At pH below 4.9 there is a single "acid form" (A), and a homogeneous "basic form" (B) exists at pH above 6.0. Between pH 5 and 6, BmPBP exists as a mixture of A and B in slow exchange on the NMR chemical shift time scale, with the transition midpoint at pH 5.4. The form B has a well-dispersed NMR spectrum, indicating that it represents a more structured, "closed" conformation than form A, which has a significantly narrower chemical shift dispersion. Conformational transitions of the kind observed here may explain heterogeneity reported for a variety of odorant-binding proteins, and it will be of interest to further investigate possible correlations with pH-dependent regulation of ligand binding and release in the biological function of this class of proteins.     

Identification and characterization of a glycosaminoglycan recognition element of the C chemokine lymphotactin

Chemokine-mediated recruitment of leukocytes in vivo depends on interactions with cell surface glycosaminoglycans. Lymphotactin, the unique member of the "C" chemokine subclass, is a highly basic protein that binds heparin, a glycosaminoglycan, with high affinity (approximately 10 nm). We detected lymphotactin-heparin binding by NMR and mapped this interaction to a narrow surface that wraps around the protein. Substitutions in and around this binding site and surface plasmon resonance analysis of heparin binding affinity identified two arginine residues of lymphotactin as critical for glycosaminoglycan binding. Both arginine mutant proteins and the combined double mutant had dramatically diminished in vivo activity in a leukocyte recruitment assay, suggesting that the lymphotactin-glycosaminoglycan interactions detected in vitro are important for the function of this chemokine. Our results demonstrate that like other chemokines, lymphotactin utilizes highly specific glycosaminoglycan-binding sites that represent potential targets for drug development.     

An approach to global fold determination using limited NMR data from larger proteins selectively protonated at specific residue types

Summary A combination of calculation and experiment is used to demonstrate that the global fold of larger proteins can be rapidly determined using limited NMR data. The approach involves a combination of heteronuclear triple resonance NMR experiments with protonation of selected residue types in an otherwise completely deuterated protein. This method of labelling produces proteins with -specific deuteration in the protonated residues, and the results suggest that this will improve the sensitivity of experiments involving correlation of side-chain (¹H and ¹³C) and backbone (¹H and ¹⁵N) amide resonances. It will allow the rapid assignment of backbone resonances with high sensitivity and the determination of a reasonable structural model of a protein based on limited NOE restraints, an application that is of increasing importance as data from the large number of genome sequencing projects accumulates. The method that we propose should also be of utility in extending the use of NMR spectroscopy to determine the structures of larger proteins.The first two authors contributed equally to this work.     

Nmrglue: an open source Python package for the analysis of multidimensional NMR data

Nmrglue, an open source Python package for working with multidimensional NMR data, is described. When used in combination with other Python scientific libraries, nmrglue provides a highly flexible and robust environment for spectral processing, analysis and visualization and includes a number of common utilities such as linear prediction, peak picking and lineshape fitting. The package also enables existing NMR software programs to be readily tied together, currently facilitating the reading, writing and conversion of data stored in Bruker, Agilent/Varian, NMRPipe, Sparky, SIMPSON, and Rowland NMR Toolkit file formats. In addition to standard applications, the versatility offered by nmrglue makes the package particularly suitable for tasks that include manipulating raw spectrometer data files, automated quantitative analysis of multidimensional NMR spectra with irregular lineshapes such as those frequently encountered in the context of biomacromolecular solid-state NMR, and rapid implementation and development of unconventional data processing methods such as covariance NMR and other non-Fourier approaches. Detailed documentation, install files and source code for nmrglue are freely available at http://nmrglue.com. The source code can be redistributed and modified under the New BSD license.     

Characterization of the carboxylate delivery module of transcarboxylase: following spontaneous decarboxylation of the 1.3S-CO2- subunit by NMR and FTIR spectroscopies

Transcarboxylase (TC) is a multisubunit enzyme that catalyzes the transfer of a carboxylate group from methylmalonyl-CoA (MMCoA) to pyruvate. The CO2- group is shuttled between the MMCoA and pyruvate binding sites by a biotin cofactor, covalently linked to the 1.3S subunit. Fully carboxylated 1.3S can be prepared in vitro using 1.3S, MMCoA, and catalytic amounts of the TC's MMCoA-binding subunit. The 1.3S-CO2- intermediate decarboxylates spontaneously over a period of hours, and this process was characterized by 1D and 2D NMR and FTIR spectroscopies. The NMR data yielded a first-order kinetic constant of 1.4 x 10(-3) min(-1) for the spontaneous decarboxylation. This rate was calculated from the 1D NMR spectrum by measuring the reappearance of biotin's ureido NH protons and the disappearance of peaks at 6.99 and 7.67 ppm assigned to Asn-8 and/or Asn-24 from the 1.3S's N-terminus. The latter peaks are absent in the 1D spectrum of non-carboxylated 1.3S due to exchange between two or more conformations within the N-terminus causing line broadening. It is proposed that interactions between the biotin-CO2- and the N-terminal amino acids perturb this conformational equilibrium causing some N-terminal residues to appear in the 1D NMR spectrum of the carboxylated form. Further details are apparent from a comparison of the 2D spectra of the 1.3S-CO2- and 1.3S proteins, where carboxylation causes several peaks from the C-terminal half to shift as well as the appearance of resonances due to some residues located at the N-terminal half of the protein. FTIR difference spectra were used also to follow spontaneous decarboxylation of the 1.3S-CO2-. For the carboxylated 1.3S, the difference spectra provided the vibrational signature of the CO2- on the biotin ring. A doublet was identified at 1695 and 1699 cm(-1) that increased in intensity with increasing t. This is assigned to an antisymmetric stretching vibration of the CO2- group bound to biotin on the 1.3S protein. Its position and profile provide further evidence for interactions occurring between the biotin-CO2- group and the 1.3S protein. These studies demonstrate the highly mobile, "poised" nature of the 1.3S protein engineered for its role as a CO2- translocator.     

Evaluation of Metabolomic Changes as a Biomarker of Chondrogenic Differentiation in 3D-cultured Human Mesenchymal Stem Cells Using Proton (1H) Nuclear Magnetic Resonance Spectroscopy

Purpose

The purpose of this study was to evaluate the metabolomic changes in 3D-cultured human mesenchymal stem cells (hMSCs) in alginate beads, so as to identify biomarkers during chondrogenesis using ¹H nuclear magnetic resonance (NMR) spectroscopy.

Materials and Methods

hMSCs (2×10⁶ cells/mL) were seeded into alginate beads, and chondrogenesis was allowed to progress for 15 days. NMR spectra of the chondrogenic hMSCs were obtained at 4, 7, 11, and 15 days using a 14.1-T (600-MHz) NMR with the water suppression sequence, zgpr. Real-Time polymerase chain reaction (PCR) was performed to confirm that that the hMSCs differentiated into chondrocytes and to analyze the metabolomic changes indicated by the NMR spectra.

Results

During chondrogenesis, changes were detected in several metabolomes as hMSC chondrogenesis biomarkers, e.g., fatty acids, alanine, glutamate, and phosphocholine. The metabolomic changes were compared with the Real-Time PCR results, and significant differences were determined using statistical analysis. We found that changes in metabolomes were closely related to biological reactions that occurred during the chondrogenesis of hMSCs.

Conclusions

In this study, we confirm that metabolomic changes detected by ¹H-NMR spectroscopy during chondrogenic differentiation of 3D-cultured hMSCs in alginate beads can be considered as biomarkers of stem cell differentiation.     

Measuring protein self-association using pulsed-field-gradient NMR spectroscopy: Application to myosin light chain 2

Summary At the millimolar concentrations required for structural studies, NMR spectra of the calcium-binding protein myosin light chain 2 (MLC2) showed resonance line widths indicative of extensive self-association. Pulsed-field-gradient (PFG) NMR spectroscopy was used to examine whether MLC2 aggregation could be prevented by the zwitterionic bile salt derivative 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS). PFG NMR measurements indicated that CHAPS was capable of preventing MLC2 self-association, but only at concentrations well above the critical micelle concentration of 7.5 mM. CHAPS was most effective at a concentration of 22.5 mM, where the apparent molecular mass of MLC2 correponded to a protein monomer plus seven molecules of bound detergent. The resolution and sensitivity of 2D ¹⁵N-¹H HSQC spectra of MLC2 were markedly improved by the addition of 25 mM CHAPS, consistent with a reduction in aggregation following addition of the detergent. The average amide nitrogen T₂ value for MLC2 increased from 30 ms in the absence of CHAPS to 56 ms in the presence of 25 mM CHAPS. The results of this study lead us to propose that PFG NMR spectroscopy can be used as a facile alternative to conventional techniques such as analytical ultracentrifugation for examining the self-association of biological macromolecules.     

Simulations of molecular diffusion in lattices of cells: insights for NMR of red blood cells

The pulsed field-gradient spin-echo (PGSE) nuclear magnetic resonance (NMR) experiment, conducted on a suspension of red blood cells (RBC) in a strong magnetic field yields a q-space plot consisting of a series of maxima and minima. This is mathematically analogous to a classical optical diffraction pattern. The method provides a noninvasive and novel means of characterizing cell suspensions that is sensitive to changes in cell shape and packing density. The positions of the features in a q-space plot characterize the rate of exchange across the membrane, cell dimensions, and packing density. A diffusion tensor, containing information regarding the diffusion anisotropy of the system, can also be derived from the PGSE NMR data. In this study, we carried out Monte Carlo simulations of diffusion in suspensions of "virtual" cells that had either biconcave disc (as in RBC) or oblate spheroid geometry. The simulations were performed in a PGSE NMR context thus enabling predictions of q-space and diffusion tensor data. The simulated data were compared with those from real PGSE NMR diffusion experiments on RBC suspensions that had a range of hematocrit values. Methods that facilitate the processing of q-space data were also developed.