A graphical method for analyzing distance restraints using residual dipolar couplings for structure determination of symmetric protein homo-oligomers

High-resolution structure determination of homo-oligomeric protein complexes remains a daunting task for NMR spectroscopists. Although isotope-filtered experiments allow separation of intermolecular NOEs from intramolecular NOEs and determination of the structure of each subunit within the oligomeric state, degenerate chemical shifts of equivalent nuclei from different subunits make it difficult to assign intermolecular NOEs to nuclei from specific pairs of subunits with certainty, hindering structural analysis of the oligomeric state. Here, we introduce a graphical method, DISCO, for the analysis of intermolecular distance restraints and structure determination of symmetric homo-oligomers using residual dipolar couplings. Based on knowledge that the symmetry axis of an oligomeric complex must be parallel to an eigenvector of the alignment tensor of residual dipolar couplings, we can represent distance restraints as annuli in a plane encoding the parameters of the symmetry axis. Oligomeric protein structures with the best restraint satisfaction correspond to regions of this plane with the greatest number of overlapping annuli. This graphical analysis yields a technique to characterize the complete set of oligomeric structures satisfying the distance restraints and to quantitatively evaluate the contribution of each distance restraint. We demonstrate our method for the trimeric E. coli diacylglycerol kinase, addressing the challenges in obtaining subunit assignments for distance restraints. We also demonstrate our method on a dimeric mutant of the immunoglobulin-binding domain B1 of streptococcal protein G to show the resilience of our method to ambiguous atom assignments. In both studies, DISCO computed oligomer structures with high accuracy despite using ambiguously assigned distance restraints.     

Partial alignment and measurement of residual dipolar couplings of proteins under high hydrostatic pressure

High-pressure NMR spectroscopy has emerged as a complementary approach for investigating various structural and thermodynamic properties of macromolecules. Noticeably absent from the array of experimental restraints that have been employed to characterize protein structures at high hydrostatic pressure is the residual dipolar coupling, which requires the partial alignment of the macromolecule of interest. Here we examine five alignment media that are commonly used at ambient pressure for this purpose. We find that the spontaneous alignment of Pf1 phage, d(GpG) and a C12E5/n-hexnanol mixture in a magnetic field is preserved under high hydrostatic pressure. However, DMPC/DHPC bicelles and collagen gel are found to be unsuitable. Evidence is presented to demonstrate that pressure-induced structural changes can be identified using the residual dipolar coupling.     

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.     

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.     

Major groove width variations in RNA structures determined by NMR and impact of <Superscript>13</Superscript>C residual chemical shift anisotropy and <Superscript>1</Superscript>H–<Superscript>13</Superscript>C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by ¹H– ¹H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic ¹H– ¹³C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.     

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.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

NMR solution structure of calerythrin, an EF-hand calcium-binding protein from Saccharopolyspora erythraea.

The structure of calerythrin, a prokaryotic 20 kDa calcium-binding protein has been determined by solution NMR spectroscopy. Distance, dihedral angle, J coupling, secondary chemical shift, residual dipolar coupling and radius of gyration restraints reveal four EF-hand motifs arranged in a compact globular structure. A tight turn in the middle of the amino acid sequence brings the two halves, each comprising a pair of EF-hands, close together. The structural similarity between calerythrin and the eukaryotic sarcoplasmic calcium-binding proteins is notable.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

The catalytic domain of MMP-1 studied through tagged lanthanides

Pseudocontact shifts (pcs) and paramagnetic residual dipolar couplings (rdc) provide structural information that can be used to assess the adequacy of a crystallographic structure to represent the solution structure of a protein. This can be done by attaching a lanthanide binding tag to the protein. There are cases in which only local rearrangements are sufficient to match the NMR data and cases where significant secondary structure or domain rearrangements from the solid state to the solution state are needed. We show that the two cases are easily distinguishable. Whereas the use of solution restraints in the latter case is described in the literature, here we deal with how to obtain a better model of the solution structure in a case (the catalytic domain of the matrix metalloproteinase MMP-1) of the former class.     

Use of residual dipolar couplings as restraints in ab initio protein structure prediction

NMR residual dipolar couplings (RDCs), in the form of the projection angles between the respective internuclear bond vectors, are used as structural restraints in the ab initio structure prediction of a test set of six proteins. The restraints are applied using a recently developed SICHO (SIde-CHain-Only) lattice protein model that employs a replica exchange Monte Carlo (MC) algorithm to search conformational space. Using a small number of RDC restraints, the quality of the predicted structures is improved as reflected by lower RMSD/dRMSD (root mean square deviation/distance root mean square deviation) values from the corresponding native structures and by the higher correlation of the most cooperative mode of motion of each predicted structure with that of the native structure. The latter, in particular, has possible implications for the structure-based functional analysis of predicted structures.     

