DNA-nanotube-induced alignment of membrane proteins for NMR structure determination

Membrane proteins are encoded by 20-35% of genes but represent <1% of known protein structures to date. Thus, improved methods for membrane-protein structure determination are of critical importance. Residual dipolar couplings (RDCs), commonly measured for biological macromolecules weakly aligned by liquid-crystalline media, are important global angular restraints for NMR structure determination. For alpha-helical membrane proteins >15 kDa in size, Nuclear-Overhauser effect-derived distance restraints are difficult to obtain, and RDCs could serve as the main reliable source of NMR structural information. In many of these cases, RDCs would enable full structure determination that otherwise would be impossible. However, none of the existing liquid-crystalline media used to align water-soluble proteins are compatible with the detergents required to solubilize membrane proteins. We report the design and construction of a detergent-resistant liquid crystal of 0.8-microm-long DNA-nanotubes that can be used to induce weak alignment of membrane proteins. The nanotubes are heterodimers of 0.4-microm-long six-helix bundles each self-assembled from a 7.3-kb scaffold strand and >170 short oligonucleotide staple strands. We show that the DNA-nanotube liquid crystal enables the accurate measurement of backbone N(H) and C(alpha)H(alpha) RDCs for the detergent-reconstituted zeta-zeta transmembrane domain of the T cell receptor. The measured RDCs validate the high-resolution structure of this transmembrane dimer. We anticipate that this medium will extend the advantages of weak alignment to NMR structure determination of a broad range of detergent-solubilized membrane proteins.     

Statistical coil model of the unfolded state: resolving the reconciliation problem

An unfolded state ensemble is generated by using a self-avoiding statistical coil model that is based on backbone conformational frequencies in a coil library, a subset of the Protein Data Bank. The model reproduces two apparently contradicting behaviors observed in the chemically denatured state for a variety of proteins, random coil scaling of the radius of gyration and the presence of significant amounts of local backbone structure (NMR residual dipolar couplings). The most stretched members of our unfolded ensemble dominate the residual dipolar coupling signal, whereas the uniformity of the sign of the couplings follows from the preponderance of polyproline II and beta conformers in the coil library. Agreement with the NMR data substantially improves when the backbone conformational preferences include correlations arising from the chemical and conformational identity of neighboring residues. Although the unfolded ensembles match the experimental observables, they do not display evidence of native-like topology. By providing an accurate representation of the unfolded state, our statistical coil model can be used to improve thermodynamic and kinetic modeling of protein folding.     

New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary ¹³C shifts, ¹H shifts, T₁/T₂ ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.     

Relation between native ensembles and experimental structures of proteins

Different experimental structures of the same protein or of proteins with high sequence similarity contain many small variations. Here we construct ensembles of "high-sequence similarity Protein Data Bank" (HSP) structures and consider the extent to which such ensembles represent the structural heterogeneity of the native state in solution. We find that different NMR measurements probing structure and dynamics of given proteins in solution, including order parameters, scalar couplings, and residual dipolar couplings, are remarkably well reproduced by their respective high-sequence similarity Protein Data Bank ensembles; moreover, we show that the effects of uncertainties in structure determination are insufficient to explain the results. These results highlight the importance of accounting for native-state protein dynamics in making comparisons with ensemble-averaged experimental data and suggest that even a modest number of structures of a protein determined under different conditions, or with small variations in sequence, capture a representative subset of the true native-state ensemble.     

Quaternary structure of hemoglobin in solution

Many important proteins perform their physiological functions under allosteric control, whereby the binding of a ligand at a specific site influences the binding affinity at a different site. Allosteric regulation usually involves a switch in protein conformation upon ligand binding. The energies of the corresponding structures are comparable, and, therefore, the possibility that a structure determined by x-ray diffraction in the crystalline state is influenced by its intermolecular contacts, and thus differs from the solution structure, cannot be excluded. Here, we demonstrate that the quaternary structure of tetrameric human normal adult carbonmonoxy-hemoglobin can readily be determined in solution at near-physiological conditions of pH, ionic strength, and temperature by NMR measurement of (15)N-(1)H residual dipolar couplings in weakly oriented samples. The structure is found to be a dynamic intermediate between two previously solved crystal structures, known as the R and R2 states. Exchange broadening at the subunit interface points to a rapid equilibrium between different structures that presumably include the crystallographically observed states.     

