Solution structure of diastereomeric thymidine 3'-O(N,N-diisopropyl-O-beta-cyanoethyl-phosphoramidothioate)s by NMR spectroscopy

The conformational properties of a diastereomeric nucleoside-phosphorothioate pair is reported as reflected by different NMR parameters. Configurational assessment is based partly on the different NOE (nuclear overhauser enhancement) effects of the individual isomers and on the trends observed in other NMR parameters. Vicinal carbon-phosphorus couplings reflect the predominance of the epsilon(-) conformation for the Sp isomer and the epsilon(t) conformation for the Rp isomer. The effects of solvent and temperature on these vicinal couplings are described and the results are interpreted in terms of conformational shift towards the preferred epsilon(t) conformation.     

Orienting domains in proteins using dipolar couplings measured by liquid-state NMR: differences in solution and crystal forms of maltodextrin binding protein loaded with beta-cyclodextrin

Protein function is often regulated by conformational changes that occur in response to ligand binding or covalent modification such as phosphorylation. In many multidomain proteins these conformational changes involve reorientation of domains within the protein. Although X-ray crystallography can be used to determine the relative orientation of domains, the crystal-state conformation can reflect the effect of crystal packing forces and therefore may differ from the physiologically relevant form existing in solution. Here we demonstrate that the solution-state conformation of a multidomain protein can be obtained from its X-ray structure using an extensive set of dipolar couplings measured by triple-resonance multidimensional NMR spectroscopy in weakly aligning solvent. The solution-state conformation of the 370-residue maltodextrin-binding protein (MBP) loaded with beta-cyclodextrin has been determined on the basis of one-bond (15)N-H(N), (15)N-(13)C', (13)C(alpha)-(13)C', two-bond (13)C'-H(N), and three-bond (13)C(alpha)-H(N) dipolar couplings measured for 280, 262, 276, 262, and 276 residues, respectively. This conformation was generated by applying hinge rotations to various X-ray structures of MBP seeking to minimize the difference between the experimentally measured and calculated dipolar couplings. Consistent structures have been derived in this manner starting from four different crystal forms of MBP. The analysis has revealed substantial differences between the resulting solution-state conformation and its crystal-state counterpart (Protein Data Bank accession code 1DMB) with the solution structure characterized by an 11(+/-1) degrees domain closure. We have demonstrated that the precision achieved in these analyses is most likely limited by small uncertainties in the intradomain structure of the protein (ca 5 degrees uncertainty in orientation of internuclear vectors within domains). In addition, potential effects of interdomain motion have been considered using a number of different models and it was found that the structures derived on the basis of dipolar couplings accurately represent the effective average conformation of the protein.     

Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.     

Sucrose in aqueous solution revisited, Part 1: molecular dynamics simulations and direct and indirect dipolar coupling analysis

Although the crystal structure of the disaccharide sucrose was solved more than 30 years ago, its conformational distribution in aqueous solution is still a matter of debate. We report here a variety of molecular dynamics simulations (mostly of 100 ns) using the GLYCAM06 force field and various water models, paying particular attention to comparisons to NMR measurements of residual dipolar couplings and electron-mediated spin-spin couplings. We focus on the glycosidic linkage conformation, the puckering phase angle of the fructose ring, and intramolecular hydrogen bonds between the two sugars. Our results show that sucrose is indeed a dynamic molecule, but the crystal conformation is qualitatively the dominant one in dilute solution. A second conformational basin, populated in many force fields, is probably overstabilized in the calculations.     

Global structure of a DNA three-way junction by solution NMR: towards prediction of 3H fold

Three-way junctions (3H) are the simplest and most commonly occurring branched nucleic acids. They consist of three double helical arms (A to C), connected at the junction point, with or without a number of unpaired bases in one or more of the three different strands. Three-way junctions with two unpaired bases in one strand (3HS2) have a high tendency to adopt either of two alternative stacked conformations in which two of the three arms A, B and C are coaxially stacked, i.e. A/B-stacked or A/C-stacked. Empirical stacking rules, which successfully predict for DNA 3HS2 A/B-stacking preference from sequence, have been extended to A/C-stacked conformations. Three novel DNA 3HS2 sequences were designed to test the validity of these extended stacking rules and their conformational behavior was studied by solution NMR. All three show the predicted A/C-stacking preference even in the absence of multivalent cations. The stacking preference for both classes of DNA 3HS2 can thus be predicted from sequence. The high-resolution NMR solution structure for one of the stacked 3HS2 is also reported. It shows a well-defined local and global structure defined by an extensive set of classical NMR restraints and residual dipolar couplings. Analysis of its global conformation and that of other representatives of the 3H family, shows that the relative orientations of the stacked and non-stacked arms, are restricted to narrow regions of conformational space, which can be understood from geometric considerations. Together, these findings open up the possibility of full prediction of 3HS2 conformation (stacking and global fold) directly from sequence.     

