Determination of the packing mode of the coiled-coil domain of cGMP-dependent protein kinase Ialpha in solution using charge-predicted dipolar couplings

Coiled-coil motifs are ubiquitous in biology and play essential roles in protein assembly and molecular recognition. Here, we show that the relative orientation and stoichiometry of coiled-coil proteins in solution can be determined by comparison of residual dipolar couplings (RDCs) measured in charged liquid-crystalline medium with values predicted from the three-dimensional charge distribution of the protein. Comparison of charge-predicted RDCs with a small set of one-bond 1DNH dipolar couplings, measured in the negatively charged liquid-crystalline Pf1 bacteriophage medium, identified the coiled-coil region of the cGMP-dependent protein kinase I as a parallel homodimer in solution and ruled out an antiparallel dimeric or monomeric state. The method is very rapid, applicable to a wide variety of liquid crystals used in biological NMR to date, and can be applied to coiled-coil structures and other proteins with higher order assembly.     

Protein alignment by a coexpressed lanthanide-binding tag for the measurement of residual dipolar couplings

A protein fusion construct of human ubiquitin with an N-terminal lanthanide binding tag (LBT) enables observation of long-range orientational restraints in solution NMR from residual dipolar couplings (RDCs) due to paramagnetic alignment of the protein. The paramagnetic lanthanide ions Tb3+, Dy3+, and Tm3+ are shown to bind to the LBT and induce different alignment tensors, in agreement with theory. RDCs, measured relative to the diamagnetic Lu3+, range from -7.6 to 5.5 Hz for Tb3+ and -6.6 to 6.1 Hz for Dy3+, while an opposite alignment tensor is observed for Tm3+ (4.5 to -2.9 Hz) at 800 MHz. Experimental RDCs are in excellent agreement with those predicted on the basis of the X-ray structure of the protein.     

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.     

16-fold degeneracy of peptide plane orientations from residual dipolar couplings: analytical treatment and implications for protein structure determination

Residual dipolar couplings (RDCs) measured for internally rigid molecular fragments provide important information about the relative orientations of these fragments. Dependent on the symmetry of the alignment tensor and the symmetry of the molecular fragment, however, there generally exist more than one solution for the fragment orientation consistent with the measured RDCs. Analytical solutions are presented that describe the complete set of orientations of internally rigid fragments that are consistent with multiple dipolar couplings measured in a single alignment medium that is rhombic. For the first time, it is shown that, for a planar fragment such as the peptide plane, there generally exist 16 different solutions with their analytical expressions presented explicitly. The presence of these solutions is shown to be highly relevant for standard structure determination protocols using RDCs to refine molecular structures. In particular, when using standard protein structure refinement with RDCs that were measured in a single alignment medium as constraints, it is found that often more than one of the peptide plane solutions is physically viable; i.e., despite being consistent with measured RDCs, the local backbone structure can be incorrect. On the basis of experimental and simulated examples, it is rationalized why protein structures that are refined against RDCs measured in a single medium can have lower resolution (precision) than one would expect on the basis of the experimental accuracy of the RDCs. Conditions are discussed under which the correct solution can be identified.     

Giant Residual Dipolar 13C–1H Couplings in High‐Spin Organoiron Complexes: Elucidation of Their Structures in Solution by 13C NMR Spectroscopy

High‐spin Fe^II–alkyl complexes with bis(pyridylimino)isoindolato ligands were synthesized and their paramagnetic ¹H and ¹³C NMR spectra were analyzed comprehensively. The experimental ¹³C—¹H coupling values are temperature (T^?1)‐ as well as magnetic‐field (B²)‐dependent and deviate considerably from typical scalar ¹J_CH couplings constants. This deviation is attributed to residual dipolar couplings (RDCs), which arise from partial alignment of the complexes in the presence of a strong magnetic field. The analysis of the experimental RDCs allows an unambiguous assignment of all ¹³C NMR resonances and, additionally, a structural refinement of the conformation of the complexes in solution. Moreover the RDCs can be used for the analysis of the alignment tensor and hence the tensor of the anisotropy of the magnetic susceptibility.     

Simultaneous assignment of all diastereotopic protons in strychnine using RDCs: PELG as alignment medium for organic molecules

The concept of using residual dipolar couplings (RDCs) for the structure determination of organic molecules is applied to the simultaneous assignment of all diastereotopic protons in strychnine. To use this important NMR parameter the molecule has to be aligned in the magnetic field. Here we present a new alignment medium for organic substrates. The optimization of the alignment properties of mixtures of poly-gamma-ethyl-L-glutamate (PELG) and CDCl(3) are described and the alignment properties of PELG at different concentrations are evaluated. A comparison of PELG with poly-gamma-benzyl-L-glutamate (PBLG) shows considerable differences in the magnitude of alignment for strychnine in the two alignment media. PELG induces a lower degree of order and makes the measurement of residual dipolar couplings (RDCs) in strychnine possible. All one-bond C-H RDCs of strychnine in PELG were determined by using 2D heteronuclear single quantum coherence (HSQC) spectroscopy. The strategy for the extraction of RDCs for methylene groups is described in detail. The RDCs and order parameters are used to assign pairs of diastereotopic protons. This methodology can distinguish not only one pair of diastereotopic protons but it can be used to assign all pairs of diastereotopic protons simultaneously. Two different calculation approaches to achieve this task are described in detail.     