Dipole tensor-based atomic-resolution structure determination of a nanocrystalline protein by solid-state NMR

Magic-angle spinning (MAS) solid-state NMR (SSNMR) techniques have emerged in recent years for solving complete structures of uniformly labeled proteins lacking macroscopic order. Strategies used thus far have relied primarily on semiquantitative distance restraints, analogous to the nuclear Overhauser effect (NOE) routinely used in solution NMR. Here, we present a complementary approach for using relative orientations of molecular fragments, determined from dipolar line shapes. Whereas SSNMR distance restraints typically have an uncertainty of approximately 1 A, the tensor-based experiments report on relative vector (pseudobond) angles with precision of a few degrees. By using 3D techniques of this type, vector angle (VEAN) restraints were determined for the majority of the 56-residue B1 immunoglobulin binding domain of protein G [protein GB1 (a total of 47 HN-HN, 49 HN-HC, and 12 HA-HB restraints)]. By using distance restraints alone in the structure calculations, the overall backbone root-mean-square deviation (bbRMSD) was 1.01 +/- 0.13 A (1.52 +/- 0.12 A for all heavy atoms), which improved to 0.49 +/- 0.05 A (1.19 +/- 0.07 A) on the addition of empirical chemical shift [torsion angle likelihood obtained from shift and sequence similarity (TALOS)] restraints. VEAN restraints further improved the ensemble to 0.31 +/- 0.06 A bbRMSD (1.06 +/- 0.07 A); relative to the structure with distances alone, most of the improvement remained (bbRMSD 0.64 +/- 0.09 A; 1.29 +/- 0.07 A) when TALOS restraints were removed before refinement. These results represent significant progress toward atomic-resolution protein structure determination by SSNMR, capabilities that can be applied to a large range of membrane proteins and fibrils, which are often not amenable to solution NMR or x-ray crystallography.     

NMR structural and kinetic characterization of a homeodomain diffusing and hopping on nonspecific DNA

Nonspecific protein-DNA interactions are inherently dynamic and involve both diffusion of the protein along the DNA and hopping of the protein from one DNA molecule or segment to another. Understanding how gene regulatory proteins interact nonspecifically with DNA in terms of both structure and dynamics is challenging because the experimental observables are an ensemble average of many rapidly exchanging states. By using a variety of NMR spectroscopic techniques, including relaxation analysis, paramagnetic relaxation enhancement, and residual dipolar couplings, we have characterized structural and kinetic aspects of the interaction of the HoxD9 homeodomain with a nonspecific, 24-bp DNA duplex in a system in which the protein is not constrained to any particular site. The data reveal that HoxD9 binds to nonspecific DNA with the same binding mode and orientation as that observed in the specific complex. The mobility, however, of Arg side-chains contacting the DNA is increased in the nonspecific complex relative to the specific one. The kinetics of intermolecular translocation between two different nonspecific DNA molecules have also been analyzed and reveal that at high DNA concentrations (such as those present in vivo) direct transfer from one nonspecific complex to another nonspecific DNA molecule occurs without going through the intermediary of free protein. This finding provides a simple mechanism for accelerating the target search in vivo for the specific site in a sea of nonspecific sites by permitting more effective sampling of available DNA sites as the protein jumps from one segment to another.     

Three-dimensional structure of proteins determined by molecular dynamics with interproton distance restraints: application to crambin.

