A New Amide Proton R 1ρ Experiment Permits Accurate Characterization of Microsecond Time-scale Conformational Exchange

A new off-resonance spin-lock experiment to record relaxation dispersion profiles of amide protons is presented. The sensitivity-enhanced HSQC-type sequence is designed to minimize the interference from cross-relaxation effects and ensure that the dispersion profiles in the absence of μs-ms time-scale dynamics are flat. Toward this end (i) the proton background is eliminated by sample deuteration (Ishima et al., 1998), (ii) ¹H spin lock is applied to two-spin modes 2(H_xSin ϑ + H_zCos ϑ) N_z, and (iii) the tilt angle ϑ ≈ 35° is maintained throughout the series of measurements (Desvaux et al. Mol. Phys., 86 (1995) 1059). The relaxation dispersion profiles recorded in this manner sample a wide range of effective rf field strengths (up to and in excess of 20 kHz) which makes them particularly suitable for studies of motions on the time scale ≤100 μs. The new experiment has been tested on the Ca²⁺-loaded regulatory domain of cardiac troponin C. Many residues show pronounced dispersions with remarkably similar correlation times of 30 μs. Furthermore, these residues are localized in the regions that have been previously implicated in conformational changes (Spyracopoulos et al. Biochemistry, 36 (1997) 12138) Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Very large residual dipolar couplings from deuterated ubiquitin

Main-chain (1)H(N)-(15)N residual dipolar couplings (RDCs) ranging from approximately -200 to 200?Hz have been measured for ubiquitin under strong alignment conditions in Pf1 phage. This represents a ten-fold increase in the degree of alignment over the typical weakly aligned samples. The measurements are made possible by extensive proton-dilution of the sample, achieved by deuteration of the protein with partial back-substitution of labile protons from 25?% H(2)O / 75?% D(2)O buffer. The spectral quality is further improved by application of deuterium decoupling. Since standard experiments using fixed-delay INEPT elements cannot accommodate a broad range of couplings, the measurements were conducted using J-resolved and J-modulated versions of the HSQC and TROSY sequences. Due to unusually large variations in dipolar couplings, the trosy (sharp) and anti-trosy (broad) signals are often found to be interchanged in the TROSY spectra. To distinguish between the two, we have relied on their respective (15)N linewidths. This strategy ultimately allowed us to determine the signs of RDCs. The fitting of the measured RDC values to the crystallographic coordinates of ubiquitin yields the quality factor Q?=?0.16, which confirms the perturbation-free character of the Pf1 alignment. Our results demonstrate that RDC data can be successfully acquired not only in dilute liquid crystals, but also in more concentrated ones. As a general rule, the increase in liquid crystal concentration improves the stability of alignment media and makes them more tolerant to variations in sample conditions. The technical ability to measure RDCs under moderately strong alignment conditions may open the door for development of alternative alignment media, including new types of media that mimic biologically relevant systems.     

Paramagnetic relaxation enhancements in unfolded proteins: Theory and application to drkN SH3 domain

Site‐directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as ¹H^N is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean‐force potential (Smoluchowski equation); (iv) explicit‐atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the ¹H^N residue number and static magnetic field strength; the results are appreciably different from the Gillespie–Shortle model. At the same time, the Gillespie–Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between ¹H^N spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin‐labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme—using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates—lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random‐coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random‐coil protein in good solvent/poor solvent.     

Structural determinants for high-affinity binding in a Nedd4 WW3* domain-Comm PY motif complex

Interactions between the WW domains of Drosophila Nedd4 (dNedd4) and Commissureless (Comm) PY motifs promote axon crossing at the CNS midline and muscle synaptogenesis. Here we report the solution structure of the dNedd4 WW3* domain complexed to the second PY motif (227'TGLPSYDEALH237') of Comm. Unexpectedly, there are interactions between WW3* and ligand residues both N- and C-terminal to the PY motif. Residues Y232'-L236' form a helical turn, following the PPII helical PY motif. Mutagenesis and binding studies confirm the importance of these extensive contacts, not simultaneously observed in other WW domain complexes, and identify a variable loop in WW3* responsible for its high-affinity interaction. These studies expand our general understanding of the molecular determinants involved in WW domain-ligand recognition. In addition, they provide insights into the specific regulation of dNedd4-mediated ubiquitination of Comm and subsequent internalization of Comm or the Comm/Roundabout complex, critical for CNS and muscle development.     

15NH/D-SOLEXSY experiment for accurate measurement of amide solvent exchange rates: application to denatured drkN SH3

