Conformational analysis of a flexible oligosaccharide using residual dipolar couplings.

We present a new approach to the analysis of the conformational and the motional properties of an oligosaccharide, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside. The approach relies on an order matrix analysis of residual dipolar couplings in the solution state. By combining a number of different types of couplings, (1)D(CH), (2)D(CH), and D(HH), an order matrix is solved for each ring of the trimannoside. The resulting order parameters indicate the internal motion at the alpha (1,3) linkage to be limited, while significant motion is suggested at the alpha (1,6) linkage. Two structures for the trimannoside were determined by aligning the order tensor principal axes obtained from two different orienting media, bicelles and phage. The very similar conformations at the alpha (1,3) linkage of these two structures confirm that the internal motion at the alpha (1,3) linkage is small and the conformation is a good representation of a single preferred structure. The different conformations at the alpha (1,6) linkage suggest that the motional amplitudes are large and the conformations must be viewed as virtual conformers. Compared with traditional NMR methods, data acquisition is easy and data analysis is straightforward.     

On the interpretation of residual dipolar couplings as reporters of molecular dynamics

The analysis of residual dipolar couplings from an ensemble of conformations to extract molecular dynamics is intricate. The very mechanism that is necessary to perturb overall molecular tumbling to generate nonvanishing residual dipolar couplings gives rise to convoluted data. The measured values are essentially weighted averages over conformations. However, the weights are not simply the populations of conformations. Consequently, the observed order parameter is not exactly the true measure of motion. In the case of paramagnetic alignment, the apparent order parameter is expected to depend on the number of torsions that separate the locus of interest from the paramagnetic site. In the case of alignment due to steric obstruction, the uneven selection of conformations by their differing Saupe order matrices leads to a bias in the residual dipolar couplings-probed molecular dynamics.     

Simultaneous definition of high resolution protein structure and backbone conformational dynamics using NMR residual dipolar couplings.

Despite the importance of molecular dynamics for biological activity, most approaches to protein structure determination, whether based on crystallographic or solution studies, propose three-dimensional atomic representations of a single configuration that take no account of conformational fluctuation. Non-averaged anisotropic NMR interactions, such as residual dipolar couplings, that become measurable under conditions of weak alignment, provide sensitive probes of both molecular structure and dynamics. Residual dipolar couplings are becoming increasingly powerful for the study of proteins in solution. In this minireview we present their use for the simultaneous determination of protein structure and dynamics.     

Direct prediction of residual dipolar couplings of small molecules in a stretched gel by stochastic molecular dynamics simulations

Residual dipolar couplings are highly useful NMR parameters for calculating and refining molecular structures, dynamics, and interactions. For some applications, however, it is inevitable that the preferred orientation of a molecule in an alignment medium is calculated a priori. Several methods have been developed to predict molecular orientations and residual dipolar couplings. Being beneficial for macromolecules and selected small‐molecule applications, such approaches lack sufficient accuracy for a large number of organic compounds for which the fine structure and eventually the flexibility of all involved molecules have to be considered or are limited to specific, well‐studied liquid crystals. We introduce a simplified model for detailed all‐atom molecular dynamics calculations with a polymer strand lined up along the principal axis as a new approach to simulate the preferred orientation of small to medium‐sized solutes in polymer‐based, gel‐type alignment media. As is shown by a first example of strychnine in a polystyrene/CDCl₃ gel, the simulations potentially enable the accurate prediction of residual dipolar couplings taking into account structural details and dynamic averaging effects of both the polymer and the solute. Copyright © 2015 John Wiley & Sons, Ltd.     

Limited flexibility of lactose detected from residual dipolar couplings using molecular dynamics simulations and steric alignment methods

The conformational flexibility of lactose in solution has been investigated by residual dipolar couplings (RDCs). One-bond carbon-proton and proton-proton coupling constants have been measured in two oriented media and interpreted in combination with molecular dynamics simulations (MD). Two different approaches, known as PALES (Zweckstetter et al., J. Am. Chem. Soc. 2000, 122, 3791-3792) and TRAMITE (Azurmendi et al., J. Am. Chem. Soc. 2002, 124, 2426-2427), have been used to determine the alignment tensor from a shape-induced alignment model with the oriented medium. The steric alignment of the structures from several MD trajectories has provided ensemble averaged RDCs that have been compared with the experimental ones. The obtained results reveal the almost exclusive presence of a major low energy region defined as syn-phi/syn-psi (> 97%), for which sampling occurs in a dynamic manner. This result satisfactorily agrees with that determined by standard NOE-based methods.     

