Use Of Residual Dipolar Couplings Research Articles

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

TITANIA: Model-Free Interpretation of Residual Dipolar Couplings in the Context of Organic Compounds.

Residual dipolar couplings (RDCs) become increasingly important as additional NMR parameters in the structure elucidation of organic compounds but are usually used in fitting procedures to discriminate between (computed) structures that are in accordance with RDCs and others that can be ruled out. Thus, the determination of configurations requires prior structural information. The direct use of RDCs as restraints to construct structures based on RDCs has only recently begun also in organic compounds. No protocol has been published though that uses the vector and dynamics information available in multialignment data sets directly for the joint determination of conformation and configuration of organic compounds. This is proposed in the current study. We show that by employing these data, even a flat or random start structure converges into the correctly configured structure when employing multiple alignment data sets in our iterative procedure. The requirements in terms of the number of RDCs and alignment media necessary are discussed in detail.

Determination of Configuration and Conformation of a Reserpine Derivative with Seven Stereogenic Centers Using Molecular Dynamics with RDC-Derived Tensorial Constraints*.

NMR‐based determination of the configuration of complex molecules containing many stereocenters is often not possible using traditional NOE data and coupling patterns. Making use of residual dipolar couplings (RDCs), we were able to determine the relative configuration of a natural product containing seven stereocenters, including a chiral amine lacking direct RDC data. To identify the correct relative configuration out of 32 possible ones, experimental RDCs were used in three different approaches for data interpretation: by fitting experimental data based singular value decomposition (SVD) using a single alignment tensor and either (i) a single conformer or (ii) multiple conformers, or alternatively (iii) using molecular dynamics simulations with tensorial orientational constraints (MDOC). Even though in all three approaches one and the same configuration could be selected and clear discrimination between possible configurations was achieved, the experimental data was not fully satisfied by the methods based on single tensor approaches. While these two approaches are faster, only MDOC is able to fully reproduce experimental results, as the obtained conformational ensemble adequately covers the conformational space necessary to describe the molecule with inherent flexibility.

Configurational analysis by residual dipolar couplings: A critical assessment of diastereomeric differentiabilities.

Two independent statistical models for evaluating the certainties of configurational assignments of compounds based on nuclear magnetic resonance (NMR) data are evaluated and compared. Both methods yield weights or probabilities with which two or more structure models (constitutional or configurational isomers or even conformers) could be differentiated based on experimental parameters. Although this paper focusses on the use of residual dipolar couplings (RDCs) for the differentiation of diastereomers, the concept can be expanded to any set of experimental NMR-derived parameters. It is demonstrated that highly reliable configurational assignments crucially must depend on thorough statistical analysis, which is frequently neglected in the literature.

The use of residual dipolar couplings for conformational analysis of non-steroidal anti-inflammatory drugs dissolved in weakly ordering media

ABSTRACTThe three-dimensional structure and conformational equilibrium of flexible bioactive molecules in solution are of considerable importance in clarifying their role in chemical and biological activity, since their action takes place in a liquid medium. The use of residual dipolar couplings (RDCs) recently emerged as a promising and effective method to address this kind of investigation. In this brief review, the conformational analysis of several flexible non-steroidal anti-inflammatory drugs, belonging to the class of salicylates and α-arylpropionic acids (profens), is described using the RDCs measured in weakly aligning poly-γ-benzyl-l-glutamate solutions. By combining these experimental parameters, extracted by using different mono- and bi-dimensional nuclear magnetic resonance experiments, with the theoretical additive potential-direct probability description approach, accurate descriptions of the probability distributions have been obtained for all the drug molecules selected.

Use of residual dipolar couplings in conformational analysis of meta-disubstituted calix[4]arenes.

Dimercuration of tetrapropoxy calix[4]arene followed by a reaction with isoamyl nitrite afforded dinitroso regioisomers with unique substitution patterns. The unusual conformational behaviour of these inherently chiral systems was revealed by the combination of dynamic NMR and residual dipolar coupling (RDC) techniques.

