Structures of protein-protein complexes are docked using only NMR restraints from residual dipolar coupling and chemical shift perturbations

NMR structures of protein-protein and protein-ligand complexes rely heavily on intermolecular NOEs. Recent work has shown that if no significant conformational changes occur upon complex formation residual dipolar coupling can replace most of the NOE restraints in protein-protein complexes, while restraints derived from chemical shift perturbations can largely replace intermolecular NOEs in protein-ligand structures. By combining restraints from chemical shift perturbations with orientation restraints derived from measurements of residual dipolar couplings, we show that the structure of the EIN-HPr complex can be calculated without NOE restraints. The final structure, built from the crystal structures of EIN and HPr in their uncomplexed form and docked only with NMR restraints, places HPr within 2.5 A of the position determined from the mean NMR structure of the complex.     

Impact of 15N R2/R1 relaxation restraints on molecular size, shape, and bond vector orientation for NMR protein structure determination with sparse distance restraints

(15)N R(2)/R(1) relaxation data contain information on molecular shape and size as well as on bond vector orientations relative to the diffusion tensor. Since the diffusion tensor can be directly calculated from the molecular coordinates, direct inclusion of (15)N R(2)/R(1) restraints in NMR structure calculations without any a priori assumptions is possible. Here we show that (15)N R(2)/R(1) restraints are particularly valuable when only sparse distance restraints are available. Using three examples of proteins of varying size, namely, GB3 (56 residues), ubiquitin (76 residues), and the N-terminal domain of enzyme I (EIN, 249 residues), we show that incorporation of (15)N R(2)/R(1) restraints results in large and significant increases in coordinate accuracy that can make the difference between being able or unable to determine an approximate global fold. For GB3 and ubiquitin, good coordinate accuracy was obtained using only backbone hydrogen-bond restraints supplemented by (15)N R(2)/R(1) relaxation restraints. For EIN, the global fold could be determined using sparse nuclear Overhauser enhancement (NOE) distance restraints involving only NH and methyl groups in conjunction with (15)N R(2)/R(1) restraints. These results are of practical significance in the study of larger and more complex systems, where the increasing spectral complexity and number of chemical shift degeneracies reduce the number of unambiguous NOE assignments that can be readily obtained, resulting in progressively reduced NOE coverage as the size of the protein increases.     

Paramagnetic-based NMR restraints lift residual dipolar coupling degeneracy in multidomain detergent-solubilized membrane proteins

Residual dipolar couplings (RDCs) are widely used as orientation-dependent NMR restraints to improve the resolution of the NMR conformational ensemble of biomacromolecules and define the relative orientation of multidomain proteins and protein complexes. However, the interpretation of RDCs is complicated by the intrinsic degeneracy of analytical solutions and protein dynamics that lead to ill-defined orientations of the structural domains (ghost orientations). Here, we illustrate how restraints from paramagnetic relaxation enhancement (PRE) experiments lift the orientational ambiguity of multidomain membrane proteins solubilized in detergent micelles. We tested this approach on monomeric phospholamban (PLN), a 52-residue membrane protein, which is composed of two helical domains connected by a flexible loop. We show that the combination of classical solution NMR restraints (NOEs and dihedral angles) with RDC and PRE constraints resolves topological ambiguities, improving the convergence of the PLN structural ensemble and giving the depth of insertion of the protein within the micelle. The combination of RDCs with PREs will be necessary for improving the accuracy and precision of membrane protein conformational ensembles, where three-dimensional structures are dictated by interactions with the membrane-mimicking environment rather than compact tertiary folds common in globular proteins.     

Direct use of unassigned resonances in NMR structure calculations with proxy residues

We present a method that significantly enhances the robustness of (automated) NMR structure determination by allowing the NOE data corresponding to unassigned NMR resonances to be used directly in the calculations. The unassigned resonances are represented by additional atoms or groups of atoms that have no interaction with the regular protein atoms except through distance restraints. These so-called "proxy" residues can be used to generate NOE-based distance restraints in a similar fashion as for the assigned part of the protein. If sufficient NOE information is available, the restraints are expected to place the proxies at positions close to the correct atoms for the unassigned resonance, which can facilitate subsequent assignment. Convergence can be further improved by supplying additional information about the possible identities of the unassigned resonances. We have implemented this approach in the widely used automated assignment and structure calculation protocols ARIA and CANDID. We find that it significantly increases the robustness of structure calculations with regard to missing assignments and yields structures of higher quality. Our approach is still able to find correctly folded structures with up to 30% randomly missing resonance assignments, and even when only backbone and beta resonances are present! This should be of significant value to NMR-based structural proteomics initiatives.     