Lanthanide-binding peptides for NMR measurements of residual dipolar couplings and paramagnetic effects from multiple angles

Lanthanide-binding peptide tags (LBTs) containing a single cysteine residue can be attached to proteins via a disulfide bond, presenting a flexible means of tagging proteins site-specifically with a lanthanide ion. Here we show that cysteine residues placed in different positions of the LBT can be used to expose the protein to different orientations of the magnetic susceptibility anisotropy (delta chi) tensor and to generate different molecular alignments in a magnetic field. Delta chi tensors determined by nuclear magnetic resonance (NMR) spectroscopy for LBT complexes with Yb3+, Tm3+, and Er3+ suggest a rational way of producing alignment tensors with different orientations. In addition, knowledge of the delta chi tensor of LBT allows modeling of the protein-LBT structures. Despite evidence for residual mobility of the LBTs with respect to the protein, the pseudocontact shifts and residual dipolar couplings displayed by proteins disulfide-bonded to LBTs are greater than those achievable with most other lanthanide binding tags.     

Determination of the structures of symmetric protein oligomers from NMR chemical shifts and residual dipolar couplings

Symmetric protein dimers, trimers, and higher-order cyclic oligomers play key roles in many biological processes. However, structural studies of oligomeric systems by solution NMR can be difficult due to slow tumbling of the system and the difficulty in identifying NOE interactions across protein interfaces. Here, we present an automated method (RosettaOligomers) for determining the solution structures of oligomeric systems using only chemical shifts, sparse NOEs, and domain orientation restraints from residual dipolar couplings (RDCs) without a need for a previously determined structure of the monomeric subunit. The method integrates previously developed Rosetta protocols for solving the structures of monomeric proteins using sparse NMR data and for predicting the structures of both nonintertwined and intertwined symmetric oligomers. We illustrated the performance of the method using a benchmark set of nine protein dimers, one trimer, and one tetramer with available experimental data and various interface topologies. The final converged structures are found to be in good agreement with both experimental data and previously published high-resolution structures. The new approach is more readily applicable to large oligomeric systems than conventional structure-determination protocols, which often require a large number of NOEs, and will likely become increasingly relevant as more high-molecular weight systems are studied by NMR.     

Highly populated turn conformations in natively unfolded tau protein identified from residual dipolar couplings and molecular simulation

Tau, a natively unstructured protein that regulates the organization of neuronal microtubules, is also found in high concentrations in neurofibrillary tangles of Alzheimer's disease and other neurodegenerative disorders. The conformational transition between these vastly different healthy and pathological forms remains poorly understood. We have measured residual dipolar couplings (RDCs), J-couplings, and nuclear Overhauser enhancement (NOE) in construct K18 of tau, containing all four repeat domains R1-R4. NHN RDCs were compared with prediction on the basis of a statistical model describing the intrinsic conformational sampling of unfolded proteins in solution. While local variation and relative amplitude of RDCs agrees with propensity-based prediction for most of the protein, homologous sequences in each repeat domain (DLKN, DLSN, DLSK, and DKFD in repeats R1-R4) show strong disagreement characterized by inversion of the sign of the central couplings. Accelerated molecular dynamic simulations (AMD) in explicit solvent revealed strong tendencies to form turns, identified as type I beta-turns for repeats R1-R3. Incorporation of the backbone dihedral sampling resulting from AMD into the statistical coil model closely reproduces experimental RDC values. These localized sequence-dependent conformational tendencies interrupt the propensity to sample more extended conformations in adjacent strands and are remarkably resistant to local environmental factors, as demonstrated by the persistence of the RDC signature even under harsh denaturing conditions (8 M urea). The role that this specific conformational behavior may play in the transition to the pathological form is discussed.     