Complex glycosaminoglycans: profiling substitution patterns by two-dimensional nuclear magnetic resonance spectroscopy

Biological and pharmacological interactions of heparin and structurally related glycosaminoglycans (GAGs) such as heparan sulfate (HS) involve complex sequences of variously sulfated uronic acid and aminosugar residues. Due to their structural microheterogeneity, these sequences are usually characterized in statistical terms, by high-performance liquid chromatographic analysis of fragments obtained by enzymatic or chemical degradation. Nuclear magnetic resonance (NMR) spectroscopy is also currently used for structural characterization of GAGs. However, the use of monodimensional NMR analysis of complex GAGs is often limited by severe signal overlap that does not allow reliable quantitative measurements. Using magnetically equivalent signals, the higher resolution achieved by two-dimensional NMR methods could be also exploited for quantitative applications. In this work, heteronuclear single quantum coherence (HSQC) spectroscopy has been evaluated to determine variously substituted monosaccharide components of HS and HS mimics obtained by chemical modification of the Escherichia coli K5 polysaccharide (K5-PS) structurally related to the common biosynthetic precursor of heparin and HS. Heparin was used as a model for assessing the influence of 1H-13C spin-spin couplings on "volumes" of the corresponding signals. For major signals, the HSQC approach permitted quantification of additional structural features both in heparins and in a typical HS. The method was applied to profile the substitution patterns of K5-PS derivatives involving different degrees of N,O-sulfation and N-acetylation, including O-sulfated heparosans bearing free amino groups.     

Conformation and dynamics of an RNA internal loop

The conformation and the dynamics of an RNA oligonucleotide (26 nucleotides) which is a model for loop E in eukaryotic 5S RNA have been investigated by one- and two-dimensional NMR. The central portion of the oligonucleotide contains two G A oppositions, a common feature of ribosomal RNAs. The exchangeable proton spectrum indicates that an internal loop separates two stems of four and five base pairs. This observation is not consistent with structures for loop E containing mismatched G.A base pairs proposed from chemical and enzymatic studies on Xenopus laevis 5S RNA. The nonexchangeable proton spectrum has been assigned by two-dimensional NMR. Scalar couplings from correlated experiments and interproton distances from NOESY experiments at short mixing times have been used to determine glycosidic angles, sugar puckers, and other conformational features. The conformation of the stems is very close to standard A-form RNA, and extensive base stacking continues into the internal loop. This result provides a structural basis for the large favorable enthalpy of duplex formation determined in thermodynamic studies. Unusual structural and dynamic features are localized in the nucleotides connecting the loop to the stems.     

Protein backbone motions viewed by intraresidue and sequential HN–Hα residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.     

Structures of heparin-derived disaccharide bound to cobra cardiotoxins: context-dependent conformational change of heparin upon binding to the rigid core of the three-fingered toxin

Glycosaminoglycans (GAGs) have been suggested to be a potential target for cobra cardiotoxin (CTX) with high affinity and specificity via a cationic belt at the concave surface of the polypeptide. The interaction of GAGs, such as high-molecular weight heparin, with CTXs not only can induce aggregation of CTX molecules but also can enhance their penetration into membranes. The binding of short chain heparin, such as a heparin-derived disaccharide [DeltaUA2S(1-->4)-alpha-D-GlcNS6S], to CTX A3 from Taiwan cobra (Naja atra), however, will not induce aggregation and was, therefore, investigated by high-resolution (1)H NMR. A novel heparin binding site on the convex side of the CTX, near the rigid disulfide bond-tightened core region of Cys38, was identified due to the observation of intermolecular NOEs between the protein and carbohydrate. The derived carbohydrate conformation using complete relaxation and conformational exchange matrix analysis (CORCEMA) of NOEs indicated that the glycosidic linkage conformation and the ring conformation of the unsaturated uronic acid in the bound state depended significantly on the charge context of CTX molecules near the binding site. Specifically, comparative binding studies of several heparin disaccharide homologues with two CTX homologues (CTX Tgamma from Naja nigricollis and CTX A3) indicated that the electrostatic interaction of N-sulfate of glucosamine with NH(3)(+)zeta of Lys12 and of the 2-O-sulfate of the unsaturated uronic acid with NH(3)(+)zeta of Lys5 played an important role. These results also suggest a model on how the CTX-heparin interaction may regulate heparin-induced aggregation of the toxin via the second heparin binding site.     