Increasing the accuracy of solution NMR structures of membrane proteins by application of residual dipolar couplings. High-resolution structure of outer membrane protein A

The structure determination of membrane proteins is one of the most challenging applications of solution NMR spectroscopy. The paucity of distance information available from the highly deuterated proteins employed requires new approaches in structure determination. Here we demonstrate that significant improvement in the structure accuracy of the membrane protein OmpA can be achieved by refinement with residual dipolar couplings (RDCs). The application of charged polyacrylamide gels allowed us to obtain two alignments and accurately measure numerous heteronuclear dipolar couplings. Furthermore, we have demonstrated that using a large set of RDCs in the refinement can yield a structure with 1 A rms deviation to the backbone of the high-resolution crystal structure. Our simulations with various data sets indicate that dipolar couplings will be critical for obtaining accurate structures of membrane proteins.     

Ensemble approach for NMR structure refinement against (1)H paramagnetic relaxation enhancement data arising from a flexible paramagnetic group attached to a macromolecule

Paramagnetic relaxation enhancement (PRE) measurements on (1)H nuclei have the potential to play an important role in NMR structure determination of macromolecules by providing unique long-range (10-35 A) distance information. Recent methodological advances for covalently attaching paramagnetic groups at specific sites on both proteins and nucleic acids have permitted the application of the PRE to various biological macromolecules. However, because artificially introduced paramagnetic groups are exposed to solvent and linked to the macromolecule by several freely rotatable bonds, they are intrinsically flexible. This renders conventional back-calculation of the (1)H-PRE using a single-point representation inaccurate, thereby severely limiting the utility of the (1)H-PRE as a tool for structure refinement. To circumvent these limitations, we have developed a theoretical framework and computational strategy with which to accurately back-calculate (1)H-PREs arising from flexible paramagnetic groups attached to macromolecules. In this scheme, the (1)H-PRE is calculated using a modified Solomon-Bloembergen equation incorporating a "model-free" formalism, based on a multiple-structure representation of the paramagnetic group in simulated annealing calculations. The ensemble approach for (1)H-PRE back-calculation was examined using several SRY/DNA complexes incorporating dT-EDTA-Mn(2+) at three distinct sites in the DNA, permitting a large data set comprising 435 experimental backbone and side-chain (1)H-PREs to be obtained in a straightforward manner from 2D through-bond correlation experiments. Calculations employing complete cross-validation demonstrate that the ensemble representation provides a means to accurately utilize backbone and side-chain (1)H-PRE data arising from a flexible paramagnetic group in structure refinement. The results of (1)H-PRE based refinement, in conjunction with previously obtained NMR restraints, indicate that significant gains in accuracy can be readily obtained. This is particularly significant in the case of macromolecular complexes where intermolecular translational restraints derived from nuclear Overhauser enhancement data may be limited.     

Completely automated, highly error-tolerant macromolecular structure determination from multidimensional nuclear overhauser enhancement spectra and chemical shift assignments

The major rate-limiting step in high-throughput NMR protein structure determination involves the calculation of a reliable initial fold, the elimination of incorrect nuclear Overhauser enhancement (NOE) assignments, and the resolution of NOE assignment ambiguities. We present a robust approach to automatically calculate structures with a backbone coordinate accuracy of 1.0-1.5 A from datasets in which as much as 80% of the long-range NOE information (i.e., between residues separated by more than five positions in the sequence) is incorrect. The current algorithm differs from previously published methods in that it has been expressly designed to ensure that the results from successive cycles are not biased by the global fold of structures generated in preceding cycles. Consequently, the method is highly error tolerant and is not easily funnelled down an incorrect path in either three-dimensional structure or NOE assignment space. The algorithm incorporates three main features: a linear energy function representation of the NOE restraints to allow maximization of the number of simultaneously satisfied restraints during the course of simulated annealing; a method for handling the presence of multiple possible assignments for each NOE cross-peak which avoids local minima by treating each possible assignment as if it were an independent restraint; and a probabilistic method to permit both inactivation and reactivation of all NOE restraints on the fly during the course of simulated annealing. NOE restraints are never removed permanently, thereby significantly reducing the likelihood of becoming trapped in a false minimum of NOE assignment space. The effectiveness of the algorithm is demonstrated using completely automatically peak-picked experimental NOE data from two proteins: interleukin-4 (136 residues) and cyanovirin-N (101 residues). The limits of the method are explored using simulated data on the 56-residue B1 domain of Streptococcal protein G.     

