Uncoiling collagen: a multidimensional mass spectrometry study.

Mass spectrometry can be used to determine structural information about ions by activating precursors and analysing the resulting series of fragments. Two-dimensional Fourier transform ion cyclotron resonance mass spectrometry (2D FT-ICR MS) is a technique that correlates the mass-to-charge (m/z) ratio of fragment and precursor ions in a single spectrum. 2D FT-ICR MS records the fragmentation of all ions in a sample without the need for isolation. To analyse specific precursors, horizontal cross-sections of the spectrum (fragment ion scans) are taken, providing an alternative to conventional tandem mass spectrometry (MS/MS) experiments. In this work, 2D FT-ICR MS has been used to study the tryptic digest of type I collagen, a large protein. Fragment ion scans have been extracted from the 2D FT-ICR MS spectrum for precursor m/z ratios: 951.81, 850.41, 634.34, and 659.34, and 2D FT-ICR MS spectra are compared with a set of 1D MS/MS spectra using different fragmentation methods. The results show that two-dimensional mass spectrometry excells at MS/MS of complex mixtures, simplifying spectra by eliminating contaminant peaks, and aiding the identification of species in the sample. Currently, with desktop computers, 2D FT-ICR MS is limited by data processing power, a limitation which should be alleviated using cluster parallel computing. In order to explore 2D FT-ICR MS for collagen, with reasonable computing time, the resolution in the fragment ion dimension is limited to 256k data points (compared to 4M data points in 1D MS/MS spectra), but the vertical precursor ion dimension has 4096 lines, so the total data set is 1G data points (4 Gbytes). The fragment ion coverage obtained with a blind, unoptimized 2D FT-ICR MS experiment was lower than conventional MS/MS, but MS/MS information is obtained for all ions in the sample regardless of selection and isolation. Finally, although all 2D FT-ICR MS peak assignments were made with the aid of 1D FT-ICR MS data, these results demonstrate the promise of 2D FT-ICR MS as a technique for studying complex protein digest mixtures.

Duplex formation and secondary structure of γ-PNA observed by NMR and CD.

Peptide nucleic acids (PNAs) are non-natural oligonucleotides mimics, wherein the phosphoribose backbone has been replaced by a peptidic moiety (N-(2-aminoethyl)glycine). This peptidic backbone lends itself to substitution and the γ-position has proven to yield oligomers with enhanced hybridization properties. In this study, we use Nuclear Magnetic Resonance (NMR) and Circular Dichroism (CD) to explore the properties of the supramolecular duplexes formed by these species. We show that standard Watson-Crick base pair as well as non-standard ones are formed in solution. The duplexes thus formed present marked melting transition temperatures substantially higher than their nucleic acid homologs. Moreover, the presence of a chiral group on the γ-peptidic backbone increases further this transition temperature, leading to very stable duplexes. PNA duplexes with a chiral backbone present a marked chiral secondary structure, observed by CD, and showing a common folding pattern for all studied structures. Nevertheless small differences are observed depending on the details of the nucleobase sequence.

¹9F nuclear magnetic resonance screening of glucokinase activators.

Human hexokinase enzyme IV (EC 2.7.1.1) catalyzes the phosphorylation of glucose and regulates the level of glucose. This enzyme exhibits strong positive cooperativity due to an allosteric transition between an inactive form and a closed active form. This form can be stabilized by activators and, thus, can increase its turnover by a kinetic memory effect characterized by a slow decay to the inactive state. The structural details of this kinetic allostery are known. Several synthetic activators have been reported. We present a preliminary nuclear magnetic resonance (NMR) screening of a chemical library in search of molecules with some affinity for glucokinase (GK). The library, composed of eight molecules with known activity as well as molecules that display no interaction, has been tested using the FAXS (fluorine chemical shift anisotropy and exchange for screening) method, based on monitoring the R2 relaxation of the (19)F spin. To ensure a valid interaction measurement, the enzyme was placed in the presence of glucose and magnesium. The binding signal of one known fluorinated ligand was measured by determining the displacement of the known ligand. This simple measure of the (19)F signal intensity after an 80-ms spin echo correlates nicely with the EC50, opening a route for NMR screening of GK activators.

