Measurement of (13)C(alpha)-(13)C(beta) dipolar couplings in (15)N,(13)C,(2)H-labeled proteins: application to domain orientation in maltose binding protein

Measurement of (1)H-(15)N and (1)H-(13)C residual dipolar couplings in nucleic acids from TROSY intensities

Analogous to the recently introduced ARTSY method for measurement of one-bond (1)H-(15)N residual dipolar couplings (RDCs) in large perdeuterated proteins, we introduce methods for measurement of base (13)C-(1)H and (15)N-(1)H RDCs in protonated nucleic acids. Measurements are based on quantitative analysis of intensities in (1)H-(15)N and (13)C-(1)H TROSY-HSQC spectra, and are illustrated for a 71-nucleotide adenine riboswitch. Results compare favorably with those of conventional frequency-based measurements in terms of completeness and convenience of use. The ARTSY method derives the size of the coupling from the ratio of intensities observed in two TROSY-HSQC spectra recorded with different dephasing delays, thereby minimizing potential resonance overlap problems. Precision of the RDC measurements is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC reference spectrum, and is approximately given by 30/(S/N) Hz for (15)N-(1)H and 65/(S/N) Hz for (13)C-(1)H. The signal-to-noise ratio of both (1)H-(15)N and (1)H-(13)C spectra greatly benefits when water magnetization during the experiments is not perturbed, such that rapid magnetization transfer from bulk water to the nucleic acid, mediated by rapid amino and hydroxyl hydrogen exchange coupled with (1)H-(1)H NOE transfer, allows for fast repetition of the experiment. RDCs in the mutated helix 1 of the riboswitch are compatible with nucleotide-specifically modeled, idealized A-form geometry and a static orientation relative to the helix 2/3 pair, which differs by ca 6° relative to the X-ray structure of the native riboswitch.

One-bond and two-bond J couplings help annotate protein secondary-structure motifs: J-coupling indexing applied to human endoplasmic reticulum protein ERp18

NMR coupling constants, both direct one-bond ((1)J) and geminal two-bond ((2)J), are employed to analyze the protein secondary structure of human oxidized ERp18. Coupling constants collected and evaluated for the 18 kDa protein comprise 1268 values of (1)J(CαHα), (1)J(CαCβ), (1)J(CαC'), (1)J(C'N'), (1)J(N'Cα), (1)J(N') (HN), (2)J(CαN'), (2)J(HNCα), (2)J(C'HN), and (2)J(HαC'). Comparison with (1)J and (2)J data from reference proteins and pattern analysis on a per-residue basis permitted main-chain ϕ,ψ torsion-angle combinations of many of the 149 amino-acid residues in ERp18 to be narrowed to particular secondary-structure motifs. J-coupling indexing is here being developed on statistical criteria and used to devise a ternary grid for interpreting patterns of relative values of J. To account for the influence of the varying substituent pattern in different amino-acid sidechains, a table of residue-type specific threshold values was compiled for discriminating small, medium, and large categories of J. For the 15-residue insertion that distinguishes the ERp18 fold from that of thioredoxin, the J-coupling data hint at a succession of five isolated Type-I β turns at progressively shorter sequence intervals, in agreement with the crystal structure.

A single destabilizing mutation (F9S) promotes concerted unfolding of an entire globular domain in gammaS-crystallin

Conformational change and aggregation of native proteins are associated with many serious age-related and neurological diseases. gammaS-Crystallin is a highly stable, abundant structural component of vertebrate eye lens. A single F9S mutation in the N-terminal domain of mouse gammaS-crystallin causes the severe Opj cataract, with disruption of cellular organization and appearance of fibrillar structures in the lens. Although the mutant protein has a near-native fold at room temperature, significant increases in hydrogen/deuterium exchange rates were observed by NMR for all the well-protected beta-sheet core residues throughout the entire N-terminal domain of the mutant protein, resulting in up to a 3.5-kcal/mol reduction in the free energy of the folding/unfolding equilibrium. No difference was detected for the C-terminal domain. At a higher temperature, this effect further increases to allow for a much more uniform exchange rate among the N-terminal core residues and those of the least well-structured surface loops. This suggests a concerted unfolding intermediate of the N-terminal domain, while the C-terminal domain stays intact. Increasing concentrations of guanidinium chloride produced two transitions for the Opj mutant, with an unfolding intermediate at approximately 1 M guanidinium chloride. The consequence of this partial unfolding, whether by elevated temperature or by denaturant, is the formation of thioflavin T staining aggregates, which demonstrated fibril-like morphology by atomic force microscopy. Seeding with the already unfolded protein enhanced the formation of fibrils. The Opj mutant protein provides a model for stress-related unfolding of an essentially normally folded protein and production of aggregates with some of the characteristics of amyloid fibrils.