Model calculations are performed to evaluate the utility of molecular dynamics with NMR interproton distance restraints for determining the three-dimensional structure of proteins. The system used for testing the method is the 1.5-A resolution crystal structure of crambin (a protein of 46 residues) from which a set of 240 approximate interproton distances of less than 4 A are derived. The convergence properties of the method are examined by using different dynamics protocols and starting from two initial structures; one is a completely extended beta-strand, and the other has residues 7-19 and 23-30 in the form of alpha-helices (as in the crystal structure) with the remaining residues in the form of extended beta-strands. In both cases global and local convergence to the correct final structure is achieved with rms atomic differences between the restrained dynamics structures and the crystal structure of 1.5-2.1 A for the backbone atoms and 2.1-2.8 A for all atoms; the averaged structure has backbone and all atom rms deviations of 1.3 and 1.9 A, respectively. Further, it is shown that a restrained dynamics structure with significantly larger deviations (i.e., 5.7 A for the backbone atoms) can be characterized as incorrect, independent of a knowledge of the crystal structure.     

Paramagnetic shifts in solid-state NMR of proteins to elicit structural information

The recent observation of pseudocontact shifts (pcs) in ¹³C high-resolution solid-state NMR of paramagnetic proteins opens the way to their application as structural restraints. Here, by investigating a microcrystalline sample of cobalt(II)-substituted matrix metalloproteinase 12 [CoMMP-12 (159 AA, 17.5 kDa)], it is shown that a combined strategy of protein labeling and dilution of the paramagnetic species (i.e., ¹³C-,¹⁵N-labeled CoMMP-12 diluted in unlabeled ZnMMP-12, and ¹³C-,¹⁵N-labeled ZnMMP-12 diluted in unlabeled CoMMP-12) allows one to easily separate the pcs contributions originated from the protein internal metal (intramolecular pcs) from those due to the metals in neighboring proteins in the crystal lattice (intermolecular pcs) and that both can be used for structural purposes. It is demonstrated that intramolecular pcs are significant structural restraints helpful in increasing both precision and accuracy of the structure, which is a need in solid-state structural biology nowadays. Furthermore, intermolecular pcs provide unique information on positions and orientations of neighboring protein molecules in the solid phase.     

Determination of solution structures of proteins up to 40 kDa using CS-Rosetta with sparse NMR data from deuterated samples

We have developed an approach for determining NMR structures of proteins over 20 kDa that utilizes sparse distance restraints obtained using transverse relaxation optimized spectroscopy experiments on perdeuterated samples to guide RASREC Rosetta NMR structure calculations. The method was tested on 11 proteins ranging from 15 to 40 kDa, seven of which were previously unsolved. The RASREC Rosetta models were in good agreement with models obtained using traditional NMR methods with larger restraint sets. In five cases X-ray structures were determined or were available, allowing comparison of the accuracy of the Rosetta models and conventional NMR models. In all five cases, the Rosetta models were more similar to the X-ray structures over both the backbone and side-chain conformations than the “best effort” structures determined by conventional methods. The incorporation of sparse distance restraints into RASREC Rosetta allows routine determination of high-quality solution NMR structures for proteins up to 40 kDa, and should be broadly useful in structural biology.     

Characterizations of diverse residue clusters in protein three-dimensional structures.

We present new methods for identifying and analyzing statistically significant residue clusters that occur in three-dimensional (3D) protein structures. Residue clusters of different kinds occur in many contexts. They often feature the active site (e.g., in substrate binding), the interface between polypeptide units of protein complexes, regions of protein-protein and protein-nucleic acid interactions, or regions of metal ion coordination. The methods are illustrated with 3D clusters centering on four themes. (i) Acidic or histidine-acidic clusters associated with metal ions. (ii) Cysteine clusters including coordination of metals such as zinc or iron-sulfur structures, cysteine knots prominent in growth factors, multiple sets of buried disulfide pairings that putatively nucleate the hydrophobic core, or cysteine clusters of mostly exposed disulfide bridges. (iii) Iron-sulfur proteins and charge clusters. (iv) 3D environments of multiple histidine residues. Study of diverse 3D residue clusters offers a new perspective on protein structure and function. The algorithms can aid in rapid identification of distinctive sites, suggest correlations among protein structures, and serve as a tool in the analysis of new structures.     

Steroid diamine-nucleic acid interactions: partial insertion of dipyrandium between unstacked base pairs of the poly(dA-dT) duplex in solution.