A new NMR technique to probe protein-ligand interaction.

Non covalent grafting of proteins on affinity phases is a very common approach for isolation, purification and re-concentration of tagged proteins. Many biophysical studies are conducted on these grafted proteins (surface plasmon resonance, quartz crystal microbalance, etc.) showing that the integrity and function of the protein is usually maintained. However, NMR studies of such samples were not undertaken so far, due to the broadening observed on this kind of heterogeneous samples. We present here the use of the HR-MAS technology to obtain 2D NMR spectra of the MAGI-1 PDZ2/6 protein domain, C13-labeled, tagged with a His-tag and grafted on a Nickel affinity resin. We optimized the C13 Methyl SOFAST HMQC experiment allowing important gains in terms of signal-to-noise. The gain comes from the gathering of proton magnetization from the resin material to the protein under study. Several methyl signals from the unstructured C-terminal tail, which is involved in the binding of the PDZ domain to C-terminal peptides of its partners, were observed and measured. The interaction of the bound PDZ domain with cognate peptides was monitored using <500µg of protein sample. A response proportional to the peptide Kd is obtained, indicating that the method can be used to rapidly and efficiently monitor protein-ligand interactions.

A general algorithm for peak-tracking in multi-dimensional NMR experiments.

We present an algorithmic method allowing automatic tracking of NMR peaks in a series of spectra. It consists in a two phase analysis. The first phase is a local modeling of the peak displacement between two consecutive experiments using distance matrices. Then, from the coefficients of these matrices, a value graph containing the a priori set of possible paths used by these peaks is generated. On this set, the minimization under constraint of the target function by a heuristic approach provides a solution to the peak-tracking problem. This approach has been named GAPT, standing for General Algorithm for NMR Peak Tracking. It has been validated in numerous simulations resembling those encountered in NMR spectroscopy. We show the robustness and limits of the method for situations with many peak-picking errors, and presenting a high local density of peaks. It is then applied to the case of a temperature study of the NMR spectrum of the Lipid Transfer Protein (LTP).

NMRb: a web-site repository for raw NMR datasets.

UNLABELLED: The development of NMR in structural proteomics requires the availability of automatic structure determination methods. Many researchers are commonly confronted with the lack of raw datasets during the validation step of such methods. In order to increase test possibilities, the NMRb web-site offers a database of NMR raw datasets, ordered by spectral characteristics. AVAILABILITY: NMRb is available from: http://nmrb.cbs.cnrs.fr. SUPPLEMENTARY INFORMATION: General organization of NMRb figure, relational model organization, and XML structure files are available from http://nmrb.cbs.cnrs.fr/nmrb-doc.html.

The use of sample rotation for minimizing convection effects in self-diffusion NMR measurements.

Undesirable temperature gradients in a NMR sample tube are usually generated by an inappropriate temperature regulation system. We have shown that such convection effects can greatly distort the measurement of translational self-diffusion coefficients. The use of sample spinning helps to minimize such undesirable effects by disruption of convection fluxes due to resulting Coriolis forces that have a strongly stabilizing effect on the conducting state of the system (J. Lounila et al., J. Magn. Reson. A 118, 50 (1996)). This simple trick allows the accurate measurement of diffusion coefficients for a wide range of temperatures and solvents without the need for a convection-compensated NMR pulse sequences or more sophisticated temperature control units. Experimental data obtained for some target compounds dissolved in several common deuterated solvents at different temperatures are reported and discussed.

Direct recognition of foreign MHC determinants by naive T cells mobilizes specific Vbeta families without skewing of the complementarity-determining region 3 length distribution.

