NMR investigations of allosteric processes in a two-domain Thermus thermophilus Hsp70 molecular chaperone

Hsp70 chaperones are two-domain proteins that assist in intra-cellular protein (re) folding processes in all species. The protein folding activity of the substrate binding domain of the Hsp70s is regulated by nucleotide binding at the nucleotide-binding domain through an as yet undefined heterotropic allosteric mechanism. The available structures of the isolated domains of Hsp70s have given very limited indications of nucleotide-induced conformational changes that could modulate the affinity for substrate proteins. Here, we present a multi-dimensional NMR study of a prokaryotic Hsp70 homolog, Thermus thermophilus DnaK, using a 54kDa construct containing both nucleotide binding domain and most of the substrate binding domain. It is determined that the nucleotide binding domain and substrate binding domain are closely associated in all ligand states studied. Comparison of the assigned NMR spectra of the two-domain construct with those of the previously studied isolated nucleotide binding domain, allowed the identification of the nucleotide binding domain-substrate binding domain interface. A global three-dimensional structure was obtained for the two-domain construct on the basis of this information and of NMR residual dipolar couplings measurements. This is the first experimental elucidation of the relative positioning of the nucleotide binding domain and substrate binding domain for any Hsp70 chaperone. Comparisons of NMR data between various ligand states including nucleotide-free, ATP, ADP.Pi and ADP.Pi+ peptide bound, identified residues involved in the allosteric inter-domain communication. In particular, peptide binding to the substrate binding domain was found to cause conformational changes in the NBD extending to the nucleotide binding pocket. Detailed analysis suggests that the inter-domain interface becomes tighter in the (nucleotide binding domain ligation/substrate binding domain ligation) order ATP/apo, ADP.Pi/apo ADP.Pi/peptide.     

Various strategies of using residual dipolar couplings in NMR-driven protein docking: application to Lys48-linked di-ubiquitin and validation against 15N-relaxation data

When classical, Nuclear Overhauser Effect (NOE)-based approaches fail, it is possible, given high-resolution structures of the free molecules, to model the structure of a complex in solution based solely on chemical shift perturbation (CSP) data in combination with orientational restraints from residual dipolar couplings (RDCs) when available. RDCs can be incorporated into the docking following various strategies: as direct restraints and/or as intermolecular intervector projection angle restraints (Meiler et al., J Biomol NMR 2000;16:245-252). The advantage of the latter for docking is that they directly define the relative orientation of the molecules. A combined protocol in which RDCs are first introduced as intervector projection angle restraints and at a later stage as direct restraints is shown here to give the best performance. This approach, implemented in our information-driven docking approach HADDOCK (Dominguez et al., J Am Chem Soc 2003;125:1731-1737), is used to determine the solution structure of the Lys48-linked di-ubiquitin, for which chemical shift mapping, RDCs, and (15)N-relaxation data have been previously obtained (Varadan et al., J Mol Biol 2002;324:637-647). The resulting structures, derived from CSP and RDC data, are cross-validated using (15)N-relaxation data. The solution structure differs from the crystal structure by a 20 degrees rotation of the two ubiquitin units relative to each other.     

SOFAST-HMQC Experiments for Recording Two-dimensional Deteronuclear Correlation Spectra of Proteins within a Few Seconds

Fast multidimensional NMR with a time resolution of a few seconds provides a new tool for high throughput screening and site-resolved real-time studies of kinetic molecular processes by NMR. Recently we have demonstrated the feasibility to record protein ¹H–¹⁵N correlation spectra in a few seconds of acquisition time using a new SOFAST-HMQC experiment (Schanda and Brutscher (2005) J. Am. Chem. Soc. 127, 8014). Here, we investigate in detail the performance of SOFAST-HMQC to record ¹H–¹⁵N and ¹H−¹³C correlation spectra of proteins of different size and at different magnetic field strengths. Compared to standard ¹H–¹⁵N correlation experiments SOFAST-HMQC provides a significant gain in sensitivity, especially for fast repetition rates. Guidelines are provided on how to set up SOFAST-HMQC experiments for a given protein sample. In addition, an alternative pulse scheme, IPAP-SOFAST-HMQC is presented that allows application on NMR spectrometers equipped with cryogenic probes, and fast measurement of one-bond ¹H–¹³C and ¹H–¹⁵N scalar and residual dipolar coupling constants.     

Determination of helical membrane protein topology using residual dipolar couplings and exhaustive search algorithm: application to phospholamban

Dipolar waves are distinct hallmarks of both the secondary and tertiary structures of alpha-helical proteins that are immobilized in membrane bilayers or embedded in anisotropic media. We present a simple, semi-empirical approach that exploits the modulation of the amplitude and average of dipolar waves to determine the topology of alpha-helical proteins. Moreover, we describe the application of this method for the structural determination of a detergent solubilized membrane protein, phospholamban (PLB) that is involved in calcium regulation of cardiac muscle. When combined with high-resolution solid-state NMR data, this method can serve as a fast route for determining the topology of helical membrane proteins solubilized in detergent micelles.     