Variation in protein C(alpha)-related one-bond J couplings

Four types of polypeptide (1)J(C alpha X) couplings are examined, involving the main-chain carbon C(alpha) and either of four possible substituents. A total 3105 values of (1)J(C alpha H alpha), (1)J(C alpha C beta), (1)J(C alpha C'), and (1)J(C alpha N') were collected from six proteins, averaging 143.4 +/- 3.3, 34.9 +/- 2.5, 52.6 +/- 0.9, and 10.7 +/- 1.2 Hz, respectively. Analysis of variances (ANOVA) reveals a variety of factors impacting on (1)J and ranks their relative statistical significance and importance to biomolecular NMR structure refinement. Accordingly, the spread in the (1)J values is attributed, in equal proportions, to amino-acid specific substituent patterns and to polypeptide-chain geometry, specifically torsions phi, psi, and chi(1) circumjacent to C(alpha). The (1)J coupling constants correlate with protein secondary structure. For alpha-helical phi, psi combinations, (1)J(C alpha H alpha) is elevated by more than one standard deviation (147.8 Hz), while both (1)J(C alpha N') and (1)J(C alpha C beta) fall short of their grand means (9.5 and 33.7 Hz). Rare positive phi torsion angles in proteins exhibit concomitant small (1)J(C alpha H alpha) and (1)J(C alpha N') (138.4 and 9.6 Hz) and large (1)J(C alpha C beta) (39.9 Hz) values. The (1)J(C alpha N') coupling varies monotonously over the phi torsion range typical of beta-sheet secondary structure and is largest (13.3 Hz) for phi around -160 degrees. All four coupling types depend on psi and thus help determine a torsion that is notoriously difficult to assess by traditional approaches using (3)J. Influences on (1)J stemming from protein secondary structure and other factors, such as amino-acid composition, are largely independent.

Conformational analysis of a tetrasaccharide based on NMR spectroscopy and molecular dynamics simulations

The conformational preference of the human milk oligosaccharide lacto-N-neotetraose, beta-d-Galp-(1 --> 4)-beta-d-GlcpNAc-(1 --> 3)-beta-d-Galp-(1 --> 4)-d-Glcp, has been analyzed using (1)H,(1)H T-ROESY and (1)H,(13)C trans-glycosidic J coupling experiments in isotropic solution and (1)H,(13)C residual dipolar couplings (RDCs) obtained in lyotropic liquid crystalline media. Molecular dynamics simulations of the tetrasaccharide with explicit water as the solvent revealed that two conformational states are significantly populated at the psi glycosidic torsion angle, defined by C(anomeric)-O-C-H, of the (1 --> 3)-linkage. Calculation of order parameters, related to the molecular shape, were based on the inertia tensor and fitting of experimental RDCs to different conformational states showed that psi(+) > 0 degrees is the major and psi(-) < 0 degrees is the minor conformation in solution, in complete agreement with a two-state analysis based on the T-ROESY data. Attention was also given to the effect of salt (200 mM NaCl) in the anisotropic medium, which was a ternary mixture of n-octyl-penta(ethylene glycol), n-octanol, and D(2)O.

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.

Solution structure of (gamma)S-crystallin by molecular fragment replacement NMR

The solution structure of murine gammaS-crystallin (gammaS) has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of dipolar couplings, recorded in different alignment media, and supplemented by a small number of NOE distance restraints. gammaS consists of two topologically similar domains, arranged with an approximate twofold symmetry, and each domain shows close structural homology to closely related (approximately 50% sequence identity) domains found in other members of the gamma-crystallin family. Each domain consists of two four-strand "Greek key" beta-sheets. Although the domains are tightly anchored to one another by the hydrophobic surfaces of the two inner Greek key motifs, the N-arm, the interdomain linker and several turn regions show unexpected flexibility and disorder in solution. This may contribute entropic stabilization to the protein in solution, but may also indicate nucleation sites for unfolding or other structural transitions. The method used for solving the gammaS structure relies on the recently introduced molecular fragment replacement method, which capitalizes on the large database of protein structures previously solved by X-ray crystallography and NMR.

Structural Comparison of the Unstable drkN SH3 Domain and a Stable Mutant†,‡

The N-terminal SH3 domain of the Drosophila adapter protein Drk (drkN SH3 domain) is marginally stable (DeltaG(U) = 1 kcal/mol) and exists in equilibrium between folded and highly populated unfolded states. The single substitution T22G, however, completely stabilizes the protein (DeltaG(U) = 4.0 kcal/mol). To probe the causes of instability of the wild-type (WT) protein and the dramatic stabilization of the mutant, we determined and compared nuclear magnetic resonance structures of the folded WT and mutant drkN SH3 domains. Residual dipolar coupling (RDC) and carbonyl chemical-shift anisotropy (C'-CSA) restraints measured for the WT and T22G domains were used for calculating the structures. The structures for the WT and mutant are highly similar. Thr22 of the WT and Gly22 of the mutant are at the i + 2 position of the diverging, type-II beta-turn. Interestingly, not only Gly22 but also Thr22 successfully adopt an alpha(L) conformation, required at this position of the turn, despite the fact that positive phi values are energetically unfavorable and normally disallowed for threonine residues. Forcing the Thr22 residue into this unnatural conformation increases the free energy of the folded state of the WT domain relative to its T22G mutant. Evidence for residual helix formation in the diverging turn region has been previously reported for the unfolded state of the WT drkN SH3 domain, and this, in addition to other residual structure, has been proposed to play a role in decreasing the free energy of the unfolded state of the protein. Together these data provide evidence that both increasing the free energy of the folded state and decreasing the free energy of the unfolded state of the protein contribute to instability of the WT drkN SH3 domain.

Accurate determination of leucine and valine side-chain conformations using U-[15N/13C/2H]/[1H-(methine/methyl)-Leu/Val] isotope labeling, NOE pattern recognition, and methine Cgamma-Hgamma/Cbeta-Hbeta residual dipolar couplings: application to the 34-kDa enzyme IIA(chitobiose)

An isotope labeling scheme is described in which specific protonation of methine and methyl protons of leucine and valine is obtained on a 15N/13C labeled background with uniform deuteration of all other non-exchangeable protons. The presence of a protonated methine group has little effect on the favorable relaxation properties of the methyl protons of Leu and Val. This labeling scheme permits the rotameric state of leucine side-chains to be readily determined by simple inspection of the pattern of Hgamma(i)-H(N)(i) and Hgamma(i)-H(N)(i+1) NOEs in a 3D 15N-separated NOE spectrum free of complications arising from spectral overlap and spin-diffusion. In addition, one-bond residual dipolar couplings for the methine 13C-1H bond vectors of Leu and Val can be accurately determined from an intensity J-modulated constant-time HCCH-COSY experiment and used to accurately orient the side-chains of Leu and Val. Incorporation of these data into structure refinement improves the accuracy with which the conformations of Leu and Val side-chains can be established. This is important to ensure optimal packing both within the protein core and at intermolecular interfaces. The impact of the method on protein structure determination is illustrated by application to enzyme IIA(Chitobiose), a 34 kDa homotrimeric phosphotransferase protein.

Efficient and precise measurement of H(alpha)-C(alpha), C(alpha)-C', C(alpha)-C(beta) and H(N)-N residual dipolar couplings from 2D H(N)-N correlation spectra

A suite of experiments are presented for the measurement of H(alpha)-C(alpha), C(alpha)-C', C(alpha)-C(beta) and H(N)-N couplings from uniformly 15N, 13C labeled proteins. Couplings are obtained from a series of intensity modulated two-dimensional H(N)-N spectra equivalent to the common 1H-15N-HSQC spectra, alleviating many overlap and assignment issues associated with other techniques. To illustrate the efficiency of this method, H(alpha)-C(alpha), C(alpha)-C', and H(N)-N isotropic scalar couplings were determined for ubiquitin from data collected in less than 4.5 h, C(alpha)-C(beta) data collection required 10 h. The resulting couplings were measured with an average error of +/-0.06, +/-0.05, +/-0.04 and +/-0.10 Hz, respectively. This study also shows H(alpha)-C(alpha) and C(alpha)-C(beta) couplings, valuable because they provide orientation of bond vectors outside the peptide plane, can be measured in a uniform and precise way. Superior accuracy and precision to existing 3D measurements for C(alpha)-C' couplings and increased precision compared to IPAP measurements for H(N)-N couplings are demonstrated. Minor modifications allow for acquisition of modulated H(N)-C' 2D spectra, which can yield additional well resolved peaks and significantly increase the number of measured RDCs for proteins with crowded 1H-15N resonances.

Structure Determination of Protein⧸RNA Complexes by NMR

Structure determination of protein?RNA complexes in solution provides unique insights into factors that are involved in protein/RNA recognition. Here, we review the methodology used in our laboratory to overcome the challenges of protein?RNA structure determination by nuclear magnetic resonance (NMR). We use as two examples complexes recently solved in our laboratory, the nucleolin RBD12/b2NRE and Rnt1p dsRBD/snR47h complexes. Topics covered are protein and RNA preparation, complex formation, identification of the protein/RNA interface, protein and RNA resonance assignment, intermolecular NOE assignment, and structure calculation and analysis.

Measurements of side-chain 13C-13C residual dipolar couplings in uniformly deuterated proteins

13C-only spectroscopy was used to measure multiple residual (13)C-(13)C dipolar couplings (RDCs) in uniformly deuterated and (13)C-labeled proteins. We demonstrate that (13)C-start and (13)C-observe spectra can be routinely used to measure an extensive set of the side-chain residual (13)C-(13)C dipolar couplings upon partial alignment of human ubiquitin in the presence of bacteriophages Pf1. We establish that, among different broadband polarization transfer schemes, the FLOPSY family can be used to exchange magnetization between a J coupled network of spins while largely decoupling dipolar interactions between these spins. An excellent correlation between measured RDCs and the 3D structure of the protein was observed, indicating a potential use of the (13)C-(13)C RDCs in the structure determination of perdeuterated proteins.

NMR analysis of methyl groups at 100-500 kDa: model systems and Arp2/3 complex

Large macromolecular machines are among the most important and challenging targets for structural and mechanistic analyses. Consequently, there is great interest in development of NMR methods for the study of multicomponent systems in the 50-500 kDa range. Biochemical methods also must be developed in concert to produce such systems in selectively labeled form. Here, we present (1)H/(13)C-HSQC spectra of protonated methyl groups in a model system that mimics molecular weights up to approximately 560 kDa. Signals from side chain methyl groups of Ile, Leu, and Val residues are clearly detectable at correlation times up to approximately 330 ns. We have also developed a biochemical procedure to produce the 240 kDa, heteroheptameric Arp2/3 actin nucleation complex selectively labeled at one subunit and obtained (1)H/(13)C-HSQC spectra of this assembly. Sensitivity in spectra of both the Arp2/3 complex and the model system indicate that methyl groups will be useful sources of information in nonsymmetric systems with molecular weights greater than 600 kDa at concentrations less than 100 microM. Methyl analyses will complement TROSY and CRINEPT analyses of amides in NMR studies of structure and molecular interactions of extremely large macromolecules and assemblies.

Partial alignment of biomolecules: an aid to NMR characterization

Partial alignment of biomolecules in solution has added a new dimension to structural investigation by high-resolution NMR methods. Applications to proteins, nucleic acids and carbohydrates now abound. Limitations initially associated with compatibility of biomolecules with the liquid-crystal media commonly used to achieve alignment have begun to disappear. This is, in part, a result of the introduction of a wide variety of new media. Future applications to biologically important problems such as the structural organization of multi-domain proteins and multi-protein assemblies look very promising.

Ligand-induced structural changes to maltodextrin-binding protein as studied by solution NMR spectroscopy

Solution NMR studies on the physiologically relevant ligand-free and maltotriose-bound states of maltodextrin-binding protein (MBP) are presented. Together with existing data on MBP in complex with beta-cyclodextrin (non-physiological, inactive ligand), these new results provide valuable information on changes in local structure, dynamics and global fold that occur upon ligand binding to this two-domain protein. By measuring a large number of different one-bond residual dipolar couplings, the domain conformations, critical for biological function, were investigated for all three states of MBP. Structural models of the solution conformation of MBP in a number of different forms were generated from the experimental dipolar coupling data and X-ray crystal structures using a quasi-rigid-body domain orientation algorithm implemented in the structure calculation program CNS. Excellent agreement between relative domain orientations in ligand-free and maltotriose-bound solution conformations and the corresponding crystal structures is observed. These results are in contrast to those obtained for the MBP/beta-cyclodextrin complex where the solution state is found to be approximately 10 degrees more closed than the crystalline state. The present study highlights the utility of residual dipolar couplings for orienting protein domains or macromolecules with respect to each other.