The capacity of T cells to interact with nonself-APC, also referred to as direct allorecognition, is an essential feature of the cellular response involved in graft rejection. However, there is no study on TCR repertoire biases associated with direct restricted T cell activation. In this paper, we have addressed the impact of direct recognition on the whole naive T cell repertoire, using a new approach that provides, for the first time, an integrated depiction of the quantitative and qualitative alterations in the TCR Vbeta transcriptome. This method can differentiate resting patterns from polyclonally activated ones, as evidenced by superantigen usage. According to this new readout, we show that direct recognition of nonself-MHC molecules triggers mRNA accumulation of several TCR Vbeta families, specific to the combination studied. Moreover, in marked contrast to the situation that prevails in indirect allorecognition, T cell activation through the direct presentation pathway was not associated with skewing of the complementarity determining region (CDR) 3 length distribution. Altogether, these data argue for the significance of TCR contacts with the MHC framework in direct allorecognition. In addition, the TCR diversity mobilized by this interaction and the massive TCRbeta mRNA accumulation observed after a few days of culture suggest that a significant proportion of naive T cells receive a signal leading to TCRbeta transcriptional activation even though only a few of them engage in mitosis.

NMR monitoring of accumulation and folding of 15N-labeled protein overexpressed in Pichia pastoris.

Postgenomic studies have led to an increasing demand for isotope-labeled proteins. We present a method for producing large quantities of truly native (15)N-labeled protein. Based on the secretion capabilities of the yeast Pichia pastoris, the recombinant protein is easily purified in a single step as it is secreted. Control of all nitrogen sources permits very high labeling yields. As a result, accumulation and folding of the recombinant protein can be monitored by heteronuclear NMR without purification. Comparison of sample spectra with the spectrum of the purified recombinant protein allows detection of the secreted protein in the culture and monitoring of its folding, from the start of the induction phase. The detection limit for a (15)N-labeled protein is estimated as 20 microM and corresponds, for a 10-kDa protein, to a load of 40 mg/liter in the fermentor. This concentration is reached by most reported preparations in P. pastoris. Further concentration by ultrafiltration would compensate for lower production. This procedure may be useful in many structural genomics and combinatorial chemistry screening projects where most protein productions meet the requirements for this method.

Characterization and molecular basis of the oligomeric structure of HIV-1 nef protein.

The Nef protein of human immunodeficiency virus type I (HIV-1) is an important determinant for the onset of AIDS disease. The self-association properties of HIV-1 Nef are analyzed by chemical cross-linking, dynamic light scattering, equilibrium analytical ultracentrifugation, and NMR spectroscopy. The experimental data show that the HIV-1 Nef core domain forms stable homo-dimers and trimers in solution, but not higher oligomers. These Nef homomers are not covalently linked by disulfide bridges, and the equilibrium between these forms is dependent on the Nef concentration. We further provide the molecular basis for the Nef core dimers and trimers obtained by analysis of crystallographic models. Oligomerization of biological polypeptides is a common tool used to trigger events in cellular signaling and endocytosis, both of which are targeted by Nef. The quaternary structure of Nef may be of physiological importance and may help to connect its cellular targets or to increase affinity of the viral molecule for its ligands. The herein described models for Nef dimers and trimers will allow further mutational studies to elucidate their role in vivo. These results provide novel insight into the structural and functional relationships of this important viral protein. Moreover, the oligomer interface may represent a novel target for the design of antiviral agents.

RESCUE: an artificial neural network tool for the NMR spectral assignment of proteins.

The assignment of the 1H spectrum of a protein or a polypeptide is the prerequisite for advanced NMR studies. We present here an assignment tool based on the artificial neural network technology, which determines the type of the amino acid from the chemical shift values observed in the 1H spectrum. Two artificial neural networks have been trained and extensively tested against a non-redundant subset of the BMRB chemical shift data bank [Seavey, B.R. et al. (1991) J. Biomol. NMR, 1, 217-236]. The most promising of the two accomplishes the analysis in two steps, grouping related amino acids together. It presents a mean rate of success above 80% on the test set. The second network tested separates down to the single amino acid; it presents a mean rate of success of 63%. This tool has been used to assist the manual assignment of peptides and proteins and can also be used as a block in an automated approach to assignment. The program has been called RESCUE and is made publicly available at the following URL: http:(/)/www.infobiosud.univ-montp1.fr/rescue.

Quantitative measurement of longitudinal and transverse cross-relaxation rates: an application to the analysis of the internal dynamics of ranalexin in water and trifluoroethanol.