Refinement of local and long-range structural order in theophylline-binding RNA using (13)C-(1)H residual dipolar couplings and restrained molecular dynamics.

13C-(1)H residual dipolar couplings (RDC) have been measured for the bases and sugars in the theophylline-binding RNA aptamer, dissolved in filamentous phage medium, and used to investigate the long-range structural and dynamic behavior of the molecule in the solution state. The orientation dependent RDC provide additional restraints to further refine the overall structure of the RNA-theophylline complex, whose long-range order was poorly defined in the NOE-based structural ensemble. Structure refinement using RDC normally assumes that molecular alignment can be characterized by a single tensor and that the molecule is essentially rigid. To address the validity of this assumption for the complex of interest, we have analyzed distinct domains of the RNA molecule separately, so that local structure and alignment tensors experienced by each region are independently determined. Alignment tensors for the stem regions of the molecule were allowed to float freely during a restrained molecular dynamics structure refinement protocol and found to converge to similar magnitudes. During the second stage of the calculation, a single alignment tensor was thus applied for the whole molecule and an average molecular conformation satisfying all experimental data was determined. Semirigid-body molecular dynamics calculations were used to reorient the refined helical regions to a relative orientation consistent with this alignment tensor, allowing determination of the global conformation of the molecule. Simultaneously, the local structure of the theophylline-binding core of the molecule was refined under the influence of this common tensor. The final ensemble has an average pairwise root mean square deviation of 1.50 +/- 0.19 A taken over all heavy atoms, compared to 3.5 +/- 1.1 A for the ensemble determined without residual dipolar coupling. This study illustrates the importance of considering both the local and long-range nature of RDC when applying these restraints to structure refinements of nucleic acids.     

Model-free approach to the dynamic interpretation of residual dipolar couplings in globular proteins.

The effects of internal motions on residual dipolar NMR couplings of proteins partially aligned in a liquid-crystalline environment are analyzed using a 10 ns molecular dynamics (MD) computer simulation of ubiquitin. For a set of alignment tensors with different orientations and rhombicities, MD-averaged dipolar couplings are determined and subsequently interpreted for different scenarios in terms of effective alignment tensors, average orientations of dipolar vectors, and intramolecular reorientational vector distributions. Analytical relationships are derived that reflect similarities and differences between motional scaling of dipolar couplings and scaling of dipolar relaxation data (NMR order parameters). Application of the self-consistent procedure presented here to dipolar coupling measurements of biomolecules aligned in different liquid-crystalline media should allow one to extract in a "model-free" way average orientations of dipolar vectors and specific aspects of their motions.     

Measurement of five dipolar couplings from a single 3D NMR multiplet applied to the study of RNA dynamics

A new three-dimensional NMR experiment is described that yields five scalar or dipolar couplings from a single cross-peak between three spins. The method is based on the E.COSY principle and is demonstrated for the H1'-C1'-C2' fragment of ribose sugars in a uniformly 13C-enriched 24-nucleotide RNA stem-loop structure, for which a complete set of couplings was obtained for all nonmodified nucleotides. The values of the isotropic J couplings and the 13C1' and 13C2' chemical shifts define the sugar pucker. Once the sugar pucker is known, the five dipolar couplings between C1'-H1', C2'-H2', H1'-H2', C1'-H2', and C2'-H1', together with C1'-C2', C3'-H3', and C4'-H4' available from standard experiments, can be used to derive the five unknowns that define the local alignment tensor, thereby simultaneously providing information on relative orientation and dynamics of the ribose units. Data indicate rather uniform alignment for all stem nucleotides in the 24-nt stem-loop structure, with only a modest reduction in order for the terminal basepair, but significantly increased mobility in part of the loop region. The method is applicable to proteins, nucleic acids, and carbohydrates, provided 13C enrichment is available.     

Weak alignment of biomacromolecules in collagen gels: an alternative way to yield residual dipolar couplings for NMR measurements