Characterization of the Interdomain Motions in Hen Lysozyme Using Residual Dipolar Couplings as Replica-Averaged Structural Restraints in Molecular Dynamics Simulations

Hen lysozyme is an enzyme characterized by the presence of two domains whose relative motions are involved in the mechanism of binding and release of the substrates. By using residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations, we characterize the breathing motions describing the interdomain fluctuations of this protein. We found that the ensemble of conformations that we determined spans the entire range of structures of hen lysozyme deposited in the Protein Data Bank, including both the free and bound states, suggesting that the thermal motions in the free state provide access to the structures populated upon binding. The approach that we present illustrates how the use of residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations makes it possible to explore conformational fluctuations of a relatively large amplitude in proteins.

Review: use of residual dipolar couplings to determine the structure of carbohydrates

Solution nuclear magnetic resonance spectroscopy is especially useful in the carbohydrate field. The measurement of residual dipolar couplings provides long-range structural information, a valuable complement for the structural study of carbohydrates either in its free form or in the bound state to proteins. They permit to deduce the geometry and the flexibility of the glycosidic linkages, which have a major influence on the conformation of carbohydrates and their overall shape. This article reviews the current application of the residual dipolar couplings methodology to carbohydrates.

MSpin‐RDC. A program for the use of residual dipolar couplings for structure elucidation of small molecules

We developed a new program, MSpin-RDC for the analysis of residual dipolar coupling data. This software, specially designed for small molecule analysis, can directly read many molecular-modeling and popular chemistry file formats and accept RDC values as a simple free-format table. Alignment tensor can then be computed by singular value decomposition, as well as predicted using inertia and gyration tensor-based methodologies. Trial structures are then ranked according to their Cornilescu's quality factor (Q) values. Analysis of multiconformational problems and fitting of RDC data to relative populations can be accomplished using the single-tensor approximation.

Solution Structure of the IIAChitobiose-HPr Complex of the N,N′-Diacetylchitobiose Branch of the Escherichia coli Phosphotransferase System

The solution structure of the complex of enzyme IIA of the N,N'-diacetylchitobiose (Chb) transporter with the histidine phosphocarrier protein HPr has been solved by NMR. The IIA(Chb)-HPr complex completes the structure elucidation of representative cytoplasmic complexes for all four sugar branches of the bacterial phosphoryl transfer system (PTS). The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state. IIA(Chb)(H89E) and HPr form a weak complex with a K(D) of ~0.7 mM. The interacting binding surfaces, concave for IIA(Chb) and convex for HPr, complement each other in terms of shape, residue type, and charge distribution, with predominantly hydrophobic residues, interspersed by some uncharged polar residues, located centrally, and polar and charged residues at the periphery. The active site histidine of HPr, His-15, is buried within the active site cleft of IIA(Chb) formed at the interface of two adjacent subunits of the IIA(Chb) trimer, thereby coming into close proximity with the active site residue, H89E, of IIA(Chb). A His89-P-His-15 pentacoordinate phosphoryl transition state can readily be modeled without necessitating any significant conformational changes, thereby facilitating rapid phosphoryl transfer. Comparison of the IIA(Chb)-HPr complex with the IIA(Chb)-IIB(Chb) complex, as well as with other cytoplasmic complexes of the PTS, highlights a unifying mechanism for recognition of structurally diverse partners. This involves generating similar binding surfaces from entirely different underlying structural elements, large interaction surfaces coupled with extensive redundancy, and side chain conformational plasticity to optimize diverse sets of intermolecular interactions.

Determination of the Absolute Configuration of 19-OH-(-)-Eburnamonine Using a Combination of Residual Dipolar Couplings, DFT Chemical Shift Predictions, and Chiroptics

19-OH-(-)-eburnamonine 1 is a new indole alkaloid isolated from Bonafusia macrocalyx. A natural derivative from (-)-eburnamonine for which absolute configuration and conformation has been determined by making use of residual dipolar couplings enhanced NMR, circular dichroism spectra and high-level computations.