We describe a quantitative processing method which gives access to the longitudinal and transverse cross-relaxation rates from off-resonance ROESY intensities. This method takes advantage of the dependence of the off-resonance ROESY experiments at any mixing time and any spin-lock angle θ on two relaxation matrices, the longitudinal and the transverse ones. This allows one to take into account multistep magnetization transfers even if the measurements are performed only at one or two mixing times. The ratio of the longitudinal to transverse cross-relaxation rates can then be used as a local indicator of the internal dynamics, without assuming a structure or a model of motion. After validation of this processing method by numerical simulations, it is applied to the analysis of the dynamics of the peptide ranalexin dissolved in pure water and in water/TFE.

Refined solution structure and backbone dynamics of 15N-labeled C12A-p8MTCP1 studied by NMR relaxation.

MTCP1 (for Mature-T-Cell Proliferation) was the first gene unequivocally identified in the group of uncommon leukemias with a mature phenotype. The three-dimensional solution structure of the human p8MTCP1 protein encoded by the MTCP1 oncogene has been previously determined by homonuclear proton two-dimensional NMR methods at 600 MHz: it consists of an original scaffold comprising three alpha-helices, associated with a new cysteine motif. Two of the helices are covalently paired by two disulfide bridges, forming an alpha-hairpin which resembles an antiparallel coiled-coil. The third helix is orientated roughly parallel to the plane defined by the alpha-antiparallel motif and appears less well defined. In order to gain more insight into the details of this new scaffold, we uniformly labeled with nitrogen-15 a mutant of this protein (C12A-p8MTCP1) in which the unbound cysteine at position 12 has been replaced by an alanine residue, thus allowing reproducibly high yields of recombinant protein. The refined structure benefits from 211 additional NOEs, extracted from 15N-edited 3D experiments, and from a nearly complete set of phi angular restraints allowing the estimation of the helical content of the structured part of the protein. Moreover, measurements of 15N spin relaxation times and heteronuclear 15N¿1H¿NOEs provided additional insights into the dynamics of the protein backbone. The analysis of the linear correlation between J(0) and J(omega) was used to interpret relaxation parameters. It appears that the apparent relative disorder seen in helix III is not simply due to a lack of experimental constraints, but associated with substantial contributions of sub-nanosecond motions in this segment.

An NMR assignment module implemented in the Gifa NMR processing program.