Collagen, consisting of glycine, proline, and hydroxyproline, is a fibrous protein that can form a rope-like left-hand triple helix structure. It is demonstrated here that the collagen gels prepared from polymerization in the magnetic field can provide weak alignment for protein. The alignment order induced by collagen gels is quite small when compared to other alignment media, but the magnitude of the dipolar couplings can be easily scaled up by increasing the initial concentration of collagen. The collagen gels showed good pH and detergent tolerance. These advantages of collagen gels make it a promising candidate for the alignment of large biomolecules or membrane protein-detergent complexes in the magnetic field.     

The influence of molecular shape and charge distribution on the molecular approach in solution. NMR dipolar relaxation studies of the many-body effects

The paramagnetic dipolar spin-lattice relaxation of nuclear spins is expressed in terms of the geometrical and dynamical parameters of the molecules and used to analyze the intra- and intermolecular microdynamics. The possible conformational changes of a TEMPOL nitroxide-labelled sugar in deuterated chloroform are studied and then the intermolecular processes are considered. The influence of the discrete polar natures of acetonitrile, chloroform, and methanol on the TEMPOL solvation dynamics is recognized. The importance of the water structure on the collisional approach of two attractive ions is emphasized. Finally, the attractive effects of the dipole-dipole electrostatic forces are shown for the water/TEMPOL pair. Diffusion and jump analytical theories, integral equations of statistical mechanics, and Monte-Carlo simulations of the diffusion are combined to treat the anisotropic interactions and the many-body effects.     

An alternative method for pucker determination in carbohydrates from residual dipolar couplings: a solution NMR study of the fructofuranosyl ring of sucrose

A simple solution NMR method is presented for pucker determination of five-membered rings using only residual dipolar couplings obtained in a single liquid crystalline medium, DMPC/DHPC bicelles (DMPC = dimyristoylphosphatidylcholine; DHPC = dihexanoylphosphatidylcholine). The method was applied to determine the pucker of the fructofuranosyl ring of sucrose. The results indicate a fructofuranosyl pucker phase in the 20 degrees - 70 degrees range. The pucker phases are in agreement with those from previous NMR and optical spectroscopic studies and, importantly, do not rely on empirically parametrized Karplus curves. Furthermore, the analysis implies more than one stable pucker phase and rapid ring interconversion in this range. The present results suggest that using residual dipolar couplings alone can reveal multiple conformations present in solution. Hence, when a sufficient number of residual dipolar coupling constants is measured, the outcome is a robust, reliable, and independent route for determining carbohydrate and nucleic acid structure by NMR spectroscopy.     

Intrinsic information content of NMR dipolar couplings: a conformational investigation of 1,3-butadiene in a nematic phase.

The conformational equilibrium of 1,3-butadiene in a condensed fluid phase is investigated by liquid-crystal NMR spectroscopy. The full set of D(HH) and D(CH) dipolar couplings is determined from the analysis of the (1)H spectra of the three 1,3-butadiene most-abundant isotopomers (i.e. the all (12)C and the two single-labeled (13)C isotopomers) for a total of 21 independent dipolar couplings. A very good starting set of spectral parameters for the analysis of the (1)H spectrum is determined in a semiautomated way by the analysis of the (N-1) (specifically, N=6, the number of 1/2 spin nuclei in the spin system) quantum refocused (5QR), and not (5Q), spectra. As an alternative approach, a Monte Carlo (MC) numerical simulation, capable of predicting the solute ordering, is tested to simulate the 5QR spectrum. The set of D(ij) couplings is very good, proving that the MC method can represent a novel, valid alternative to the existing spectral simplification procedures. The experimentally determined dipolar-coupling data set is fully compatible with the 1,3-butadiene conformational distribution reported in the literature for isolated molecules, indicating the presence of about 99 % of s-trans conformer. With regards to the remaining 1 %, in spite of the direct and very strong dependence of the observables on the molecular structure, it was not possible to discriminate between the planar s-cis and s-gauche forms, both of which produce a very good fit of the dipolar couplings. Vibrational corrections, up to the anharmonic term, were applied; the calculated geometrical parameters are in good- although not exact-agreement with those reported in the literature from experimental and theoretical investigations. This result can be considered as supporting the methodology used for obtaining the structure and conformational distribution of a flexible molecule in a liquid phase.