Sets Of Residual Dipolar Couplings Research Articles

Programming Peptide Liquid Crystal Media to Acquire Independent Sets of Residual Dipolar Couplings and Enantiodiscrimination in Multiple Solvent Systems.

Multiple independent sets of residual dipolar couplings (RDCs) acquired by relying on different alignment media show the great potential for de novo structure determination of organic compounds. However, this methodology is severely compromised by the limited availability of multialignment media. In this work, an engineering strategy was developed to program the oligopeptide amphiphiles (OPAs) to create different peptide liquid crystal (LC) media for the acquisition of independent sets of RDCs. With no need for de novo design on peptide sequences, the molecular alignment can be simply modulated by varying the length of the hydrophobic tails within OPAs. Relying on these programmed peptide LC media, five independent sets of RDCs were extracted in a highly efficient and accurate manner. Because of the similar bulk composition of OPAs, this approach offers the significant advantage in circumventing the possible incompatibilities of analytes with one or several different alignment media, therefore avoiding the analysis complication. Notably, these peptide LC media show enantiodifferentiating properties, and the enantiodiscriminating capabilities could also be optimized through the programmed strategy. Furthermore, we show that these media are compatible with different polar solvents, allowing the possible de novo structure elucidation of organic compounds with varied polarities and solubilities.

Programmable alignment media from self-assembled oligopeptide amphiphiles for the measurement of independent sets of residual dipolar couplings in organic solvents.

NMR spectroscopy in anisotropic media has emerged as a powerful technique for the structural elucidation of organic molecules. Its application requires weak alignment of analytes by means of suitable alignment media. Although a number of alignment media, that are compatible with organic solvents, have been introduced in the last 20 years, acquiring a number of independent, non-linearly related sets of anisotropic NMR data from the same organic solvent system remains a formidable challenge, which is however crucial for the alignment simulations and deriving dynamic and structural information of organic molecules unambiguously. Herein, we introduce a programmable strategy to construct several distinct peptide-based alignment media by adjusting the amino acid sequence, which allows us to measure independent sets of residual dipolar couplings (RDCs) in a highly efficient and accurate manner. This study opens a new avenue for de novo structure determination of organic compounds without requiring prior structural information.

Poly(arylisocyanides) as Versatile, Enantiodiscriminating Alignment Media for Small Molecules.

Lyotropic liquid crystalline (LLC) phases of amino acid derived polyarylisocyanides were employed as chiral alignment media for the measurement of residual dipolar couplings (RDCs) of small chiral organic molecules. Anisotropic samples in CDCl3 displayed quadrupolar splittings of the deuterium signal in the range of several hundreds of Hertz. The LLC phases showed excellent orienting properties for a broad range of analytes bearing various functional groups. The precise extraction of RDCs in the range of up to ±40 Hertz from F2-coupled HSQC spectra was possible. Additionally, the chiral environment offers the opportunity for diastereomorphous interactions with the enantiomers of chiral analytes leading to two different sets of RDCs. This differential order effect was particularly pronounced with ketones and alcohols.

Valine derived poly (acetylenes) as versatile chiral lyotropic liquid crystalline alignment media for RDC-based structure elucidations.

Anisotropic samples of lyotropic liquid crystalline (LLC) phases of valine derived polyaryl acetylenes were employed as chiral alignment media for the measurement of residual dipolar couplings (RDCs) of 12 small, chiral, organic molecules. The quadrupolar splitting of the deuterium signal of CDCl3 can be adjusted by temperature and concentration changes from 0 to 350 Hz. The LLC phases showed excellent orienting properties for all analytes bearing various functional groups. The precise extraction of RDCs in the range of up to ±30 Hz from F2-coupled HSQC spectra was possible. Additionally, the chiral environment led to diastereomorphous interactions with the enantiomers of chiral analytes leading to two different sets of RDCs. This differential order effect was particularly pronounced with H-bond donors like alcohols and 2° amines.

Inter-helical conformational preferences of HIV-1 TAR-RNA from maximum occurrence analysis of NMR data and molecular dynamics simulations.

Detecting conformational heterogeneity in biological macromolecules is a key for the understanding of their biological function. We here provide a comparison between two independent approaches to assess conformational heterogeneity: molecular dynamics simulations, performed without inclusion of any experimental data, and maximum occurrence (MaxOcc) distribution over the topologically available conformational space. The latter only reflects the extent of the averaging and identifies regions which are most compliant with the experimentally measured NMR Residual Dipolar Couplings (RDCs). The analysis was performed for the HIV-1 TAR RNA, consisting of two helical domains connected by a flexible bulge junction, for which four sets of RDCs were available as well as an 8.2 μs all-atom molecular dynamics simulation. A sample and select approach was previously applied to extract from the molecular dynamics trajectory conformational ensembles in agreement with the four sets of RDCs. The MaxOcc analysis performed here identifies the most likely sampled region in the conformational space of the system which, strikingly, overlaps well with the structures independently sampled in the molecular dynamics calculations and even better with the RDC selected ensemble.