High-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a hydroxamic acid inhibitor

The high-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a sulfonamide derivative of a hydroxamic acid compound (WAY-151693) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated for residues 7-164 by means of hybrid distance geometry-simulated annealing using a total of 3280 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures is 0.43(+/-0.05) A for the backbone atoms, 0.80(+/-0.09) A for all atoms, and 0.47(+/-0.04) A for all atoms excluding disordered side-chains. The overall structure of MMP-13 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices where its overall fold is consistent with previously solved MMP structures. A comparison of the NMR structure of MMP-13 with the published 1.6 A resolution X-ray structure indicates that the major differences between the structures is associated with loop dynamics and crystal-packing interactions. The side-chains of some active-site residues for the NMR and X-ray structures of MMP-13 adopt distinct conformations. This is attributed to the presence of unique inhibitors in the two structures that encounter distinct interactions with MMP-13. The major structural difference observed between the MMP-13 and MMP-1 NMR structures is the relative size and shape of the S1' pocket where this pocket is significantly longer for MMP-13, nearly reaching the surface of the protein. Additionally, MMP-1 and MMP-13 exhibit different dynamic properties for the active-site loop and the structural Zn-binding region. The inhibitor WAY-151693 is well defined in the MMP-13 active-site based on a total of 52 distance restraints. The binding motif of WAY-151693 in the MMP-13 complex is consistent with our previously reported MMP-1:CGS-27023A NMR structure and is similar to the MMP-13: RS-130830 X-ray structure.     

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.     

Incorporating residual dipolar couplings into the NMR solution structure determination of nucleic acids

NMR solution structures of nucleic acids are generally less well defined than similar-sized proteins. Most NMR structures of nucleic acids are defined only by short-range interactions, such as intrabase-pair or sequential nuclear Overhauser effects (NOEs), and J-coupling constants, and there are no long-range structural data on the tertiary structure. Residual dipolar couplings represent an extremely valuable source of distance and angle information for macromolecules but they average to zero in isotropic solutions. With the recent advent of general methods for partial alignment of macromolecules in solution, residual dipolar couplings are rapidly becoming indispensable constraints for solution NMR structural studies. These residual dipolar couplings give long-range global structural information and thus complement the strictly local structural data obtained from standard NOE and torsion angle constraints. Such global structural data are especially important in nucleic acids due to the more elongated, less-globular structure of many DNAs and RNAs. Here we review recent progress in application of residual dipolar couplings to structural studies of nucleic acids. We also present results showing how refinement procedures affect the final solution structures of nucleic acids.Copyright 2001 John Wiley & Sons, Inc.     

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 determination of the two DNA-binding domains in the Schizosaccharomyces pombe Abp1 protein by a combination of dipolar coupling and diffusion anisotropy restraints

We have solved the solution structure of the N-terminal region of the fission yeast centromere protein, Abp1, bound to a 21-base pair DNA fragment bearing its recognition site (Mw = 30 kDa). Although the two DNA-binding domains in the Abp1 protein were defined well by a conventional NOE-based NMR methodology, the overall structure of the Abp1 protein was poorly defined, due to the lack of interdomain distance restraints. Therefore, we additionally used residual dipolar couplings measured in a weakly aligned state, and rotational diffusion anisotropies. Neither the NH residual dipolar couplings nor the backbone ¹⁵N T ₁/T ₂ data were sufficient to determine the overall structure of the Abp1 protein, due to spectral overlap. We used a combination of these two orientational restraints (residual dipolar coupling and rotational diffusion anisotropy), which significantly improved the convergence of the overall structures. The range of the observed T ₁/T ₂ ratios was wider (20–50 for the secondary structure regions of Abp1) than the previously reported data for several globular proteins, indicating that the overall shape of the Abp1DNA complex is ellipsoid. This extended form would facilitate the recognition of the two separate sites in the relatively long DNA sequence by the DNA-binding domains of Apb1.     

Automated structure determination from NMR spectra

Automated methods for protein structure determination by NMR have increasingly gained acceptance and are now widely used for the automated assignment of distance restraints and the calculation of three-dimensional structures. This review gives an overview of the techniques for automated protein structure analysis by NMR, including both NOE-based approaches and methods relying on other experimental data such as residual dipolar couplings and chemical shifts, and presents the FLYA algorithm for the fully automated NMR structure determination of proteins that is suitable to substitute all manual spectra analysis and thus overcomes a major efficiency limitation of the NMR method for protein structure determination.