Combined Use of Residual Dipolar Couplings and Solution X-ray Scattering To Rapidly Probe Rigid-Body Conformational Transitions in a Non-phosphorylatable Active-Site Mutant of the 128 kDa Enzyme I Dimer

The first component of the bacterial phosphotransferase system, enzyme I (EI), is a multidomain 128 kDa dimer that undergoes large rigid-body conformational transitions during the course of its catalytic cycle. Here we investigate the solution structure of a non-phosphorylatable active-site mutant in which the active-site histidine is substituted by glutamine. We show that perturbations in the relative orientations and positions of the domains and subdomains can be rapidly and reliably determined by conjoined rigid-body/torsion angle/Cartesian simulated annealing calculations driven by orientational restraints from residual dipolar couplings and shape and translation information afforded by small- and wide-angle X-ray scattering. Although histidine and glutamine are isosteric, the conformational space available to a Gln side chain is larger than that for the imidazole ring of His. An additional hydrogen bond between the side chain of Gln189 located on the EIN(α/β) subdomain and an aspartate (Asp129) on the EIN(α) subdomain results in a small (∼9°) reorientation of the EIN(α) and EIN(α/β) subdomains that is in turn propagated to a larger reorientation (∼26°) of the EIN domain relative to the EIC dimerization domain, illustrating the positional sensitivity of the EIN domain and its constituent subdomains to small structural perturbations.

Structural Disorder within Sendai Virus Nucleoprotein and Phosphoprotein: Insight into the Structural Basis of Molecular Recognition

Intrinsically disordered regions of significant length are present throughout eukaryotic genomes, and are particularly prevalent in viral proteins. Due to their inherent flexibility, these proteins inhabit a conformational landscape that is too complex to be described by classical structural biology. The elucidation of the role that conformational flexibility plays in molecular function will redefine our understanding of the molecular basis of biological function, and the development of appropriate technology to achieve this aim remains one of the major challenges for the future of structural biology. NMR is the technique of choice for studying intrinsically disordered proteins, providing information about structure, flexibility and interactions at atomic resolution even in completely disordered proteins. In particular residual dipolar couplings (RDCs) are sensitive and powerful tools for determining local and long-range structural behaviour in flexible proteins. Here we describe recent applications of the use of RDCs to quantitatively describe the level of local structure in intrinsically disordered proteins involved in replication and transcription in Sendai virus.

Enhancing the Orienting Properties of Poly(γ‐benzyl‐L‐glutamate) by means of Additives

Residual dipolar couplings (RDCs) have recently become increasingly important in organic structure determination due to their unique information content. One main limitation for the use of RDCs in organic compounds is the orientation that needs to be induced to be able to measure RDCs. So far, there are very few possibilities to modulate the orientational properties of organic solutes and even less when chiral media are considered. Based on our recent findings that the critical concentration of the liquid-crystalline phase of homopolypeptides depends on their molecular weight, we sought for further ways to modulate the orienting properties. We were especially interested in seeing whether we could not only influence the induced degree of orientation, but whether we could also change the solute's preferred orientation and even enhance enantiodifferentiation. We thus tried different aprotic and protic additives and were successful in all of the above-mentioned aspects by using CCl(4) as the additive. Furthermore, we consider DMSO to be a very useful additive. The LC phase of low MW poly(gamma-benzyl-L-glutamate) (PBLG) is usually unstable when DMSO is added. The high MW PBLG used in this study, however, remained stable up to a DMSO/CDCl(3) ratio of 1:2. By using this combination of solvents, the alignment of the two enantiomers of a compound, which is insoluble in CDCl(3), namely, the HCl salt of a tryptophane ester, was possible leading to high-quality spectra. The two enantiomers of the tryptophane ester showed different couplings, thus indicating that enantiodifferentiation is taking place. Thus we were able to modulate the orienting properties (degree of orientation, preferred orientation and enantiodifferentiation) of PBLG by using additives and to increase the accessible solvent and solute range significantly.

Orientational Properties of Poly‐γ‐benzyl‐L‐glutamate: Influence of Molecular Weight and Solvent on Order Parameters of the Solute

