Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

What is the average conformation of bacteriophage T4 lysozyme in solution? A domain orientation study using dipolar couplings measured by solution NMR

Lysozyme from T4 bacteriophage is comprised of two domains that are both involved in binding substrate. Although wild-type lysozyme has been exclusively crystallized in a closed form that is similar to the peptidoglycan-bound conformation, a more open structure is thought to be required for ligand binding. To determine the relative arrangement of domains within T4 lysozyme in the solution state, dipolar couplings were measured in several different dilute liquid crystalline media by solution NMR methods. The dipolar coupling data were analyzed with a domain orientation procedure described previously that utilizes high- resolution X-ray structures. The cleft between the domains is significantly larger in the average solution structure than what is observed in the X-ray structure of the ligand-free form of the protein (approximately 17 degrees closure from solution to X-ray structures). A comparison of the solution domain orientation with X-ray-derived structures in the protein data base shows that the solution structure resembles a crystal structure obtained for the M6I mutant. Dipolar couplings were also measured on the lysozyme mutant T21C/T142C, which was oxidized to form an inter-domain disulfide bond (T4SS). In this case, the inter-domain solution structure was found to be more closed than was observed in the crystal (approximately 11 degrees). Direct refinement of lysozyme crystal structures with the measured dipolar couplings using the program CNS, establishes that this degree of closure can be accommodated whilst maintaining the inter-domain cystine bond. The differences between the average solution conformations obtained using dipolar couplings and the crystal conformations for both forms of lysozyme investigated in this study illustrate the impact that crystal packing interactions can have on the arrangement of domains within proteins and the importance of alternative methods to X-ray crystallography for evaluating inter-domain structure.     

A physical picture of atomic motions within the Dickerson DNA dodecamer in solution derived from joint ensemble refinement against NMR and large-angle X-ray scattering data