A caged lanthanide complex as a paramagnetic shift agent for protein NMR

A lanthanide complex, named CLaNP (caged lanthanide NMR probe) has been developed for the characterisation of proteins by paramagnetic NMR spectroscopy. The probe consists of a lanthanide chelated by a derivative of DTPA (diethylenetriaminepentaacetic acid) with two thiol reactive functional groups. The CLaNP molecule is attached to a protein by two engineered, surface-exposed, Cys residues in a bidentate manner. This drastically limits the dynamics of the metal relative to the protein and enables measurements of pseudocontact shifts. NMR spectroscopy experiments on a diamagnetic control and the crystal structure of the probe-protein complex demonstrate that the protein structure is not affected by probe attachment. The probe is able to induce pseudocontact shifts to at least 40 A from the metal and causes residual dipolar couplings due to alignment at a high magnetic field. The molecule exists in several isomeric forms with different paramagnetic tensors; this provides a fast way to obtain long-range distance restraints.     

Periodicity in residual dipolar couplings and nucleic acid structures

The periodicity in nucleic acid duplex structures is shown to be correlated to the periodicity in residual dipolar couplings (RDCs) in the form of an "RDC wave". This "RDC wave" is characteristic of the alignment of the duplex in the magnetic field, and hence fitting of the data allows the duplex global orientation (, Phi) to be extracted. Further, because the "RDC wave" is fit as a data set of a corresponding secondary structure element, the degeneracy problem is greatly reduced. Consequently, with the global orientation (, Phi) determined, local bond vector conformations are defined. The fit is demonstrated in the examples of the imino RDCs of the negative regulator of splicing RNA fragment (NRS23) and for the C1'H1' RDCs of the Dickerson dodecamer.     

Rigid,Highly Reactive and Stable DOTA-like Tags Containing a Thiol-Specific Phenylsulfonyl Pyridine Moiety for Protein Modification and NMR Analysis**

Site specific installation of a paramagnetic ion with magnetic anisotropy in a biomolecule generates valuable structural restraints, such as pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs). These paramagnetic effects can be used to characterize the structures, interactions and dynamics of biological macromolecules and their complexes. Two single-armed DOTA-like tags, BrPSPy-DO3M(S)A-Ln and BrPSPy-6M-DO3M(S)A-Ln, each containing a thiol-specific reacting group, that is, a phenylsulfonyl pyridine moiety, are demonstrated as rigid, reactive and stable paramagnetic tags for protein modification by formation of a reducing resistant thioether bond between the protein and the tag. The two tags present high reactivity with the solvent exposed thiol group in aqueous solution at room temperature. The introduction of Br at the meta-position in pyridine enhances the reactivity of 4-phenylsulfonyl pyridine towards the solvent exposed thiol group in a protein, whereas the ortho-methyl group in pyridine increases the rigidity of the tag in the protein conjugates. The high performance of these two tags has been demonstrated in different cysteine mutants of ubiquitin and GB1. The high reactivity and rigidity of these two tags can be added in the toolbox of paramagnetic tags suitable for the high-resolution NMR measurements of biological macromolecules and their complexes.     

Engineering encodable lanthanide-binding tags into loop regions of proteins

Lanthanide-binding tags (LBTs) are valuable tools for investigation of protein structure, function, and dynamics by NMR spectroscopy, X-ray crystallography, and luminescence studies. We have inserted LBTs into three different loop positions (denoted L, R, and S) of the model protein interleukin-1β (IL1β) and varied the length of the spacer between the LBT and the protein (denoted 1?3). Luminescence studies demonstrate that all nine constructs bind Tb3+ tightly in the low nanomolar range. No significant change in the fusion protein occurs from insertion of the LBT, as shown by two X-ray crystallographic structures of the IL1β-S1 and IL1β-L3 constructs and for the remaining constructs by comparing the 1H?15N heteronuclear single-quantum coherence NMR spectra with that of the wild-type IL1β. Additionally, binding of LBT-loop IL1β proteins to their native binding partner in vitro remains unaltered. X-ray crystallographic phasing was successful using only the signal from the bound lanthanide. Large residual dipolar couplings (RDCs) could be determined by NMR spectroscopy for all LBT-loop constructs and revealed that the LBT-2 series were rigidly incorporated into the interleukin-1β structure. The paramagnetic NMR spectra of loop-LBT mutant IL1β-R2 were assigned and the Δχ tensor components were calculated on the basis of RDCs and pseudocontact shifts. A structural model of the IL1β-R2 construct was calculated using the paramagnetic restraints. The current data provide support that encodable LBTs serve as versatile biophysical tags when inserted into loop regions of proteins of known structure or predicted via homology modeling.     