Structural characterization of proteins with an attached ATCUN motif by paramagnetic relaxation enhancement NMR spectroscopy

The use of a short, three-residue Cu(2+)-binding sequence, the ATCUN motif, is presented as an approach for extracting long-range distance restraints from relaxation enhancement NMR spectroscopy. The ATCUN motif is prepended to the N-termini of proteins and binds Cu(2+) with a very high affinity. Relaxation rates of amide protons in ATCUN-tagged protein in the presence and absence of Cu(2+) can be converted into distance restraints and used for structure refinement by using a new routine, PMAG, that has been written for the structure calculation program CNS. The utility of the approach is demonstrated with an application to ATCUN-tagged ubiquitin. Excellent agreement between measured relaxation rates and those calculated on the basis of the X-ray structure of the protein have been obtained.     

Interpreting NMR data for beta-peptides using molecular dynamics simulations

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.     

Ultraviolet resonance Raman spectroscopy of folded and unfolded states of an integral membrane protein

The vibrational structure of native anchoring tryptophan (Trp) and tyrosine residues in an integral membrane protein, bacterial outer membrane protein A (OmpA), have been investigated using UV resonance Raman (UVRR) spectroscopy for the first time. Spectra of native OmpA, a single-Trp mutant, and a Trp-less mutant were recorded in folded and unfolded states, and reveal significant changes in tryptophan structure and local environment. Salient alterations upon folding include loss of hydrogen-bonding character of indole N1H, evidenced by a shift in W17 frequency from 874 and 878 cm(-1), and growth in hydrophobicity of the local tryptophan environment, supported by increase in the ratio I1361/I1340. In addition to these site-specific changes in a single tryptophan residue, modification of the vibrational structure of the remaining native tryptophan and tyrosine amino acids is also evident. Finally, the UVRR data presented here indicate that the structures of OmpA folded in vesicle and folded in detergent may differ, and provide important foundations for ongoing studies of membrane protein folding.     

Charged gels as orienting media for measurement of residual dipolar couplings in soluble and integral membrane proteins

Measurement of residual dipolar couplings for membrane proteins will dramatically improve the quality of the structures obtainable by solution NMR spectroscopy. While there has been some success in achieving alignment of membrane-bound peptides, there has been very limited success in achieving alignment for functional membrane proteins. Herein, we demonstrate that charged polyacrylamide-based copolymers are suitable for obtaining weak alignment of membrane proteins reconstituted in detergent micelles. Varying the copolymer compositions, we prepared positively, zwitterionic, and negatively charged gels that are very stable at low concentration and can be used for obtaining weak alignment by compression in an NMR tube. Application of this method is demonstrated for the integral membrane protein OmpA in DPC micelles.     

Structure determination of a key intermediate of the enantioselective Pd complex catalyzed allylic substitution reaction

The structure of a catalytic intermediate with important implications for the interpretation of the stereochemical outcome of the palladium complex catalyzed allylic substitution with phosphino-oxazoline (PHOX) ligands is determined by liquid state NMR. The complex displays a novel structure that is highly distorted compared with other palladium eta2-olefin complexes known so far. The structure has been determined from nuclear overhauser data (NOE), scalar coupling constants, and long range projection angle restraints derived from dipole dipole cross-correlated relaxation of multiple quantum coherence. The latter restraints have been implemented into a distance geometry protocol. The projection angle restraints yield a higher precision in the determination of the relative orientation of the two molecular moieties and are essential to provide an exact structural definition of the olefinic part of the catalytic intermediate with respect to the ligand.     

Bioconjugation of Proteins with a Paramagnetic NMR and Fluorescent Tag

