Hyperdimensional NMR spectroscopy with nonlinear sampling

An approach is described for joint interleaved recording, real-time processing, and analysis of NMR data sets. The method employs multidimensional decomposition to find common information in a set of conventional triple-resonance spectra recorded in the nonlinear sampling mode, and builds a model of hyperdimensional (HD) spectrum. While preserving sensitivity per unit of measurement time and allowing for maximal spectral resolution, the approach reduces data collection time on average by 2 orders of magnitude compared to the conventional method. The 7-10 dimensional HD spectrum, which is represented as a set of deconvoluted 1D vectors, is easy to handle and amenable for automated analysis. The method is exemplified by automated assignment for two protein systems of low and high spectral complexity: ubiquitin (globular, 8 kDa) and zetacyt (naturally disordered, 13 kDa). The collection and backbone assignment of the data sets are achieved in real time after approximately 1 and 10 h, respectively. The approach removes the most critical time bottlenecks in data acquisition and analysis. Thus, it can significantly increase the value of NMR spectroscopy in structural biology, for example, in high-throughput structural genomics applications.

Targeted acquisition for real-time NMR spectroscopy

A target-oriented approach for the acquisition of information in biomolecular NMR spectroscopy is being developed. This approach combines concurrent data accumulation, processing, and monitoring of spectral quality. Real-time estimation of parameters allows acquisition to be stopped when results are complete and have a specified precision. The technique is based on multidimensional decomposition, which can process incomplete data. An incremental nonuniform sampling scheme ensures the optimization of resolution sensitivity. To validate this method, 3D HNCO spectra of three biomolecular systems (8 kDa ubiquitin, 22 kDa barstar-barnase complex, and 82 kDa malate synthase G) are processed incrementally at small acquisition time steps. The range of molecular sizes illustrates applicability in both sample- and sensitivity-limited regimes. In each case, the target was to acquire all backbone resonances in the spectra. For the three systems, the targets are achieved after 4.5 min, 1.6 h, and 22 h of acquisition time, respectively. A number of other targets that can be similarly monitored as a function of time are discussed.

Quantification of H/D isotope effects on protein hydrogen-bonds by h3JNC' and 1JNC' couplings and peptide group 15N and 13C' chemical shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.

From structure and dynamics of protein L7/L12 to molecular switching in ribosome

Based on the (1)H-(15)N NMR spectroscopy data, the three-dimensional structure and internal dynamic properties of ribosomal protein L7 from Escherichia coli were derived. The structure of L7 dimer in solution can be described as a set of three distinct domains, tumbling rather independently and linked via flexible hinge regions. The dimeric N-terminal domain (residues 1-32) consists of two antiparallel alpha-alpha-hairpins forming a symmetrical four-helical bundle, whereas the two identical C-terminal domains (residues 52-120) adopt a compact alpha/beta-fold. There is an indirect evidence of the existence of transitory helical structures at least in the first part (residues 33-43) of the hinge region. Combining structural data for the ribosomal protein L7/L12 from NMR spectroscopy and x-ray crystallography, it was suggested that its hinge region acts as a molecular switch, initiating "ratchet-like" motions of the L7/L12 stalk with respect to the ribosomal surface in response to elongation factor binding and GTP hydrolysis. This hypothesis allows an explanation of events observed during the translation cycle and provides useful insights into the role of protein L7/L12 in the functioning of the ribosome.

Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.

An intrahelical salt bridge within the trigger site stabilizes the GCN4 leucine zipper

We previously reported that a helical trigger segment within the GCN4 leucine zipper monomer is indispensable for the formation of its parallel two-stranded coiled coil. Here, we demonstrate that the intrinsic secondary structure of the trigger site is largely stabilized by an intrahelical salt bridge. Removal of this surface salt bridge by a single amino acid mutation induced only minor changes in the backbone structure of the GCN4 leucine zipper dimer as verified by nuclear magnetic resonance. The mutation, however, substantially destabilized the dimeric structure. These findings support the proposed hierarchic folding mechanism of the GCN4 coiled coil in which local helix formation within the trigger segment precedes dimerization.

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.

NMR evidence for progressive stabilization of native-like structure upon aggregation of acid-denatured LysN

The acid-denatured form of the protein LysN aggregates reversibly at pH 2.0. The strength of self-association increases with increasing Cl(-) anion concentration. At low concentrations of protein or Cl(-) anion, resonances of denatured LysN are in slow exchange with a minor form of the protein, which shows native-like NMR chemical shifts. The minor native-like resonances increase in intensity with increasing protein concentration, demonstrating that a native-like monomer fold is stabilized on aggregation of the acid-denatured protein. At high concentrations of protein or Cl(-) anion, interconversion between the major and minor resonances appears to shift from slow to intermediate exchange on the NMR timescale. NMR line-broadening is more pronounced for the major resonances of the denatured protein, which show sigmoidal decay curves with increasing Cl(-) concentration. The mid-points of the decay curves for residues in different parts of the molecule are non-coincident. We propose that differences in the NMR line-broadening transitions of individual residues reflect a stepwise stabilization of native-like structure on aggregation, starting with the segments of the protein that form the initial association interface. The resonances of the denatured protein with the greatest sensitivity to self-association correspond roughly to those that are most perturbed in the native protein on binding of the natural substrate tRNA(Lys). This suggests that the hydrophobic surfaces that promote intermolecular misfolding of acid-denatured LysN, may resemble those used for substrate binding by the native protein.