We report on an NMR investigation of steroid diamine.nucleic acid complexes as a function of phosphate-to-drug ratios in aqueous solution in order to evaluate the structural and kinetic aspects of the binding of a nonintercalative drug to a synthetic DNA in solution. The nonexchangeable proton chemical shift parameters for the dipyrandium.poly(dA-dT) complex demonstrate unstacking of base pairs and partial insertion of the steroid diamine at the complexation site. The chemical shifts and linewidths of the exchangeable protons as a function of pH demonstrate that the base pairs are intact but partially exposed to solvent at the steroid diamine binding site. The phosphorus chemical shifts suggest that the base pairs unstack upon complex formation without changes in the omega,omega' polynucleotide backbone torsion angles. The NMR line shape parameters require rapid exchange of the steroid diamine among potential binding sites and are consistent with greater segmental flexibility in the complex compared to the synthetic DNA in solution. The NMR experiments are discussed in relation to Sobell's proposed model [Sobell, H.M., Tsai, C.C., Gilbert, S.G., Jain, S.C. & Sakore, T.D. (1976) Proc. Natl. Acad. Sci USA 73, 3068-3072] for the steroid diamine-DNA complex.     

Templates are available to model nearly all complexes of structurally characterized proteins

Traditional approaches to protein-protein docking sample the binding modes with no regard to similar experimentally determined structures (templates) of protein-protein complexes. Emerging template-based docking approaches utilize such similar complexes to determine the docking predictions. The docking problem assumes the knowledge of the participating proteins' structures. Thus, it provides the possibility of aligning the structures of the proteins and the template complexes. The progress in the development of template-based docking and the vast experience in template-based modeling of individual proteins show that, generally, such approaches are more reliable than the free modeling. The key aspect of this modeling paradigm is the availability of the templates. The current common perception is that due to the difficulties in experimental structure determination of protein-protein complexes, the pool of docking templates is insignificant, and thus a broad application of template-based docking is possible only at some future time. The results of our large scale, systematic study show that, surprisingly, in spite of the limited number of protein-protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the known protein-protein interactions, provided the components themselves have a known structure or can be homology-built. About one-third of the templates are of good quality when they are compared to experimental structures in test sets extracted from the Protein Data Bank and would be useful starting points in modeling the complexes. This finding dramatically expands our ability to model protein interactions, and has far-reaching implications for the protein docking field in general.     

Rapid protein fold determination using unassigned NMR data

Experimental structure determination by x-ray crystallography and NMR spectroscopy is slow and time-consuming compared with the rate at which new protein sequences are being identified. NMR spectroscopy has the advantage of rapidly providing the structurally relevant information in the form of unassigned chemical shifts (CSs), intensities of NOESY crosspeaks [nuclear Overhauser effects (NOEs)], and residual dipolar couplings (RDCs), but use of these data are limited by the time and effort needed to assign individual resonances to specific atoms. Here, we develop a method for generating low-resolution protein structures by using unassigned NMR data that relies on the de novo protein structure prediction algorithm, rosetta [Simons, K. T., Kooperberg, C., Huang, E. & Baker, D. (1997) J. Mol. Biol. 268, 209-225] and a Monte Carlo procedure that searches for the assignment of resonances to atoms that produces the best fit of the experimental NMR data to a candidate 3D structure. A large ensemble of models is generated from sequence information alone by using rosetta, an optimal assignment is identified for each model, and the models are then ranked based on their fit with the NMR data assuming the identified assignments. The method was tested on nine protein sequences between 56 and 140 amino acids and published CS, NOE, and RDC data. The procedure yielded models with rms deviations between 3 and 6 A, and, in four of the nine cases, the partial assignments obtained by the method could be used to refine the structures to high resolution (0.6-1.8 A) by repeated cycles of structure generation guided by the partial assignments, followed by reassignment using the newly generated models.     

Insufficiently dehydrated hydrogen bonds as determinants of protein interactions

The prediction of binding sites and the understanding of interfaces associated with protein complexation remains an open problem in molecular biophysics. This work shows that a crucial factor in predicting and rationalizing protein-protein interfaces can be inferred by assessing the extent of intramolecular desolvation of backbone hydrogen bonds in monomeric structures. Our statistical analysis of native structures shows that, in the majority of soluble proteins, most backbone hydrogen bonds are thoroughly wrapped intramolecularly by nonpolar groups except for a few ones. These latter underwrapped hydrogen bonds may be dramatically stabilized by removal of water. This fact implies that packing defects are "sticky" in a way that decisively contributes to determining the binding sites for proteins, as an examination of numerous complexes demonstrates.     