Amide solvent exchange rates are regarded as a valuable source of information on structure/dynamics of unfolded (disordered) proteins. Proton-based saturation transfer experiments, normally used to measure solvent exchange, are known to meet some serious difficulties. The problems mainly arise from the need to (1) manipulate water magnetization and (2) discriminate between multiple magnetization transfer pathways that occur within the proton pool. Some of these issues are specific to unfolded proteins. For example, the compensation scheme used to cancel the Overhauser effect in the popular CLEANEX experiment is not designed for use with unfolded proteins. In this report we describe an alternative experimental strategy, where amide ¹⁵N is used as a probe of solvent exchange. The experiment is performed in 50% H₂O–50% D₂O solvent and is based on the (HACACO)NH pulse sequence. The resulting spectral map is fully equivalent to the conventional HSQC. To fulfill its purpose, the experiment monitors the conversion of deuterated species, ¹⁵N^D, into protonated species, ¹⁵N^H, as effected by the solvent exchange. Conceptually, this experiment is similar to EXSY which prompted the name of ¹⁵N^H/D-SOLEXSY (SOLvent EXchange SpectroscopY). Of note, our experimental scheme, which relies on nitrogen rather than proton to monitor solvent exchange, is free of the complications described above. The developed pulse sequence was used to measure solvent exchange rates in the chemically denatured state of the drkN SH3 domain. The results were found to correlate well with the CLEANEX-PM data, r = 0.97, thus providing a measure of validation for both techniques. When the experimentally measured exchange rates are converted into protection factors, most of the values fall in the range 0.5–2, consistent with random-coil behavior. However, elevated values, ca. 5, are obtained for residues R38 and A39, as well as the side-chain indole of W36. This is surprising, given that high protection factors imply hydrogen bonding or hydrophobic burial not expected to occur in a chemically denatured state of a protein. We, therefore, hypothesized that elevated protection factors are an artefact arising from the calculation of the reference (random-coil) exchange rates. To confirm this hypothesis, we prepared samples of several short peptides derived from the sequence of the drkN SH3 domain; these samples were used to directly measure the reference exchange rates. The revised protection factors obtained in this manner proved to be close to 1.0. These results also have implications for the more compact unfolded state of drkN SH3, which appears to be fully permeable to water as well, with no manifestations of hydrophobic burial.     

Orienting domains in proteins using dipolar couplings measured by liquid-state NMR: differences in solution and crystal forms of maltodextrin binding protein loaded with beta-cyclodextrin

Protein function is often regulated by conformational changes that occur in response to ligand binding or covalent modification such as phosphorylation. In many multidomain proteins these conformational changes involve reorientation of domains within the protein. Although X-ray crystallography can be used to determine the relative orientation of domains, the crystal-state conformation can reflect the effect of crystal packing forces and therefore may differ from the physiologically relevant form existing in solution. Here we demonstrate that the solution-state conformation of a multidomain protein can be obtained from its X-ray structure using an extensive set of dipolar couplings measured by triple-resonance multidimensional NMR spectroscopy in weakly aligning solvent. The solution-state conformation of the 370-residue maltodextrin-binding protein (MBP) loaded with beta-cyclodextrin has been determined on the basis of one-bond (15)N-H(N), (15)N-(13)C', (13)C(alpha)-(13)C', two-bond (13)C'-H(N), and three-bond (13)C(alpha)-H(N) dipolar couplings measured for 280, 262, 276, 262, and 276 residues, respectively. This conformation was generated by applying hinge rotations to various X-ray structures of MBP seeking to minimize the difference between the experimentally measured and calculated dipolar couplings. Consistent structures have been derived in this manner starting from four different crystal forms of MBP. The analysis has revealed substantial differences between the resulting solution-state conformation and its crystal-state counterpart (Protein Data Bank accession code 1DMB) with the solution structure characterized by an 11(+/-1) degrees domain closure. We have demonstrated that the precision achieved in these analyses is most likely limited by small uncertainties in the intradomain structure of the protein (ca 5 degrees uncertainty in orientation of internuclear vectors within domains). In addition, potential effects of interdomain motion have been considered using a number of different models and it was found that the structures derived on the basis of dipolar couplings accurately represent the effective average conformation of the protein.     

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.     

Regulation of phenylacetic acid degradation genes of Burkholderia cenocepacia K56-2

Background  

Metabolically versatile soil bacteria Burkholderia cepacia complex (Bcc) have emerged as opportunistic pathogens, especially of cystic fibrosis (CF). Previously, we initiated the characterization of the phenylacetic acid (PA) degradation pathway in B. cenocepacia, a member of the Bcc, and demonstrated the necessity of a functional PA catabolic pathway for full virulence in Caenorhabditis elegans. In this study, we aimed to characterize regulatory elements and nutritional requirements that control the PA catabolic genes in B. cenocepacia K56-2.     

Detection of nanosecond time scale side-chain jumps in a protein dissolved in water/glycerol solvent

In solution, the correlation time of the overall protein tumbling, τ _R , plays a role of a natural dynamics cutoff—internal motions with correlation times on the order of τ _R or longer cannot be reliably identified on the basis of spin relaxation data. It has been proposed some time ago that the ‘observation window’ of solution experiments can be expanded by changing the viscosity of solvent to raise the value of τ _R . To further explore this concept, we prepared a series of samples of α-spectrin SH3 domain in solvent with increasing concentration of glycerol. In addition to the conventional ¹⁵N labeling, the protein was labeled in the Val, Leu methyl positions (¹³CHD₂ on a deuterated background). The collected relaxation data were used in asymmetric fashion: backbone ¹⁵N relaxation rates were used to determine τ _R across the series of samples, while methyl ¹³C data were used to probe local dynamics (side-chain motions). In interpreting the results, it has been initially suggested that addition of glycerol leads only to increases in τ _R , whereas local motional parameters remain unchanged. Thus the data from multiple samples can be analyzed jointly, with τ _R playing the role of experimentally controlled variable. Based on this concept, the extended model-free model was constructed with the intent to capture the effect of ns time-scale rotameric jumps in valine and leucine side chains. Using this model, we made a positive identification of nanosecond dynamics in Val-23 where ns motions were already observed earlier. In several other cases, however, only tentative identification was possible. The lack of definitive results was due to the approximate character of the model—contrary to what has been assumed, addition of glycerol led to a gradual ‘stiffening’ of the protein. This and other observations also shed light on the interaction of the protein with glycerol, which is one of the naturally occurring osmoprotectants. In particular, it has been found that the overall protein tumbling is controlled by the bulk solvent, and not by a thin solvation layer which contains a higher proportion of water.