The structure and dynamics of the Dickerson DNA dodecamer [5'd(CGCGAATTCGCG)2] in solution have been investigated by joint simulated annealing refinement against NMR and large-angle X-ray scattering data (extending from 0.25 to 3 A-1). The NMR data comprise an extensive set of hetero- and homonuclear residual dipolar coupling and 31P chemical shift anisotropy restraints in two alignment media, supplemented by NOE and 3J coupling data. The NMR and X-ray scattering data cannot be fully ascribed to a single structure representation, indicating the presence of anisotropic motions that impact the experimental observables in different ways. Refinement with ensemble sizes (Ne) of >or=2 to represent the atomic motions reconciles all the experimental data within measurement error. Cross validation against both the dipolar coupling and X-ray scattering data suggests that the optimal ensemble size required to account for the current data is 4. The resulting ensembles permit one to obtain a detailed view of the conformational space sampled by the dodecamer in solution and permit one to analyze fluctuations in helicoidal parameters, sugar puckers, and BI-BII backbone transitions and to obtain quantitative metrics of atomic motion such as generalized order parameters and thermal B factors. The calculated order parameters are in good agreement with experimental order parameters obtained from 13C relaxation measurements. Although DNA behaves as a relatively rigid rod with a persistence length of approximately 150 bp, dynamic conformational heterogeneity at the base pair level is functionally important since it readily permits optimization of intermolecular protein-DNA interactions.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Global folds of proteins with low densities of NOEs using residual dipolar couplings: application to the 370-residue maltodextrin-binding protein

The global fold of maltose-binding protein in complex with the substrate beta-cyclodextrin was determined by solution NMR methods. The two-domain protein is comprised of a single polypeptide chain of 370 residues, with a molecular mass of 42 kDa. Distance information in the form of H(N)-H(N), H(N)-CH(3) and CH(3)-CH(3) NOEs was recorded on (15)N, (2)H and (15)N, (13)C, (2)H-labeled proteins with methyl protonation in Val, Leu, and Ile (C(delta1) only) residues. Distances to methyl protons, critical for the structure determination, comprised 77 % of the long-range restraints. Initial structures were calculated on the basis of 1943 NOEs, 48 hydrogen bond and 555 dihedral angle restraints. A global pair-wise backbone rmsd of 5.5 A was obtained for these initial structures with rmsd values for the N and C domains of 2.4 and 3.8 A, respectively. Direct refinement against one-bond (1)H(N)-(15)N, (13)C(alpha)-(13)CO, (15)N-(13)CO, two-bond (1)H(N)-(13)CO and three-bond (1)H(N)-(13)C(alpha) dipolar couplings resulted in structures with large numbers of dipolar restraint violations. As an alternative to direct refinement against measured dipolar couplings we have developed an approach where discrete orientations are calculated for each peptide plane on the basis of the dipolar couplings described above. The orientation which best matches that in initial NMR structures calculated from NOE and dihedral angle restraints exclusively is used to refine further the structures using a new module written for CNS. Modeling studies from four different proteins with diverse structural motifs establishes the utility of the methodology. When applied to experimental data recorded on MBP the precision of the family of structures generated improves from 5.5 to 2.2 A, while the rmsd with respect to the X-ray structure (1dmb) is reduced from 5.1 to 3.3 A.     

Measurement of 15N relaxation in deuterated amide groups in proteins using direct nitrogen detection

15N chemical shielding tensors contain useful structural information, and their knowledge is essential for accurate analysis of protein backbone dynamics. The anisotropic component (CSA) of 15N chemical shielding can be obtained from 15N relaxation measurements in solution. However, the predominant contribution to nitrogen relaxation from 15N-(1)H dipolar coupling in amide groups limits the sensitivity of these measurements to the actual CSA values. Here we present nitrogen-detected NMR experiments for measuring 15N relaxation in deuterated amide groups in proteins, where the dipolar contribution to 15N relaxation is significantly reduced by the deuteration. Under these conditions nitrogen spin relaxation becomes a sensitive probe for variations in 15N chemical shielding tensors. Using the nitrogen direct-detection experiments we measured the rates of longitudinal and transverse 15N relaxation for backbone amides in protein G in D(2)O at 11.7 T. The measured relaxation rates are validated by comparing the overall rotational diffusion tensor obtained from these data with that from the conventional 15N relaxation measurements in H(2)O. This analysis revealed a 17-24 degree angle between the NH-bond and the unique axis of the 15N chemical shielding tensor.     

Force field dependence of NMR-Based,restrained molecular dynamics DNA structure calculations including an analysis of the influence of residual dipolar coupling restraints

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.     

Structure and dynamics of micelle-associated human immunodeficiency virus gp41 fusion domain

The N-terminal fusion domain of the HIV-1 gp41 envelope glycoprotein is responsible for initiating the fusion of viral and cellular membranes, leading to the subsequent infection of the host cell by HIV-1. We have investigated the backbone structure and dynamics of the 30 N-terminal residues of HIV-1 gp41 in membrane-mimicking environments using NMR spectroscopy and (15)N- and (15)N,(13)C,(2)H-labeled peptides. Similar (15)N-(1)H HSQC spectra were obtained in a variety of detergents, including SDS, DPC, mixed DPC/SDS, and LPPG micelles, indicating that the peptide structure is not strongly influenced by the type of detergent used. Detailed characterization was carried out in SDS micelles, where the long-term sample stability was found to be optimal. In addition to J-coupling and NOE restraints, a nearly complete set of backbone residual dipolar coupling restraints was recorded for the fusion domain-micelle complex aligned with respect to the magnetic field using a stretched polyacrylamide gel. Backbone amide (15)N spin relaxation and amide hydrogen exchange rates with the solvent were also measured. The ensemble of NMR structures reveals an uninterrupted alpha-helix for the least mobile residues (S(2) > 0.65), Ile-4 to Met-19, with transient helical character extending up to Ala-22. A 12-residue (Ile-4 to Ala-15) segment is fully shielded from solvent, with Gly-3 and Gly-16 found at micelle-solvent interfaces. Residues external to the micelle exhibit enhanced picosecond to nanosecond time scale dynamics relative to the residues buried in the micelle, and their mobility increases with the distance from the micelle.     

Force Field Dependence of NMR-Based,Restrained Molecular Dynamics DNA Structure Calculations Including an Analysis of the Influence of Residual Dipolar Coupling Restraints

Abstract

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.     

Solution structure and backbone dynamics of the TGFbeta type II receptor extracellular domain

Isoforms of transforming growth factor beta (TGFbeta) are 25 kDa homodimeric polypeptides that signal by binding and bringing together two related, functionally distinct cell surface receptors designated as TbetaR1 and TbetaR2. Here, we report the solution structure of the 13.8 kDa extracellular domain of human TbetaR2 (ecTbetaR2) as calculated from N(N)-H(N), C(alpha)-H(alpha), and C(alpha)-C(O) residual dipolar coupling restraints in conjunction with NOE distance, dihedral angle, and scalar coupling restraints. Comparison of the free ecTbetaR2 solution structure with the TGFbeta3-bound ecTbetaR2 crystal structure reveals backbone conformations that superimpose with RMSDs of 1.0 A over the regions of regular secondary structure and 1.4 A overall. The differences in structure fall mainly in loop regions that are either poorly defined by the available NMR data or are involved in crystal contacts. The noted similarities between the NMR structure of the free form and the crystal structure of the TGFbeta-bound form are also consistent with the close correspondence, 0.16 A RMSD for regions of secondary structure and 0.51 A RMSD overall, for the crystal structure of free ecTbetaR2 as compared to the crystal structure of TGFbeta3-bound ecTbetaR2. Despite the apparent similarities between the free and the bound forms, there appears to be small but significant differences in structure involving the interfacial contact region of the receptor. Measurements of backbone (15)N relaxation times and interpretation of these by the model-free formalism with axial diffusional anisotropy further reveal significant ms to micros time scale motions centered about two of the conserved disulfide bonds and in several residues that comprise the TGFbeta binding surface. Together, these observations indicate that binding likely occurs through a mechanism with a small component of induced fit character, whereby flexibility within the receptor facilitates the transition to the TGFbeta-bound state.     

The structure and dynamics of rat apo-cellular retinol-binding protein II in solution: comparison with the X-ray structure

The structure and dynamics of rat apo-cellular retinol binding protein II (apo-CRBP II) in solution has been determined by multidimensional NMR analysis of uniformly enriched recombinant rat 13C, 15N-apo-CRBP II and 15N-apo-CRBP II. The final ensemble of 24 NMR structures has been calculated from 3274 conformational restraints or 24.4 restraints/residue. The average root-mean-square deviation of the backbone atoms for the final 24 structures relative to their mean structure is 1.06 A. Although the average solution structure is very similar to the crystal structure, it differs at the putative entrance to the binding cavity, which is formed by the helix-turn-helix motif, the betaC-betaD turn and the betaE-betaF turn. The mean coordinates of the main-chain atoms of amino acid residues 28-38 are displaced in the solution structure relative to the crystal structure. The side-chain of F58, located on the betaC-betaD turn, is reoriented such that it interacts with L37 and no longer blocks entry into the ligand-binding pocket. Residues 28-35, which form the second helix of the helix-turn-helix motif in the crystal structure, do not exhibit a helical conformation in the solution structure. The solution structure of apo-CRBP II exhibits discrete regions of backbone disorder which are most pronounced at residues 28-32, 37-38 and 73-76 in the betaE-betaF turn as evaluated by the consensus chemical shift index, the root-mean-square deviation, amide 1H exchange rates and 15N relaxation studies. These studies indicate that fluctuations in protein conformation occur on the microseconds to ms time-scale in these regions of the protein. Some of these exchange processes can be directly observed in the three-dimensional 15N-resolved NOESY spectrum. These results suggest that in solution, apo-CRBP II undergoes conformational changes on the microseconds to ms time-scale which result in increased access to the binding cavity.     

A structure refinement protocol combining NMR residual dipolar couplings and small angle scattering restraints

We present the implementation of a target function based on Small Angle Scattering data (Gabel et al. Eur Biophys J 35(4):313-327, 2006) into the Crystallography and NMR Systems (CNS) and demonstrate its utility in NMR structure calculations by simultaneous application of small angle scattering (SAS) and residual dipolar coupling (RDC) restraints. The efficiency and stability of the approach are demonstrated by reconstructing the structure of a two domain region of the 31 kDa nuclear export factor TAP (TIP-associated protein). Starting with the high resolution X-ray structures of the two individual TAP domains, the translational and orientational domain arrangement is refined simultaneously. We tested the stability of the protocol against variations of the SAS target parameters and the number of RDCs and their uncertainties. The activation of SAS restraints results in an improved translational clustering of the domain positions and lifts part of the fourfold degeneracy of their orientations (associated with a single alignment tensor). The resulting ensemble of structures reflects the conformational space that is consistent with the experimental SAS and RDC data. The SAS target function is computationally very efficient. SAS restraints can be activated at different levels of precision and only a limited SAS angular range is required. When combined with additional data from chemical shift perturbation, paramagnetic relaxation enhancement or mutational analysis the SAS refinement is an efficient approach for defining the topology of multi-domain and/or multimeric biomolecular complexes in solution based on available high resolution structures (NMR or X-ray) of the individual domains.     

Solution structure of a phage-derived peptide antagonist in complex with vascular endothelial growth factor

Vascular endothelial growth factor (VEGF) is a potent endothelial cell-specific mediator of angiogenesis and vasculogenesis. VEGF is involved pathologically in cancer, proliferative retinopathy and rheumatoid arthritis, and as such represents an important therapeutic target. Three classes of disulfide-constrained peptides that antagonize binding of the VEGF dimer to its receptors, KDR and Flt-1, were identified previously using phage display methods. NMR studies of a representative peptide from the most potent class of these peptide antagonists, v107 (GGNECDAIRMWEWECFERL), were undertaken to characterize its interactions with VEGF. v107 has no defined structure free in solution, but binding to VEGF induces folding of the peptide. The solution structure of the VEGF receptor-binding domain-v107 complex was determined using 3940 (1970 per VEGF monomer) internuclear distance and 476 (238 per VEGF monomer) dihedral angle restraints derived from NMR data obtained using samples containing either (13)C/(15)N-labeled protein plus excess unlabeled peptide or (13)C/(15)N-labeled peptide plus excess unlabeled protein. Residual dipolar coupling restraints supplemented the structure determination of the complex and were found to increase significantly both the global precision of VEGF in the complex and the agreement with available crystal structures of VEGF. The calculated ensemble of structures is of high precision and is in excellent agreement with the experimental restraints. v107 has a turn-helix conformation with hydrophobic residues partitioned to one face of the peptide and polar or charged residues at the other face. Contacts between two v107 peptides and the VEGF dimer are mediated by primarily hydrophobic side-chain interactions. The v107-binding site on VEGF overlaps partially with the binding site of KDR and is similar to that for domain 2 of Flt-1. The structure of the VEGF-v107 complex provides new insight into how binding to VEGF can be achieved that may be useful for the design of small molecule antagonists.     

1H NMR studies of human lysozyme: spectral assignment and comparison with hen lysozyme

Complete main-chain (NH and alpha CH) 1H NMR assignments are reported for the 130 residues of human lysozyme, along with extensive assignments for side-chain protons. Analysis of 2-D NOESY experiments shows that the regions of secondary structure for human lysozyme in solution are essentially identical with those found previously in a similar study of hen lysozyme and are in close accord with the structure of the protein reported previously from X-ray diffraction studies in the crystalline state. Comparison of the chemical shifts, spin-spin coupling constants, and hydrogen exchange behavior are also consistent with closely similar structures for the two proteins in solution. In a number of cases specific differences in the NMR parameters between hen and human lysozymes can be correlated with specific differences observed in the crystal structures.     

Structure of hen egg-white lysozyme solvated in TFE/water: a molecular dynamics simulation study based on NMR data

Various experimental studies of hen egg white lysozyme (HEWL) in water and TFE/water clearly indicate structural differences between the native state and TFE state of HEWL, e.g. the helical content of the protein in the TFE state is much higher than in the native state. However, the available detailed NMR studies were not sufficient to determine fully a structure of HEWL in the TFE state. Different molecular dynamics (MD) simulations, i.e. at room temperature, at increased temperature and using proton–proton distance restraints derived from NMR NOE data, have been used to generate configurational ensembles corresponding to the TFE state of HEWL. The configurational ensemble obtained at room temperature using atom-atom distance restraints measured for HEWL in TFE/water solution satisfies the experimental data and has the lowest protein energy. In this ensemble residues 50–58, which are part of the β-sheet in native HEWL, adopt fluctuating α-helical secondary structure.     

Refinement of the solution structure of the heparin-binding domain of vascular endothelial growth factor using residual dipolar couplings

Previous NMR structural studies of the heparin-binding domain of vascular endothelial growth factor (VEGF₁₆₅) revealed a novel fold comprising two subdomains, each containing two disulfide bridges and a short two-stranded antiparallel -sheet. The mutual orientation of the two subdomains was poorly defined by the NMR data. Heteronuclear relaxation data suggested that this disorder resulted from a relative lack of experimental restraints due to the limited size of the interface, rather than inherent high-frequency flexibility. Refinement of the structure using ¹H^N-¹⁵N residual dipolar coupling restraints results in significantly improved definition of the relative subdomain orientations.     

Refinement of protein structure against non-redundant carbonyl <Superscript>13</Superscript>C NMR relaxation

Carbonyl ¹³C′ relaxation is dominated by the contribution from the ¹³C′ chemical shift anisotropy (CSA). The relaxation rates provide useful and non-redundant structural information in addition to dynamic parameters. It is straightforward to acquire, and offers complimentary structural information to the ¹⁵N relaxation data. Furthermore, the non-axial nature of the ¹³C′ CSA tensor results in a T₁/T₂ value that depends on an additional angular variable even when the diffusion tensor of the protein molecule is axially symmetric. This dependence on an extra degree of freedom provides new geometrical information that is not available from the NH dipolar relaxation. A protocol that incorporates such structural restraints into NMR structure calculation was developed within the program Xplor-NIH. Its application was illustrated with the yeast Fis1 NMR structure. Refinement against the ¹³C′ T₁/T₂ improved the overall quality of the structure, as evaluated by cross-validation against the residual dipolar coupling as well as the ¹⁵N relaxation data. In addition, possible variations of the CSA tensor were addressed. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Response of native and denatured hen lysozyme to high pressure studied by (15)N/(1)H NMR spectroscopy.

High-pressure (15)N/(1)H NMR techniques were used to characterize the conformational fluctuations of hen lysozyme, in its native state and when denatured in 8 M urea, over the pressure range 30--2000 bar. Most (1)H and (15)N signals of native lysozyme show reversible shifts to low field with increasing pressure, the average pressure shifts being 0.069 +/- 0.101 p.p.m. ((1)H) and 0.51 +/- 0.36 p.p.m. ((15)N). The shifts indicate that the hydrogen bonds formed to carbonyl groups or water molecules by the backbone amides are, on average, shortened by approximately 0.02 A as a result of pressure. In native lysozyme, six residues in the beta domain or at the alpha/beta domain interface have anomalously large nonlinear (15)N and (1)H chemical-shift changes. All these residues lie close to water-containing cavities, suggesting that there are conformational changes involving these cavities, or the water molecules within them, at high pressure. The pressure-induced (1)H and (15)N shifts for lysozyme denatured in 8 M urea are much more uniform than those for native lysozyme, with average backbone amide shifts of 0.081 +/- 0.029 p.p.m. ((1)H) and 0.57 +/- 0.14 p.p.m. ((15)N). The results show that overall there are no significant variations in the local conformational properties of denatured lysozyme with pressure, although larger shifts in the vicinity of a persistent hydrophobic cluster indicate that interactions in this part of the sequence may rearrange. NMR diffusion measurements demonstrate that the effective hydrodynamic radius of denatured lysozyme, and hence the global properties of the denatured ensemble, do not change detectably at high pressure.     

Determination of dihedral Psi angles in large proteins by combining NH(N)/C(alpha)H(alpha) dipole/dipole cross-correlation and chemical shifts

We propose a strategy based on the combination of experimental NH(N)/C(alpha)H(alpha) dipole/dipole cross-correlated relaxation rates and chemical shift analysis for the determination of Psi torsion angles in proteins. The method allows the determination of a dihedral angle that is not easily accessible by nuclear magnetic resonance (NMR). The measurement of dihedral angle restraints can be used for structure calculation, which is known to improve the quality of NMR structures. The method is of particular interest in the case of large proteins, for which spectral assignment of the nuclear Overhauser effect spectra, and therefore straightforward structural determination, is out of reach. One advantage of the method is that it is reasonably simple to implement, and could be used in association with other methods aiming at obtaining structural information on complex systems, such as residual dipolar coupling measurements. An illustrative example is analyzed in the case of the 30-kDa protein 6-phosphogluconolactonase.     

Quaternary structure built from subunits combining NMR and small-angle x-ray scattering data

A new principle in constructing molecular complexes from the known high-resolution domain structures joining data from NMR and small-angle x-ray scattering (SAXS) measurements is described. Structure of calmodulin in complex with trifluoperazine was built from N- and C-terminal domains oriented based on residual dipolar couplings measured by NMR in a dilute liquid crystal, and the overall shape of the complex was derived from SAXS data. The residual dipolar coupling data serves to reduce angular degrees of freedom, and the small-angle scattering data serves to confine the translational degrees of freedom. The complex built by this method was found to be consistent with the known crystal structure. The study demonstrates how approximate tertiary structures of modular proteins or quaternary structures composed of subunits can be assembled from high-resolution structures of domains or subunits using mutually complementary NMR and SAXS data.