High-resolution molecular structure of a peptide in an amyloid fibril determined by magic angle spinning NMR spectroscopy

Amyloid fibrils are self-assembled filamentous structures associated with protein deposition conditions including Alzheimer's disease and the transmissible spongiform encephalopathies. Despite the immense medical importance of amyloid fibrils, no atomic-resolution structures are available for these materials, because the intact fibrils are insoluble and do not form diffraction-quality 3D crystals. Here we report the high-resolution structure of a peptide fragment of the amyloidogenic protein transthyretin, TTR(105-115), in its fibrillar form, determined by magic angle spinning NMR spectroscopy. The structure resolves not only the backbone fold but also the precise conformation of the side chains. Nearly complete (13)C and (15)N resonance assignments for TTR(105-115) formed the basis for the extraction of a set of distance and dihedral angle restraints. A total of 76 self-consistent experimental measurements, including 41 restraints on 19 backbone dihedral angles and 35 (13)C-(15)N distances between 3 and 6 A were obtained from 2D and 3D NMR spectra recorded on three fibril samples uniformly (13)C, (15)N-labeled in consecutive stretches of four amino acids and used to calculate an ensemble of peptide structures. Our results indicate that TTR(105-115) adopts an extended beta-strand conformation in the amyloid fibrils such that both the main- and side-chain torsion angles are close to their optimal values. Moreover, the structure of this peptide in the fibrillar form has a degree of long-range order that is generally associated only with crystalline materials. These findings provide an explanation of the unusual stability and characteristic properties of this form of polypeptide assembly.     

Identification of slow correlated motions in proteins using residual dipolar and hydrogen-bond scalar couplings

Despite their importance for biological activity, slower molecular motions beyond the nanosecond range remain poorly understood. We have assembled an unprecedented set of experimental NMR data, comprising up to 27 residual dipolar couplings per amino acid, to define the nature and amplitude of backbone motion in protein G using the Gaussian axial fluctuation model in three dimensions. Slower motions occur in the loops, and in the beta-sheet, and are absent in other regions of the molecule, including the alpha-helix. In the beta-sheet an alternating pattern of dynamics along the peptide sequence is found to form a long-range network of slow motion in the form of a standing wave extending across the beta-sheet, resulting in maximal conformational sampling at the interaction site. The alternating nodes along the sequence match the alternation of strongly hydrophobic side chains buried in the protein core. Confirmation of the motion is provided through extensive cross-validation and by independent hydrogen-bond scalar coupling analysis that shows this motion to be correlated. These observations strongly suggest that dynamical information can be transmitted across hydrogen bonds and have important implications for understanding collective motions and long-range information transfer in proteins.     

Consistent blind protein structure generation from NMR chemical shift data

Protein NMR chemical shifts are highly sensitive to local structure. A robust protocol is described that exploits this relation for de novo protein structure generation, using as input experimental parameters the (13)C(alpha), (13)C(beta), (13)C', (15)N, (1)H(alpha) and (1)H(N) NMR chemical shifts. These shifts are generally available at the early stage of the traditional NMR structure determination process, before the collection and analysis of structural restraints. The chemical shift based structure determination protocol uses an empirically optimized procedure to select protein fragments from the Protein Data Bank, in conjunction with the standard ROSETTA Monte Carlo assembly and relaxation methods. Evaluation of 16 proteins, varying in size from 56 to 129 residues, yielded full-atom models that have 0.7-1.8 A root mean square deviations for the backbone atoms relative to the experimentally determined x-ray or NMR structures. The strategy also has been successfully applied in a blind manner to nine protein targets with molecular masses up to 15.4 kDa, whose conventional NMR structure determination was conducted in parallel by the Northeast Structural Genomics Consortium. This protocol potentially provides a new direction for high-throughput NMR structure determination.