The use of dipolar couplings for determining the solution structure of rat apo-S100B(betabeta)

The relative orientations of adjacent structural elements without many well-defined NOE contacts between them are typically poorly defined in NMR structures. For apo-S100B(betabeta) and the structurally homologous protein calcyclin, the solution structures determined by conventional NMR exhibited considerable differences and made it impossible to draw unambiguous conclusions regarding the Ca2+-induced conformational change required for target protein binding. The structure of rat apo-S100B(betabeta) was recalculated using a large number of constraints derived from dipolar couplings that were measured in a dilute liquid crystalline phase. The dipolar couplings orient bond vectors relative to a single-axis system, and thereby remove much of the uncertainty in NOE-based structures. The structure of apo-S100B(betabeta) indicates a minimal change in the first, pseudo-EF-hand Ca2+ binding site, but a large reorientation of helix 3 in the second, classical EF-hand upon Ca2+ binding.     

The closed and compact domain organization of the 70-kDa human cytochrome P450 reductase in its oxidized state as revealed by NMR

The NADPH cytochrome P450 reductase (CPR), a diflavin enzyme, catalyzes the electron transfer (ET) from NADPH to the substrate P450. The crystal structures of mammalian and yeast CPRs show a compact organization for the two domains containing FMN (flavin mononucleotide) and FAD (flavin adenine dinucleotide), with a short interflavin distance consistent with fast ET from the NADPH-reduced FAD to the second flavin FMN. This conformation, referred as "closed", contrasts with the alternative opened or extended domain arrangements recently described for partially reduced or mutant CPR. Internal domain flexibility in this enzyme is indeed necessary to account for the apparently conflicting requirements of having FMN flavin accessible to both the FAD and the substrate P450 at the same interface. However, how interdomain dynamics influence internal and external ETs in CPR is still largely unknown. Here, we used NMR techniques to explore the global, domain-specific and residue-specific structural and dynamic properties of the nucleotide-free human CPR in solution in its oxidized state. Based on the backbone resonance assignment of this 70-kDa protein, we collected residue-specific (15)N relaxation and (1)H-(15)N residual dipolar couplings. Surprisingly and in contrast with previous studies, the analysis of these NMR data revealed that the CPR exists in a unique and predominant conformation that highly resembles the closed conformation observed in the crystalline state. Based on our findings and the previous observations of conformational equilibria of the CPR in partially reduced states, we propose that the large-scale conformational transitions of the CPR during the catalytic cycle are tightly controlled to ensure optimal electron delivery.     

Physicochemical studies of the protein-lipid interactions in melittin-containing micelles

Complexes of melittin with detergents and phospholipids have been characterized by fluorescence, circular dichroism, ultracentrifugation, quasi-elastic light scattering and 1H nuclear magnetic resonance (NMR) experiments. By ultracentrifugation and quasi-elastic light-scattering measurements it is shown that melittin forms stoichiometrically well-defined complexes with dodecylphosphocholine micelles consisting of one melittin molecule and approximately forty detergent molecules. Evidence from fluorescence, circular dichroism and 1H nuclear magnetic resonance experiments indicates that the conformation of melittin bound to micelles of various detergents or of diheptanoyl phosphatidylcholine is largely independent of the type of lipid and furthermore appears to be quite closely related to the conformation of melittin bound to phosphatidylcholine bilayers. 1H NMR is used to investigate the conformation of micelle-bound melittin in more detail and to compare certain aspects of the melittin conformation in the micelles with the spatial structures of monomeric and self-aggregated tetrameric melittin in aqueous solution. The experience gained with this system demonstrates that high resolution NMR of complexes of membrane proteins with micelles provides a viable method for conformational studies of membrane proteins.     