Residual dipolar couplings (RDCs) have recently become increasingly important in organic structure determination due to their unique information content. One main limitation for the use of RDCs in organic compounds, however, is that the compound in question needs to be oriented with respect to the magnetic field in order to measure RDCs. So far, there are very few possibilities for modulating the induced degree of orientation. The situation is even worse when chiral orienting media are considered, which could allow absolute configuration determination in the future. We have conducted a systematic investigation into modulating the orientation induced by one chiral orienting medium, namely organic solutions of PBLG (poly-gamma-benzyl-L-glutamate), as a function of its molecular weight and the organic co-solvent used, and have obtained significant insights into factors that influence the order induced. With increasing molecular weight of the polypeptide the orientation of the solutes decreases, leading to well-resolved spectra with improved line shapes. This can be attributed exclusively to the fact that the critical concentration of the liquid-crystalline phase decreases with increasing molecular weight (pure dilution effect). Any influence of increasing flexibility on the orientation can be ruled out.

Use of RDCs in rigid organic compounds and some practical considerations concerning alignment media

AbstractA molecule that is introduced into a liquid crystal or alignment medium forfeits its isotropic tumbling motion. The resulting partial alignment leads to observable anisotropic interactions, dipolar couplings being one. To control the degree of alignment that the molecule experiences is essential to retain the quality of high‐resolution solution‐state spectra. An overview of alignment media is given. The problems associated with very large residual dipolar couplings (RDCs) are pointed out, and means to reduce residual dipolar couplings (RDCs) for strychnine in substituted polyglutamate liquid crystals are described. The use of RDCs for the assignment of diastereotopic protons and for the determination of the relative configuration in rigid compounds is then described in detail, also pointing out some pitfalls. © 2007 Wiley Periodicals, Inc. Concepts Magn Reson Part A 30A: 65–80, 2007.

Simultaneous Determination of Protein Backbone Structure and Dynamics from Residual Dipolar Couplings

Determination of protein structure classically results in a single average configuration that takes no account of conformational fluctuation. Dynamics are, however, inherently linked to structure and crucial to our understanding of biological function. In this study we have used analytical descriptions of dynamic averaging of residual dipolar couplings (RDCs) to simultaneously determine the backbone structure and dynamics of protein GB3. RDCs alone are used to determine an ultrahigh-resolution structure that compares very closely with a refined X-ray structure (rmsd of 0.34 A overall backbone residues). Dynamic amplitudes reporting on motions up to the millisecond time scale reproduce the main characteristics of dynamics previously determined in conjunction with the crystal structure. The use of RDCs alone allows a bias-free comparison with a purely static approach to structure determination. Extensive cross validation clearly demonstrates that the dynamic description is superior to the static approximation. The demonstration that this level of structural resolution and dynamic detail can be extracted from RDCs supports previous indications that these parameters contain extremely precise information about biomolecular conformational sampling.

Solution Structure of a Post-transition State Analog of the Phosphotransfer Reaction between the A and B Cytoplasmic Domains of the Mannitol Transporter IIMannitol of the Escherichia coli Phosphotransferase System

The solution structure of the post-transition state complex between the isolated cytoplasmic A (IIAMtl) and phosphorylated B (phospho-IIBMtl) domains of the mannitol transporter of the Escherichia coli phosphotransferase system has been solved by NMR. The active site His-554 of IIAMtl was mutated to glutamine to block phosphoryl transfer activity, and the active site Cys-384 of IIBMtl (residues of IIBMtl are denoted in italic type) was substituted by serine to permit the formation of a stable phosphorylated form of IIBMtl. The two complementary interaction surfaces are predominantly hydrophobic, and two methionines on IIBMtl, Met-388 and Met-393, serve as anchors by interacting with two deep pockets on the surface of IIAMtl. With the exception of a salt bridge between the conserved Arg-538 of IIAMtl and the phosphoryl group of phospho-IIBMtl, electrostatic interactions between the two proteins are limited to the outer edges of the interface, are few in number, and appear to be weak. This accounts for the low affinity of the complex (Kd approximately 3.7 mm), which is optimally tuned to the intact biological system in which the A and B domains are expressed as a single polypeptide connected by a flexible 21-residue linker. The phosphoryl transition state can readily be modeled with no change in protein-protein orientation and minimal perturbations in both the backbone immediately adjacent to His-554 and Cys-384 and the side chains in close proximity to the phosphoryl group. Comparison with the previously solved structure of the IIAMtl-HPr complex reveals how IIAMtl uses the same interaction surface to recognize two structurally unrelated proteins and explains the much higher affinity of IIAMtl for HPr than IIBMtl.