NMR hydrogen exchange of the OB-fold protein LysN as a function of denaturant: the most conserved elements of structure are the most stable to unfolding

The structure of LysN contains an OB-fold motif composed of a structurally conserved five-stranded beta-barrel capped by a poorly conserved alpha-helix between strands beta3 and beta4. Two additional alpha-helices, unique to the LysN structure, flank the N terminus of the OB-fold. The stability of LysN to unfolding has been investigated with NMR native state hydrogen exchange measurements as a function of guanidinium hydrochloride concentration, and equilibrium unfolding transitions monitored by ellipticity at 222 nm and fluorescence at 350 nm. The spectrophotometric measurements suggest an apparent two-state unfolding transition with DeltaGu(0) approximately 6 kcal/mol and m approximately 3 kcal/(molM). By contrast, NMR hydrogen exchange measurements manifest a distribution of DeltaGu(0) and m values which indicate that the protein can undergo subglobal unfolding. The largest DeltaGu(0) values from hydrogen exchange are for residues in the beta-sheet of the protein. These values, which reflect complete unfolding of the protein, are between 3 and 4 kcal/mol higher than those obtained from circular dichroism or fluorescence. This discrepancy may be due to the comparison of NMR hydrogen exchange parameters measured at residue-level resolution, with spectrophotometric parameters that reflect an unresolved super position of unfolding transitions of the alpha-helices and beta-strands. The largest DeltaGu(0) values obtained from hydrogen exchange for the subset of residues in the alpha-helices of the protein, agree with the DeltaGu(0) values obtained from circular dichroism or fluorescence. Based on the hydrogen exchange data, however, the three alpha-helices of LysN are on average 3 kcal/mol less stable than the beta-sheet. Consistent with the subglobal unfolding of LysN evinced by hydrogen exchange, a deletion mutant that lacks the first alpha-helix of the protein retains a cooperatively folded structure. Taken together with previous results on the OB-fold proteins SN and CspA, the present results for LysN suggest that the most conserved elements of structure in the OB-fold motif are the most resistant to denaturation. In all three proteins, stability to denaturation correlates with sequence hydrophobicity.

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.

Three-dimensional structure of toxin OSK1 from Orthochirus scrobiculosus scorpion venom

A 600 MHz 1H NMR study of toxin OSK1, blocker of small-conductance Ca2+-activated K+ channels, is presented. The unambiguous sequential assignment of all the protons of the toxin was obtained using TOCSY, DQF-COSY, and NOESY experiments at pH 3.0 (10, 30, and 45 degrees C) in aqueous solution. 3J(N alpha), 3J(alphabeta) vicinal spin coupling constants were determined in high-resolution spectra. The cross-peak volumes in NOESY spectra and the coupling constants were used to define the local structure of the protein by the program HABAS and to generate torsion angle and interproton distance constraints for the program DIANA. Hydrogen-deuterium exchange rates of amide protons showed possible locations of hydrogen bonds. The hydrogen bond acceptors and disulfide bridges between residues 8-28, 14-33, and 18-35 were determined when analyzing distance distribution in preliminary DIANA structures. All constraints were used to obtain a set of 30 structures by DIANA. The resulting rms deviations over 30 structures are 1.30 A for the heavy atoms and 0.42 A for the backbone heavy atoms. The structures were refined by constrained energy minimization using the SYBYL program. Their analysis indicated the existence of the alpha-helix (residues 10-21) slightly distorted at the Cys14 residue, two main strands of the antiparallel beta-sheet (24-29, 32-38), and the extended fragment (2-6). The motif is stabilized by the disulfide bridges in the way, common to all known scorpion toxins. Using the fine spatial toxin structure, alignment of the homologues, mutagenesis analysis, and comparison of scorpion toxin family functions, we delineate some differences significant for the toxin specificity.

Amino acid residue: is it structural or functional?

A new approach is suggested for delineating the structural and functional amino acid residues in proteins with known three-dimensional structure, basing on the involvement of residues in intramolecular hydrophobic and hydrophilic interactions and additional information about the conservativity of the residues. The approach is applied to the families of homologous neurotoxins and cardiotoxins. The results obtained concerning the role of amino acid residues in both families of toxins accord well with the similarity of their fold, but different mechanisms of action. Current approach can be used for detailed characterization of protein spatial structures, as well as for rational protein engineering.