Contribution of NMR spectroscopy to the study of structure-function relations of proteins and peptides

NMR spectroscopy provides a unique means to study molecular conformation, mechanisms of action and structure-function relationships for peptides and proteins in solution under conditions approaching those of their physiological environment. Development of NMR techniques, especially directed to the peptide and protein conformational analysis, is considered under the topics of two-level signal assignment and structural significance of homo- and heteronuclear spin-spin couplings. The results of NMR conformational analysis are presented for solution spatial structure of valinomy cin and gramicidin A antibiotics, honey-bee neurotoxin apamin, scorpion insectotoxin I5A and snake venom neurotoxins of "short" and "long" types. The structure-function relationships are discussed for these biologically active molecules.     

Sign determination of dipolar couplings in field-oriented bicelles by variable angle sample spinning (VASS)

Residual dipolar couplings are being increasingly used as structural constraints for NMR studies of biomolecules. A problem arises when dipolar coupling contributions are larger than scalar contributions for a given spin pair, as is commonly observed in solid state NMR studies, in that signs of dipolar couplings cannot easily be determined. Here the sign ambiguities of dipolar couplings in field-oriented bicelles are resolved by variable angle sample spinning (VASS) techniques. The director behavior of field-oriented bicelles (DMPC/DHPC, DMPC/CHAPSO) in VASS is studied by ³¹P NMR. A stable configuration occurs when the spinning angle is smaller than the magic angle, 54.7°, and the director (or bicelle normal) of the disks is mainly distributed in a plane perpendicular to the rotation axis. Since the dipolar couplings depend on how the bicelles are oriented with respect to the magnetic field, it is shown that the dipolar interaction can be scaled to the same order as the J-coupling by moving the spinning axis from 0° toward 54.7°. Thus the relative sign of dipolar and scalar couplings can be determined.     

Concerted motions in HIV-1 TAR RNA may allow access to bound state conformations: RNA dynamics from NMR residual dipolar couplings.

Ground-state dynamics in RNA is a critical precursor for structural adaptation observed ubiquitously in protein-RNA recognition. A tertiary conformational analysis of the stem-loop structural element in the transactivation response element (TAR) from human immunodeficiency virus type 1 (HIV-I) RNA is presented using recently introduced NMR methods that rely on the measurement of residual dipolar couplings (RDC) in partially oriented systems. Order matrix analysis of RDC data provides evidence for inter-helical motions that are of amplitude 46(+/-4) degrees, of random directional character, and that are executed about an average conformation with an inter-helical angle between 44 degrees and 54 degrees. The generated ensemble of TAR conformations have different organizations of functional groups responsible for interaction with the trans-activator protein Tat, including conformations similar to the previously characterized bound-state conformation. These results demonstrate the utility of RDC-NMR for simultaneously characterizing RNA tertiary dynamics and average conformation, and indicate an avenue for TAR complex formation involving tertiary structure capture.     

Conformational studies of Lewis X and Lewis A trisaccharides using NMR residual dipolar couplings.

The conformations of the histo-blood group carbohydrate antigens Lewis X (Le(x)) and Lewis A (Le(a)) were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. A strategy for rapid calculation of the difference between theoretical and experimental dipolar couplings of a large number of model structures generated by computer simulations was developed, resulting in an accurate model structure for the compounds. Monte Carlo simulations were used to generate models for the trisaccharides, and orientations of each model were sought that could reproduce the experimental residual dipolar coupling values. For both, Le(a) and Le(x), single low energy models giving excellent agreement with experiment were found, implying a compact rigidly folded conformation for both trisaccharides. The new approach was also applied to the pentasaccharides lacto-N-fucopentaose 2 (LNF-2) and lacto-N-fucopentaose 3 (LNF-3) proving its consistency and robustness. For describing the conformation of tightly folded oligosaccharides, a definition for characterization of ring planes in pyranoside chairs is proposed and applied to the analysis of the relation between the fucose and galactose residues in the epitopes, revealing the structural similarity between them.     

The 1,3‐diyne linker as a rigid “i,i+7” staple for α‐helix stabilization: Stereochemistry at work