Site‐specific labeling of proteins with lanthanide ions offers great opportunities for investigating the structure, function, and dynamics of proteins by virtue of the unique properties of lanthanides. Lanthanide‐tagged proteins can be studied by NMR, X‐ray, fluorescence, and EPR spectroscopy. However, the rigidity of a lanthanide tag in labeling of proteins plays a key role in the determination of protein structures and interactions. Pseudocontact shift (PCS) and paramagnetic relaxation enhancement (PRE) are valuable long‐range structure restraints in structural‐biology NMR spectroscopy. Generation of these paramagnetic restraints generally relies on site‐specific tagging of the target proteins with paramagnetic species. To avoid nonspecific interaction between the target protein and paramagnetic tag and achieve reliable paramagnetic effects, the rigidity, stability, and size of lanthanide tag is highly important in paramagnetic labeling of proteins. Here 4′‐mercapto‐2,2′: 6′,2′′‐terpyridine‐6,6′′‐dicarboxylic acid (4MTDA) is introduced as a a rigid paramagnetic and fluorescent tag which can be site‐specifically attached to a protein by formation of a disulfide bond. 4MTDA can be readily immobilized by coordination of the protein side chain to the lanthanide ion. Large PCSs and RDCs were observed for 4MTDA‐tagged proteins in complexes with paramagnetic lanthanide ions. At an excitation wavelength of 340 nm, the complex formed by protein–4MTDA and Tb³⁺ produces high fluorescence with the main emission at 545 nm. These interesting features of 4MTDA make it a very promising tag that can be exploited in NMR, fluorescence, and EPR spectroscopic studies on protein structure, interaction, and dynamics.     

Stereospecific assignments of protein NMR resonances based on the tertiary structure and 2D/3D NOE data

In many cases of protein structure determination by NMR a high-quality structure is required. An important contribution to structural precision is stereospecific assignment of magnetically nonequivalent prochiral methylene and methyl groups, eliminating the need for introducing pseudoatoms and pseudoatom corrections in distance restraint lists. Here, we introduce the stereospecific assignment program that uses the resonance assignment, a preliminary 3D structure and 2D and/or 3D nuclear Overhauser effect spectroscopy peak lists for stereospecific assignment. For each prochiral group the algorithm automatically calculates a score for the two different stereospecific assignment possibilities, taking into account the presence and intensity of the nuclear Overhauser effect (NOE) peaks that are expected from the local environment of each prochiral group (i.e., the close neighbors). The performance of the algorithm has been tested and used on NMR data of alpha-helical and beta-sheet proteins using homology models and/or X-ray structures. The program produced no erroneous stereospecific assignments provided the NOEs were carefully picked and the 3D model was sufficiently accurate. The set of NOE distance restraints produced by nmr2st using the results of the SSA module was superior in generating good-quality ensembles of NMR structures (low deviations from upper limits in conjunction with low root-mean-square-deviation values) in the first round of structure calculations. The program uses a novel approach that employs the entire 3D structure of the protein to obtain stereospecific assignment; it can be used to speed up the NMR structure refinement and to increase the quality of the final NMR ensemble even when no scalar or residual dipolar coupling information is available.     

Conformation and absolute configuration of nocathiacin I determined by NMR spectroscopy and chiral capillary electrophoresis

Nocathiacin I (BMS-249524) is a highly cross-linked thiazolyl peptide that displays potent activity against Gram-positive bacteria, including a number of antibiotic-resistant strains. This natural product contains 10 chiral centers. NMR studies have been performed to characterize the solution structure of nocathiacin I. A uniformly 13C,15N-labeled sample was used to obtain NMR assignments. Restrained simulated annealing calculations were performed by using accurately determined NOE distance restraints. All of the chiral centers were allowed to float during the simulated annealing protocol. Two clusters of structures were obtained that satisfy the NOE restraints very well and that are reasonably consistent with vicinal J-coupling constants. Within each cluster, all 10 chiral centers are uniquely defined. The two clusters are effectively mirror images of each other: all chiral centers that have the R(S) configuration in one cluster have the S(R) configuration in the other. The single threonine residue in nocathiacin I was subsequently determined to be l-threonine by chiral capillary electrophoresis, allowing the absolute configurations of all 10 chiral centers to be defined.     

3D structure determination of the Crh protein from highly ambiguous solid-state NMR restraints