Side Chain Conformational Distributions of a Small Protein Derived from Model-Free Analysis of a Large Set of Residual Dipolar Couplings.

Accurate quantitative measurement of structural dispersion in proteins remains a prime challenge to both X-ray crystallography and NMR spectroscopy. Here we use a model-free approach based on measurement of many residual dipolar couplings (RDCs) in differentially orienting aqueous liquid crystalline solutions to obtain the side chain χ1 distribution sampled by each residue in solution. Applied to the small well-ordered model protein GB3, our approach reveals that the RDC data are compatible with a single narrow distribution of side chain χ1 angles for only about 40% of the residues. For more than half of the residues, populations greater than 10% for a second rotamer are observed, and four residues require sampling of three rotameric states to fit the RDC data. In virtually all cases, sampled χ1 values are found to center closely around ideal g–, g+ and t rotameric angles, even though no rotamer restraint is used when deriving the sampled angles. The root-mean-square difference between experimental 3JHαHβ couplings and those predicted by the Haasnoot-parametrized, motion-adjusted Karplus equation reduces from 2.05 to 0.75 Hz when using the new rotamer analysis instead of the 1.1-Å X-ray structure as input for the dihedral angles. A comparison between observed and predicted 3JHαHβ values suggests that the root-mean-square amplitude of χ1 angle fluctuations within a given rotamer well is ca. 20°. The quantitatively defined side chain rotamer equilibria obtained from our study set new benchmarks for evaluating improved molecular dynamics force fields, and also will enable further development of quantitative relations between side chain chemical shift and structure.

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings.

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

High accuracy of Karplus equations for relating three-bond J couplings to protein backbone torsion angles.

(3) JC'C' and (3) JHNHα couplings are related to the intervening backbone torsion angle ${\varphi }$ by standard Karplus equations. Although these couplings are known to be affected by parameters other than ${\varphi }$, including H-bonding, valence angles and residue type, experimental results and quantum calculations indicate that the impact of these latter parameters is typically very small. The solution NMR structure of protein GB3, newly refined by using extensive sets of residual dipolar couplings, yields 50-60 % better Karplus equation agreement between ${\varphi }$ angles and experimental (3) JC'C' and (3) JHNHα values than does the high-resolution X-ray structure. In intrinsically disordered proteins, (3) JC'C' and (3) JHNHα couplings can be measured at even higher accuracy, and the impact of factors other than the intervening torsion angle on (3) J will be smaller than in folded proteins, making these couplings exceptionally valuable reporters on the ensemble of ${\varphi }$ angles sampled by each residue.

A Comparison of the Conformational Distributions of the Achiral Symmetric Liquid Crystal Dimer CB7CB in the Achiral Nematic and Chiral Twist-Bend Nematic Phases

The sets of residual dipolar couplings between carbon and hydrogen nuclei obtained from the proton-encoded (13)C 2D NMR experiment are used to investigate the conformational changes which occur when the achiral symmetric liquid crystal dimer CB7CB changes from the achiral nematic to the chiral twist-bend nematic phase. It is found that these changes are a consequence of the chirality of the twist-bend nematic phase, rather than being the driving force for the stability of this phase.

Modulating alignment of membrane proteins in liquid-crystalline and oriented gel media by changing the size and charge of phospholipid bicelles

We demonstrate that alignment of a structured peptide or small protein solubilized in mixed phospholipid:detergent micelles or bicelles, when embedded in a compressed gel or liquid crystalline medium, can be altered by either changing the phospholipid aggregate shape, charge, or both together. For the hemagglutinin fusion peptide solubilized in bicelles, we show that bicelle shape and charge do not change its helical hairpin structure but impact its alignment relative to the alignment medium, both in charged compressed acrylamide gel and in liquid crystalline d(GpG). The method can be used to generate sets of residual dipolar couplings that correspond to orthogonal alignment tensors, and holds promise for high-resolution structural refinement and dynamic mapping of membrane proteins.

Do the molecules which form discotic liquid crystals have disc‐like structures? The conformation of a simple model compound, 1,2‐dihydroxydiacetylbenzene, determined from the NMR spectra of samples dissolved in liquid crystalline solvents