MOTIVATION: Peptide and protein structures are determined daily using NMR spectroscopy. Assignment of the NMR spectra is an important step within the procedure and is usually the limiting one. Computer-aided assignment tools should be user friendly with open architecture to communicate with other programs involved in the structure determination. RESULTS: Here we present an interactive NMR assignment module which provides numerous graphic tools for the user. The module is composed of a database management system-handling peaks, spins and spin-systems. The assignment information is maintained as a set of interrelated associative arrays, which serve as generic high-level data structures. The module is developed in the macro language embedded in the Gifa NMR processing program (Pons et al. , J. Biomol. NMR, 8 , 445-452, 1996). This provides the user with a consistent interface, a set of sophisticated tools, and an easily extendible and customizable environment. AVAILABILITY: The program is available on request from the authors. The Gifa package can be accessed at: ((http://www.cbs. univ-montp1.fr/GIFA)) CONTACT: Marc-Andre.Delsuc@cbs.univ-montp1.fr

Structural and functional roles of a conserved proline residue in the alpha2 helix of Escherichia coli thioredoxin.

Proline 40 in Escherichia coli thioredoxin is located close to the redox active site (Cys32-Cys35) within the alpha2 helix. The conservation of this residue among most of the thioredoxins suggests that it could play an important role in the structure and/or function of this protein. We have substituted Pro40 for Ala by using site-directed mutagenesis and expressed the mutant P40A in E.coli. The effects of the mutation on the biophysical and biological properties of thioredoxin have been analyzed and compared with molecular dynamics simulations. Modeling predicted that the replacement of Pro40 by Ala induced a displacement of the active site which exposes Trp31 to the solvent and opens a cleft located between helices alpha2 and alpha3. The solvation free energy (SFE) calculation also indicated that P40A became more hydrophobic as W31 became more accessible. These predictions were totally in agreement with the experimental results. The mutant P40A exhibited chromatographic behavior and fluorescence properties very different from those of the wild-type (WT) protein, in relationship with the displacement of W31. The determination of the free energy of unfolding of P40A showed that the mutant was globally destabilized by 2.9 kcal/mol. However, the effect of the mutation on the transition curve was highly unusual as the midpoint of the unfolding transition increased, indicating that some local structures were actually stabilized by the mutation. Despite these structural modifications, neither the ability of the protein to reduce a chloroplastic enzyme nor its reactivity with the bacterial reductase decreased. The only functional difference was the higher stability of P40A in light activation of NADP-malate dehydrogenase under air, which suggests that the mutant was less rapidly re-oxidized than WT. Therefore, it can be concluded that Pro40 is not essential for maintaining the redox function of thioredoxin but rather is required for the stability of the protein.

Gifa V. 4: A complete package for NMR data set processing.

The Gifa program is designed for processing, displaying and analysing 1D, 2D and 3D NMR data sets. It has been constructed in a modular fashion, based on three independent modules: a set of commands that perform all the basic processing operations such as apodisation functions, a complete set of Fourier Transforms, phasing and baseline correction, peak-picking and line fitting, linear prediction and maximum entropy processing; a set of command language primitives that permit the execution of complex macro commands; and a set of graphic commands that permit to build a complete graphic user interface, allowing the user to interact easily with the program. We have tried to create a versatile program that can be easily extended according to the user's requirements and that is adapted to a novice as well as an experienced user. The program runs on any UNIX computer, with or without graphic display, in interactive or batch mode.

Use of the Cadzow procedure in 2D NMR for the reduction of t(1) noise.

A data processing approach is proposed for reducing the t(1) noise observed in multidimensional NMR spectra. This method is based on the use of the Cadzow procedure [Cadzow, J.A. (1988) IEEE Trans. Acous. Speech Signal Proc., 36, 49-62], and is demonstrated to be efficient for simulated cases as well as real experiments.

An estimate of spin diffusion in a spin subset: Application to iterative distance calculation from 3D (15)N NOESY-HMQC.

A method for quantification of distances between amide hydrogens using only the 3D NOESY-HMQC experiment recorded on a (15)N-labelled protein is presented. This method is based on an approximate expression of the NOE intensities between amide hydrogens obtained from continuum modelling of the non-amide spins; this expression is used in a distance calculation algorithm. The algorithm has been named CROWD, standing for Continuum approximation of Relaxati On path Ways between Dilute spins. This approximation as well as the CROWD algorithm are tested on a simulated case; the CROWD algorithm is then applied to experimental data, measured on a fragment of bacteriorhodopsin.

The DNA-binding domain of the yeast Saccharomyces cerevisiae CYP1(HAP1) transcription factor possesses two zinc ions which are complexed in a zinc cluster.

Various fragments of the N-terminal, DNA-binding domain of the yeast Saccharomyces cerevisiae transcriptional activator CYP1(HAP1) have been cloned and expressed in Escherichia coli. The corresponding polypeptides have been analysed biochemically and we have undertaken a more extensive physical study of a fragment consisting of amino acids 49-126 [CYP1(49-126)]. We show that this CYP1(49-126) peptide requires zinc or cadmium in the growth medium in order to maintain a stable structure. A method to purify CYP1(49-126) is presented. We demonstrate that the purified CYP1(49-126) fragment contains two zinc ions/fragment or two cadmium ions/fragment, which are necessary for DNA binding. 113Cd one-dimensional NMR data suggest that CYP1(HAP1) has a tetrahedral coordination, and that it forms a zinc-cluster complex like GAL4.

Accurate estimation of inter-atomic distances in large proteins by NMR.

Recently a method was proposed which permits the extraction of the exact interatomic distance information from the measurement of the evolution of a single cross-peak relative to the mixing time in a NOESY experiment. This is performed through a careful multi-exponential analysis allowing the extraction of the relaxation parameter, and, consequently, the inter-proton distance. We investigate in the present paper whether this technique, already evaluated theoretically, can be used in a real experimental case. We have recorded and analyzed a set of 56 NOESY experiments on a lysozyme sample. Some 81 nOe build-up curves obtained from these data were analyzed in terms of distance. It is shown that the correlation between the measured distances and the reference distances obtained from crystallographic studies, is quite good. An accuracy of the order of 10% is obtained.