In a wide variety of proteins, insolubility presents a challenge to structural biology, as X-ray crystallography and liquid-state NMR are unsuitable. Indeed, no general approach is available as of today for studying the three-dimensional structures of membrane proteins and protein fibrils. We here demonstrate, at the example of the microcrystalline model protein Crh, how high-resolution 3D structures can be derived from magic-angle spinning solid-state NMR distance restraints for fully labeled protein samples. First, we show that proton-mediated rare-spin correlation spectra, as well as carbon-13 spin diffusion experiments, provide enough short, medium, and long-range structural restraints to obtain high-resolution structures of this 2 x 10.4 kDa dimeric protein. Nevertheless, the large number of 13C/15N spins present in this protein, combined with solid-state NMR line widths of about 0.5-1 ppm, induces substantial ambiguities in resonance assignments, preventing 3D structure determination by using distance restraints uniquely assigned on the basis of their chemical shifts. In the second part, we thus demonstrate that an automated iterative assignment algorithm implemented in a dedicated solid-state NMR version of the program ARIA permits to resolve the majority of ambiguities and to calculate a de novo 3D structure from highly ambiguous solid-state NMR data, using a unique fully labeled protein sample. We present, using distance restraints obtained through the iterative assignment process, as well as dihedral angle restraints predicted from chemical shifts, the 3D structure of the fully labeled Crh dimer refined at a root-mean-square deviation of 1.33 A.     

Hydrogen exchange study on the hydroxyl groups of serine and threonine residues in proteins and structure refinement using NOE restraints with polar side-chain groups

We recently developed new NMR methods for monitoring the hydrogen exchange rates of tyrosine hydroxyl (Tyr-OH) and cysteine sulfhydryl (Cys-SH) groups in proteins. These methods facilitate the identification of slowly exchanging polar side-chain protons in proteins, which serve as sources of NOE restraints for protein structure refinement. Here, we have extended the methods for monitoring the hydrogen exchange rates of the OH groups of serine (Ser) and threonine (Thr) residues in an 18.2 kDa protein, EPPIb, and thus demonstrated the usefulness of NOE restraints with slowly exchanging OH protons for refining the protein structure. The slowly exchanging Ser/Thr-OH groups were readily identified by monitoring the (13)C(β)-NMR signals in an H(2)O/D(2)O (1:1) mixture, for the protein containing Ser/Thr residues with (13)C, (2)H-double labels at their β carbons. Under these circumstances, the OH groups exist in equilibrium between the protonated and deuterated isotopomers, and the (13)C(β) peaks of the two species are resolved when their exchange rate is slower than the time scale of the isotope shift effect. In the case of EPPIb dissolved in 50 mM sodium phosphate buffer (pH 7.5) at 40 °C, one Ser and four Thr residues were found to have slowly exchanging hydroxyl groups (k(ex) < ~40 s(-1)). With the information for the slowly exchanging Ser/Thr-OH groups in hand, we could collect additional NOE restraints for EPPIb, thereby making a unique and important contribution toward defining the spatial positions of the OH protons, and thus the hydrogen-bonding acceptor atoms.     

Modeling errors in NOE data with a log-normal distribution improves the quality of NMR structures

The distribution of the deviation of calculated from measured nuclear Overhauser effect (NOE) intensities is a priori unknown. The use of a log-normal distribution to describe these deviations permits the direct calculation of a structure from the measured intensities without first converting them into distance bounds. We show that the log-normal distribution is a natural choice for describing errors in NOE data and that it improves the accuracy, precision, and quality of the calculated structures compared to the usual bounds representation.     

Protein structure determination using long-distance constraints from double-quantum coherence ESR: study of T4 lysozyme

We report the use of a novel pulsed ESR technique for distance measurement, based on the detection of double quantum coherence (DQC), which yields high quality dipolar spectra, to significantly extend the range of measurable distances in proteins using nitroxide spin-labels. Eight T4 lysozyme (T4L) mutants, doubly labeled with methanethiosulfonate spin-label (MTSSL), have been studied using DQC-ESR at 9 and 17 GHz. The distances span the range from 20 A for the 65/76 mutant to 47 A for the 61/135 mutant. The high quality of the dipolar spectra also allows the determination of the distance distributions, the width of which can be used to set upper and lower bounds in future computational strategy. It is also demonstrated that the shape of these distributions can reveal the presence of multiple conformations of the spin-label, an issue of critical relevance to the structural interpretation of the distances. The distances and distributions found in this study are readily rationalized in terms of the known crystal structure, the characteristic conformers of the nitroxide side chains, and molecular modeling. This study sets the stage for the use of DQC-ESR for determining the tertiary structure of large proteins with just a small number of long-distance constraints.