Many discotic mesogens are molecules with a central aromatic ring with adjacent alkylcarboxylate substituents. The simplest such molecule, 1,2‐dihydroxydiacetylbenzene, which is not mesogenic, is studied by NMR spectroscopy as a solute in a nematic solvent. The spectra are analysed to give sets of residual dipolar couplings, Dij , which are then used to test models for the conformation adopted by the acetate side groups. The conformations and geometry of an isolated molecule are calculated by the ab initio MP2/6‐311G method and also by the DFT approach using the B3LYP functional with the 6‐311++G** basis set. The quantum chemical calculations find that the minimum energy conformer has the acetate groups rotated in opposite directions out of the ring plane, and this kind of structure is also consistent with the NMR data.

Insight into Disulfide Bond Catalysis in Chlamydia from the Structure and Function of DsbH, a Novel Oxidoreductase

The Chlamydia family of human pathogens uses outer envelope proteins that are highly cross-linked by disulfide bonds but nevertheless keeps an unusually high number of unpaired cysteines in its secreted proteins. To gain insight into chlamydial disulfide bond catalysis, the structure, function, and substrate interaction of a novel periplasmic oxidoreductase, termed DsbH, were determined. The structure of DsbH, its redox potential of -269 mV, and its functional properties are similar to thioredoxin and the C-terminal domain of DsbD, i.e. characteristic of a disulfide reductase. As compared with these proteins, the two central residues of the DsbH catalytic motif (CMWC) shield the catalytic disulfide bond and are selectively perturbed by a peptide ligand. This shows that these oxidoreductase family characteristic residues are not only important in determining the redox potential of the catalytic disulfide bond but also in influencing substrate interactions. For DsbH, three functional roles are conceivable; that is, reducing intermolecular disulfides between proteins and small molecules, keeping a specific subset of exported proteins reduced, or maintaining the periplasm of Chlamydia in a generally reducing state.

Visualizing spatially correlated dynamics that directs RNA conformational transitions

RNAs fold into three-dimensional (3D) structures that subsequently undergo large, functionally important, conformational transitions in response to a variety of cellular signals. RNA structures are believed to encode spatially tuned flexibility that can direct transitions along specific conformational pathways. However, this hypothesis has proved difficult to examine directly because atomic movements in complex biomolecules cannot be visualized in 3D by using current experimental methods. Here we report the successful implementation of a strategy using NMR that has allowed us to visualize, with complete 3D rotational sensitivity, the dynamics between two RNA helices that are linked by a functionally important trinucleotide bulge over timescales extending up to milliseconds. The key to our approach is to anchor NMR frames of reference onto each helix and thereby directly measure their dynamics, one relative to the other, using 'relativistic' sets of residual dipolar couplings (RDCs). Using this approach, we uncovered super-large amplitude helix motions that trace out a surprisingly structured and spatially correlated 3D dynamic trajectory. The two helices twist around their individual axes by approximately 53 degrees and 110 degrees in a highly correlated manner (R = 0.97) while simultaneously (R = 0.81-0.92) bending by about 94 degrees. Remarkably, the 3D dynamic trajectory is dotted at various positions by seven distinct ligand-bound conformations of the RNA. Thus even partly unstructured RNAs can undergo structured dynamics that directs ligand-induced transitions along specific predefined conformational pathways.

Mapping the Conformational Landscape of Urea-Denatured Ubiquitin Using Residual Dipolar Couplings

Characterization of the unfolded state is a fundamental prerequisite for understanding protein stability and folding. We have investigated local conformational sampling in urea-denatured ubiquitin at pH 2.5 by measuring an extensive set of residual dipolar couplings (RDCs) under conditions of partial molecular alignment. Seven experimental RDCs per peptide unit, including complementary fixed-geometry and interproton (1H(N)-1H(N) and 1H(N)-1H(alpha)) couplings, were used to investigate the structural properties of the peptide chain. Amino-acid-specific potentials that simultaneously reproduce all RDCs in the molecule are found to sample more extended conformations than the standard statistical coil description. Analysis of 3J(HNHalpha) scalar couplings measured under the same conditions suggests that neither polyproline II nor extended beta-regions dominate this additional sampling of extended conformations. Using this approach we propose a model of the conformational landscape throughout the peptide chain of urea-denatured ubiquitin, providing an essential description for understanding the unfolded state.

Limits on Variations in Protein Backbone Dynamics from Precise Measurements of Scalar Couplings