Compensating bends in a 16-base-pair DNA oligomer containing a T(3)A(3) segment: A NMR study of global DNA curvature

In-phase ligated DNA containing T(n)A(n) segments fail to exhibit the retarded polyacrylamide gel electrophoresis (PAGE) migration observed for in-phase ligated A(n)T(n) segments, a behavior thought to be correlated with macroscopic DNA curvature. The lack of macroscopic curvature in ligated T(n)A(n) segments is thought to be due to cancellation of bending in regions flanking the TpA steps. To address this issue, solution-state NMR, including residual dipolar coupling (RDC) restraints, was used to determine a high-resolution structure of [d(CGAGGTTTAAACCTCG)2], a DNA oligomer containing a T3A3 tract. The overall magnitude and direction of bending, including the regions flanking the central TpA step, was measured using a radius of curvature, Rc, analysis. The Rc for the overall molecule indicated a small magnitude of global bending (Rc = 138 +/- 23 nm) towards the major groove, whereas the Rc for the two halves (72 +/- 33 nm and 69 +/- 14 nm) indicated greater localized bending into the minor groove. The direction of bending in the regions flanking the TpA step is in partial opposition (109 degrees), contributing to cancellation of bending. The cancellation of bending did not correlate with a pattern of roll values at the TpA step, or at the 5' and 3' junctions, of the T3A3 segment, suggesting a simple junction/roll model is insufficient to predict cancellation of DNA bending in all T(n)A(n) junction sequence contexts. Importantly, Rc analysis of structures refined without RDC restraints lacked the precision and accuracy needed to reliably measure bending.     

A target function for quaternary structural refinement from small angle scattering and NMR orientational restraints

We present a novel target function based on atomic coordinates that permits quaternary structural refinement of multi-domain protein–protein or protein–RNA complexes. It requires that the high-resolution structures of the individual domains are known and that small angle scattering (SAS) data as well as NMR orientational restraints from residual dipolar couplings (RDCs) of the complex are available. We show that, when used in combination, the translational and rotational restraints contained in SAS intensities and RDCs, respectively, define a target potential function that permits to determine the overall topology of complexes made up of domains with low internal symmetry. We apply the target function on a modestly anisotropic model system, the Barnase/Barstar complex, and discuss factors that influence the structural refinement such as data errors and the geometrical properties of the individual domains.     

Influence of the fluctuations of the alignment tensor on the analysis of the structure and dynamics of proteins using residual dipolar couplings

It has been suggested that the fluctuations of the alignment tensor can affect the results of procedures for characterizing the structure and the dynamics of proteins using residual dipolar couplings. We show here that the very significant fluctuations of the steric alignment tensor caused by the dynamics of proteins can be safely ignored when they do not correlate with those of the bond vectors. A detailed analysis of these correlations in the protein ubiquitin reveals that their effects are negligible for the analysis of backbone motions within secondary structure elements, but also that they may be significant in turns, loops and side chains, especially for bond vectors that have small residual dipolar couplings. Our results suggest that methods that explicitly consider the motions of the alignment tensor will be needed to study the large-scale structural fluctuations that take place on the millisecond timescale, which are often important for the biological function of proteins, from residual dipolar coupling measurements.     

A tool for the prediction of structures of complex sugars

In two recent back to back articles(Xia et al., J Chem Theory Comput 3:1620–1628 and 1629–1643, 2007a, b) we have started to address the problem of complex oligosaccharide conformation and folding. The scheme previously presented was based on exhaustive searches in configuration space in conjunction with Nuclear Overhauser Effect (NOE) calculations and the use of a complex rotameric library that takes branching into account. NOEs are extremely useful for structural determination but only provide information about short range interactions and ordering. Instead, the measurement of residual dipolar couplings (RDC), yields information about molecular ordering or folding that is long range in nature. In this article we show the results obtained by incorporation RDC calculations into our prediction scheme. Using this new approach we are able to accurately predict the structure of six human milk sugars: LNF-1, LND-1, LNF-2, LNF-3, LNnT and LNT. Our exhaustive search in dihedral configuration space combined with RDC and NOE calculations allows for highly accurate structural predictions that, because of the non-ergodic nature of these molecules on a time scale compatible with molecular dynamics simulations, are extremely hard to obtain otherwise (Almond et al., Biochemistry 43:5853–5863, 2004). Molecular dynamics simulations in explicit solvent using as initial configurations the structures predicted by our algorithm show that the histo-blood group epitopes in these sugars are relatively rigid and that the whole family of oligosaccharides derives its conformational variability almost exclusively from their common linkage (β-d-GlcNAc-(1→3)-β-d-Gal) which can exist in two distinct conformational states. A population analysis based on the conformational variability of this flexible glycosidic link indicates that the relative population of the two distinct states varies for different human milk oligosaccharides. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.