Conformational Ensemble of Disordered Proteins Probed by Solvent Paramagnetic Relaxation Enhancement (sPRE)

Characterization of the conformational ensemble of disordered proteins is highly important for understanding protein folding and aggregation mechanisms, but remains a computational and experimental challenge owing to the dynamic nature of these proteins. New observables that can provide unique insights into transient residual structures in disordered proteins are needed. Here using denatured ubiquitin as a model system, NMR solvent paramagnetic relaxation enhancement (sPRE) measurements provide an accurate and highly sensitive probe for detecting low populations of residual structure in a disordered protein. Furthermore, a new ensemble calculation approach based on sPRE restraints in conjunction with residual dipolar couplings (RDCs) and small‐angle X‐ray scattering (SAXS) is used to define the conformational ensemble of disordered proteins at atomic resolution. The approach presented should be applicable to a wide range of dynamic macromolecules.     

CHARMM36 all‐atom additive protein force field: Validation based on comparison to NMR data

Protein structure and dynamics can be characterized on the atomistic level with both nuclear magnetic resonance (NMR) experiments and molecular dynamics (MD) simulations. Here, we quantify the ability of the recently presented CHARMM36 (C36) force field (FF) to reproduce various NMR observables using MD simulations. The studied NMR properties include backbone scalar couplings across hydrogen bonds, residual dipolar couplings (RDCs) and relaxation order parameter, as well as scalar couplings, RDCs, and order parameters for side‐chain amino‐ and methyl‐containing groups. It is shown that the C36 FF leads to better correlation with experimental data compared to the CHARMM22/CMAP FF and suggest using C36 in protein simulations. Although both CHARMM FFs contains the same nonbond parameters, our results show how the changes in the internal parameters associated with the peptide backbone via CMAP and the χ₁ and χ₂ dihedral parameters leads to improved treatment of the analyzed nonbond interactions. This highlights the importance of proper treatment of the internal covalent components in modeling nonbond interactions with molecular mechanics FFs. © 2013 Wiley Periodicals, Inc.     

Dipolar waves as NMR maps of protein structure

The anisotropy of nuclear spin interactions results in a unique mapping of structure to the resonance frequencies and split tings observed in NMR spectra, however, the determination of molecular structure from experimentally measured spectral parameters is complicated by angular ambiguities resulting from the symmetry properties of dipole-dipole and chemical shift interactions. This issue can be addressed through the periodicity inherent in secondary structure elements, which can be used as an index of topology. Distinctive wheel-like patterns are observed in two-dimensional 1H-15N heteronuclear dipolar/15N chemical shift PISEMA (polarization inversion spin-exchange at the magic angle) spectra of helical membrane proteins in highly aligned lipid bilayer samples. One-dimensional dipolar waves are an extension of two-dimensional PISA (polarity index slant angle) wheels to map protein structure in NMR spectra of both highly and weakly aligned samples. Dipolar waves describe the periodic wavelike variations of the magnitudes of the static heteronuclear dipolar couplings as a function of residue number in the absence of chemical shift effects. Weakly aligned samples of proteins display these same effects, primarily as residual dipolar couplings (RDCs), in solution NMR spectra. The corresponding properties of the RDCs in solution NMR spectra of weakly aligned helices represent a convergence of solid-state and solution NMR approaches to structure determination.     

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.     

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.     

Structure determination of protein-ligand complexes by transferred paramagnetic shifts

Rational drug design depends on the knowledge of the three-dimensional (3D) structure of complexes between proteins and lead compounds of low molecular weight. A novel nuclear magnetic resonance (NMR) spectroscopy strategy based on the paramagnetic effects from lanthanide ions allows the rapid determination of the 3D structure of a small ligand molecule bound to its protein target in solution and, simultaneously, its location and orientation with respect to the protein. The method relies on the presence of a lanthanide ion in the protein target and on fast exchange between bound and free ligand. The binding affinity of the ligand and the paramagnetic effects experienced in the bound state are derived from concentration-dependent (1)H and (13)C spectra of the ligand at natural isotopic abundance. Combined with prior knowledge of the crystal or solution structure of the protein and of the magnetic susceptibility tensor of the lanthanide ion, the paramagnetic data define the location and orientation of the bound ligand molecule with respect to the protein from simple 1D NMR spectra. The method was verified with the ternary 30 kDa complex between the lanthanide-labeled N-terminal domain of the epsilon exonuclease subunit from the Escherichia coli DNA polymerase III, the subunit theta, and thymidine. The binding mode of thymidine was found to be very similar to that of thymidine monophosphate present in the crystal structure.     

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.     

Calculation of residual dipolar couplings from disordered state ensembles using local alignment

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.