3JHN,Halpha, 3JHN,Cbeta, and 3JHN,C' couplings, all related to the backbone torsion angle phi, were measured for the third immunoglobulin binding domain of protein G, or GB3. Measurements were carried out using both previously published methods and novel sequences based on the multiple-quantum principle, which limit attenuation of experimental couplings caused by finite lifetimes of the spin states of passive spins. High reproducibility between the multiple-quantum and conventional approaches confirms the accuracy of the measurements. With few exceptions, close agreement between 3JHN,Halpha, 3JHN,Cbeta, and 3JHN,C' and values predicted by their respective Karplus equations is observed. For the three types of couplings, up to 20% better agreement is obtained when fitting the experimental couplings to a dynamic ensemble NMR structure, which has a phi angle root-mean-square spread of 9 +/- 4 degrees and was previously calculated on the basis of a very extensive set of residual dipolar couplings, than for any single static NMR structure. Fits of 3J couplings to a 1.1-A X-ray structure, with hydrogens added in idealized positions, are 40-90% worse. Approximately half of the improvement when fitting to the NMR structures relates to the amide proton deviating from its idealized, in-peptide-plane position, indicating that the positioning of hydrogens relative to the backbone atoms is one of the factors limiting the accuracy at which the backbone torsion angle phi can be extracted from 3J couplings. Introducing an additional, residue-specific variable for the amplitude of phi angle fluctuations does not yield a statistically significant improvement when fitting to a set of dynamic Karplus curves, pointing to a homogeneous behavior of these amplitudes.

Validating the use of database potentials in protein structure determination by NMR

The refinement of protein structures determined by nuclear magnetic resonance spectroscopy against database potentials of mean force allows for the exclusion of unfavourable conformations of the protein backbone during a structure calculation, resulting in protein structures with a marked improvement in Ramachandran statistics. In this communication, we use multiple sets of residual dipolar couplings as quality assessment criteria for several proteins and show that not only do the Ramachandran and structural quality statistics improve, but a significant improvement in the accuracy of structures is achieved upon refinement.

Measurement of Ribose Carbon Chemical Shift Tensors for A-form RNA by Liquid Crystal NMR Spectroscopy

Incomplete motional averaging of chemical shift anisotropy upon weak alignment of nucleic acids and proteins in a magnetic field results in small changes in chemical shift. Knowledge of nucleus-specific chemical shift (CS) tensor magnitudes and orientations is necessary to take full advantage of these measurements in biomolecular structure determination. We report the determination by liquid crystal NMR of the CS tensors for all ribose carbons in A-form helical RNA, using a series of novel 3D NMR pulse sequences for accurate and resolved measurement of the ribose (13)C chemical shifts. The orientation of the riboses relative to the rhombic alignment tensor of the molecule studied, a stem-loop sequence corresponding to helix-35 of 23S rRNA, is known from an extensive set of residual dipolar couplings (RDC), previously used to refine its structure. Singular-value-decomposition fits of the chemical shift changes to this structure, or alternatively to a database of helical RNA X-ray structures, provide the CS tensor for each type of carbon. Quantum chemical calculations complement the experimental results and confirm that the most shielded tensor component lies approximately along the local carbon-oxygen bond axis in all cases and that shielding anisotropy for C3' and C4' is much larger than for C1' and C2', with C5' being intermediate.

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings--an assessment of the interrelation of NMR restraints.

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.

Overall structure and sugar dynamics of a DNA dodecamer from homo- and heteronuclear dipolar couplings and 31P chemical shift anisotropy.

The solution structure of d(CGCGAATTCGCG)(2) has been determined on the basis of an exceptionally large set of residual dipolar couplings. In addition to the heteronuclear (13)C-(1)H and (15)N-(1)H and qualitative homonuclear (1)H-(1)H dipolar couplings, previously measured in bicelle medium, more than 300 quantitative (1)H-(1)H and 22 (31)P-(1)H dipolar restraints were obtained in liquid crystalline Pf1 medium, and 22 (31)P chemical shift anisotropy restraints. High quality DNA structures can be obtained solely on the basis of these new restraints, and these structures are in close agreement with those calculated previously on the basis of (13)C-(1)H and (15)N-(1)H dipolar couplings. In the newly calculated structures, (31)P-(1)H dipolar and (3)JsubH3(')Psub couplings and (31)P CSA data restrain the phosphodiester backbone torsion angles. The final structure represents a quite regular B-form helix with a modest bending of approximately 10 degrees, which is essentially independent of whether or not electrostatic terms are used in the calculation. Combined, the number of homo- and heteronuclear dipolar couplings significantly exceeds the number of degrees of freedom in the system. Results indicate that the dipolar coupling data cannot be fit by a single structure, but are compatible with the presence of rapid equilibria between C2(')-endo and C3(')-endo deoxyribose puckers (sugar switching). The C2(')-H2(')/H2(") dipolar couplings in B-form DNA are particularly sensitive to sugar pucker and yield the largest discrepancies when fit to a single structure. To resolve these discrepancies, we suggest a simplified dipolar coupling analysis that yields N/S equilibria for the ribose sugar puckers, which are in good agreement with previous analyses of NMR J(HH) couplings, with a population of the minor C3(')-endo form higher for pyrimidines than for purines.