Three-dimensional Solution Structure of the Cytoplasmic B Domain of the Mannitol Transporter IIMannitol of the Escherichia coli Phosphotransferase System

The solution structure of the cytoplasmic B domain of the mannitol (Mtl) transporter (II(Mtl)) from the mannitol branch of the Escherichia coli phosphotransferase system has been solved by multidimensional NMR spectroscopy with extensive use of residual dipolar couplings. The ordered IIB(Mtl) domain (residues 375-471 of II(Mtl)) consists of a four-stranded parallel beta-sheet flanked by two helices (alpha(1) and alpha(3)) on one face and helix alpha(2) on the opposite face with a characteristic Rossmann fold comprising two right-handed beta(1)alpha(1)beta(2) and beta(3)alpha(2)beta(4) motifs. The active site loop is structurally very similar to that of the eukaryotic protein tyrosine phosphatases, with the active site cysteine (Cys-384) primed in the thiolate state (pK(a) < 5.6) for nucleophilic attack at the phosphorylated histidine (His-554) of the IIA(Mtl) domain through stabilization by hydrogen bonding interactions with neighboring backbone amide groups at positions i + 2/3/4 from Cys-384 and with the hydroxyl group of Ser-391 at position i + 7. Modeling of the phosphorylated state of IIB(Mtl) suggests that the phosphoryl group can be readily stabilized by hydrogen bonding interactions with backbone amides in the i + 2/4/5/6/7 positions as well as with the hydroxyl group of Ser390 at position i + 6. Despite the absence of any significant sequence identity, the structure of IIB(Mtl) is remarkably similar to the structures of bovine protein tyrosine phosphatase (which contains two long insertions relative to IIB(Mtl)) and the cytoplasmic B component of enzyme II(Chb), which fulfills an analogous role to IIB(Mtl) in the N,N'-diacetylchitobiose branch of the phosphotransferase system. All three proteins utilize a cysteine residue in the nucleophilic attack of a phosphoryl group covalently bound to another protein.

High-temperature solution NMR structure of TmCsp.

Cold shock proteins (Csps) are assumed to play a central role in the regulation of gene expression under cold shock conditions. Acting as single-stranded nucleic acid-binding proteins, they trigger the translation process and are therefore involved in the compensation of the influence of low temperatures (cold shock) upon the cell metabolism. However, it is unknown so far how Csps are switched on and off as a function of temperature. The aim of the present study is the study of possible structural changes responsible for this switching process. (1)H-(15)N HSQC spectra recorded at different temperatures and chemical-shift analysis have indicated subtle conformational changes for the cold-shock protein from the hyperthermophilic bacterium Thermotoga maritima (TmCsp) when the temperature is elevated from 303 K to its physiological temperature (343 K). The three-dimensional structure of TmCsp was determined by nuclear magnetic resonance (NMR) spectroscopy at 343 K to obtain quantitative information concerning these structural changes. By use of residual dipolar couplings, the loss of NOE information at high temperature could be compensated successfully. Most pronounced conformational changes compared with room-temperature conditions are observed for amino acid residues closely neighbored to two characteristic beta-bulges and a well-defined loop region of the protein. Because the residues shown to be responsible for the interaction of TmCsp with single-stranded nucleic acids can almost exclusively be found within these regions, nucleic acid-binding activity might be down-regulated with increasing temperature by the described conformational changes.