Short alphahelical peptide sequences were stabilized through Glaser‐Hay couplings of propargylated l ‐ and/or d ‐serine residues at positions i and i+7. NMR analysis confirmed a full stabilization of the helical structure when a d ‐Ser (i), l ‐Ser (i+7) combination was applied. In case two l ‐Ser residues were involved in the cyclization, the helical conformation is disrupted outside the peptide's macrocycle.     

Non-Watson-Crick structures in oligodeoxynucleotides: self-association of d(TpCpGpA) stabilized at acidic pH

The 1H NMR spectrum of the tetradeoxynucleotide d(TpCpGpA) was examined as a function of temperature, pH, and concentration. At pH 7 and above the solution conformation for this oligodeoxynucleotide appears to be a mixture of random coil and Watson-Crick duplex. At 25 degrees C, a pH titration of d(TpCpGpA) shows that distinct conformational changes occur as the pH is lowered below 7.0. These conformational changes are reversible upon readjusting the pH to neutrality, indicating the presence of a pH-dependent set of conformational equilibria. At 25 degrees C, the various conformational states in the mixture are in rapid exchange on the NMR time scale. Examination of the titration curve shows the presence of distinct conformational states at pH greater than 7, and between pH 4 and pH 5. At pH less than 4, a third conformational state is present. When the pH titration is repeated at 5 degrees C, the conformational equilibria are in slow exchange on the NMR time scale; distinct signals from each conformational state are observable. The stable conformational state present between pH 4 and pH 5 represents an ordered conformation of d(TpCpGpA) which dissociates to a less ordered structure upon raising the temperature. This ordered conformation does not result from an intramolecular rearrangement, as is shown by by spectra obtained by varying oligodeoxynucleotide concentration at constant pH. The ordered conformation differs from the Watson-Crick helix, as is shown from nuclear Overhauser enhancement experiments, as well as chemical shift data. An ordered conformation for d(TpCpGpA) was previously reported [Reid, D. G., Salisbury, S. A., Brown, T., & Williams, D. H. (1985) Biochemistry 24, 4325-4332].(ABSTRACT TRUNCATED AT 250 WORDS)     

Graphical representation of the salient conformational features of protein residues.

A composite plot for depicting in two dimensions the conformation and the secondary structural features of protein residues has been developed. Instead of presenting the exact values of the main- and side-chain torsion angles (φ, psi and chi(1)), it indicates the region in the three-dimensional conformational space to which a residue belongs. Other structural aspects, like the presence of a cis peptide bond and disulfide linkages, are also displayed. The plot may be used to recognize patterns in the backbone and side-chain conformation along a polypeptide chain and to compare protein structures derived from X-ray crystallography, NMR spectroscopy or molecular modelling studies and also to highlight the effect of mutation on structure.     

Virtual and solution conformations of oligosaccharides

The possibility that observed nuclear Overhauser enhancements and bulk longitudinal relaxation times, parameters measured by 1H NMR and often employed in determining the preferred solution conformation of biologically important molecules, are the result of averaging over many conformational states is quantitatively evaluated. Of particular interest was to ascertain whether certain 1H NMR determined conformations are "virtual" in nature; i.e., the fraction of the population of molecules actually found at any time within the subset of conformational space defined as the "solution conformation" is vanishingly small. A statistical mechanics approach was utilized to calculate an ensemble average relaxation matrix from which (NOE)'s and (T1)'s are calculated. Model glycosidic linkages in four oligosaccharides were studied. The solution conformation at any glycosidic linkage is properly represented by a normalized, Boltzmann distribution of conformers generated from an appropriate potential energy surface. The nature of the resultant population distributions is such that 50% of the molecular population is found within 1% of available microstates, while 99% of the molecular population occupies about 10% of the ensemble microstates, a number roughly equal to that sterically allowed. From this analysis we conclude that in many cases quantitative interpretation of NMR relaxation data, which attempts to define a single set of allowable torsion angle values consistent with the observed data, will lead to solution conformations that are either virtual or reflect torsion angle values possessed by a minority of the molecular population. On the other hand, calculation of ensemble average NMR relaxation data yields values in agreement with experimental results. Observed values of NMR relaxation data are the result of the complex interdependence of the population distribution and NOE (or T1) surfaces in conformational space. In conformational analyses, NMR data can therefore be used to test different population distributions calculated from empirical potential energy functions.