31P NMR spectra of oligodeoxyribonucleotide duplex lac operator-repressor headpiece complexes: importance of phosphate ester backbone flexibility in protein-DNA recognition.

The 31P NMR spectra of various 14-base-pair lac operators bound to both wild-type and mutant lac repressor headpiece proteins were analyzed to provide information on the backbone conformation in the complexes. The 31P NMR spectrum of a wild-type symmetrical operator, d(TGTGAGCGCTCACA)2, bound to the N-terminal 56-residue headpiece fragment of a Y7I mutant repressor was nearly identical to the spectrum of the same operator bound to the wild-type repressor headpiece. In contrast, the 31P NMR spectrum of the mutant operator, d(TATAGAGCGCTCATA)2, wild-type headpiece complex was significantly perturbed relative to the wild-type repressor-operator complex. The 31P chemical shifts of the phosphates of a second mutant operator, d(TGTGTGCGCACACA)2, showed small but specific changes upon complexation with either the wild-type or mutant headpiece. The 31P chemical shifts of the phosphates of a third mutant operator, d(TCTGAGCGCTCAGA)2, showed no perturbations upon addition of the wild-type headpiece. The 31P NMR results provide further evidence for predominant recognition of the 5'-strand of the 5'-TGTGA/3'-ACACT binding site in a 2:1 protein to headpiece complex. It is proposed that specific, strong-binding operator-protein complexes retain the inherent phosphate ester conformational flexibility of the operator itself, whereas the phosphate esters are conformationally restricted in the weak-binding operator-protein complexes. This retention of backbone torsional freedom in strong complexes is entropically favorable and provides a new (and speculative) mechanism for protein discrimination of different operator binding sites. It demonstrates the potential importance of phosphate geometry and flexibility on protein recognition and binding.     

Assignment of the 1H-NMR spectrum of a lac repressor headpiece-operator complex in H2O and identification of NOEs. Consequences for protein-DNA interaction

A complex between the headpiece amino-terminal residues 1-56 of lac repressor (HP56) and an 11-bp lac operator fragment was studied by 1H NMR. The sequence specific assignment of the exchangeable and non-exchangeable protons has been accomplished. Several protons have favourable chemical shifts in the complex, therefore new intraprotein NOEs could be found that had not been unambigously identified in the free protein. By comparison, most of these intraprotein NOEs are also present in the spectra of the free headpiece but some are different. Furthermore, several new proteins DNA NOEs could be identified. The NOE between the side-chain amide protons of Gln18 and C5H of C7 confirms the specific contact between these residues which was proposed from genetic experiments [Ebright, R. M. (1985) J. Biomol. Struct. & Dyn. 3, 281-297]. The implications of the new data for the interaction between the lac repressor headpiece and its operator are discussed.     

lac Repressor headpiece binds specifically to half of the lac operator: a proton nuclear magnetic resonance study

The complex formation of the N-terminal domain (headpiece) of the Escherichia coli lac repressor and a synthetic 14-base-pair lac operator fragment has been investigated by 1H NMR. Titration shifts in the imino-proton region of the DNA spectrum and in the aromatic region of the headpiece spectrum are examined in detail and interpreted where possible. The assignment of the resonances in the complex follows in part from the titration data and is completed by nuclear Overhauser measurements. The shift of the His-29 C-2 resonance has been used to assess the binding strength of the complex. Evidence is presented for the presence of a high-affinity site on the lac operator fragment (KD less than or equal to 2 X 10(-5) M), which shows features in common with one of the specific binding sites on the complete lac operator, and for the presence of a second, nonspecific binding site with lower affinity. The influence of this second site on the interpretation of the binding data is discussed.     

31P NMR study of daunorubicin-d(CGTACG) complex in solution. Evidence of the intercalation sites

The interaction of daunorubicin with the self-complementary DNA fragment d(CGTACG) was studied by 31P NMR spectroscopy. The individual phosphates have been assigned for the nucleotide and the complex and signals from bound and free species in slow exchange at 19 degrees C were detected. In solution, the hexanucleotide binds two molecules of daunorubicin, which intercalate in the d(CG) sequence at both ends of the helix. Evidence for local deformations of the backbone at the sites of C5pG6, C1pG2 and G2pT3 phosphates is given. The binding constants for the stepwise equilibrium and the rate of dissociation of the intercalated duplex were also determined.     

Two-dimensional 1H, 15N NMR investigation of uniformly 15N-labeled lac repressor headpiece.

15N uniformly labeled lac repressor and lac repressor headpiece were prepared. 15N NMR spectra of lac repressor were shown resolution inadequate for detailed study while the data showed that the 15N labeled N-terminal part of the protein is quite suitable for this type of study allowing future investigation of the specific interaction of the lac repressor headpiece with the lac operator. We report here the total assignment of proton 1H and nitrogen 15NH backbone resonances of this headpiece in the free state. Assignments of the 15N resonances of the protein were obtained in a sequential manner using heteronuclear multiple quantum coherence (HMQC), relayed HMQC nuclear Overhauser and relayed HMQC-HOHAHA spectroscopy. More than 80 per cent of residues were assigned by their 15NH(i)-N1H(i + 1) and 15NH(i)-N1H(i - 1) connectivities. Values of the 3JNH alpha splitting for 39 of the 51 residues of the headpiece were extracted from HMQC and HMQC-J. The observed 15NH(i)-C beta H cross peaks and the 3JNH alpha coupling constants values are in agreement with the three alpha-helices previously described [Zuiderweg, E.R.P., Scheek, R.M., Boelens, R., van Gunsteren, W.F. and Kaptein, R., Biochimie 67, 707 (1985)]. The 3JNH alpha coupling constants can be now used for a more confident determination of the lac repressor headpiece. From these values it is shown that the geometry of the ends of the second and third alpha-helices exhibit deviation from the canonical alpha-helix structure. On the basis of NOEs and 3JNH alpha values, the geometry of the turn of the helix-turn-helix motif is discussed.     

NMR study of the interaction between the lac repressor and the lac operator

Binding of the lac repressor headpiece, the N-terminal region of the lac repressor, to the lac operator of Escherichia coli was studied by 1H-NMR spectroscopy. Two DNA fragments, of 51 base pairs and 62 base pairs, containing the lac operator region, were investigated. The signals of their hydrogen-bonded imino protons were well resolved in the 500-MHz NMR spectra. The spectra of the free lac operator DNA are similar to those obtained from ring-current-shift calculations for a B-DNA structure. Complex formation with the headpiece led to small but nevertheless characteristic changes in the spectra. The fact that very few imino resonances shifted upon addition of headpiece, as well as the variety in direction and size of these chemical shifts, indicate the formation of a specific complex between the lac repressor and the lac operator. The observed changes in the resonance positions exclude the intercalation of tyrosine residues of the headpiece between adjacent base pairs of the lac operator as well as the formation of a cruciform structure. They rather reflect a small conformational transition in the DNA itself, caused for example by an alteration in the tilt of a few base pairs or a shift of the keto-enol tautomeric equilibrium of the bases towards the enolic form.     

Two-dimensional 1H and 31P NMR spectra of a decamer oligodeoxyribonucleotide duplex and a quinoxaline ((MeCys3, MeCys7]TANDEM) drug duplex complex

Assignment of the 1H and 31P NMR spectra of a decamer oligodeoxyribonucleotide duplex, d(CCCGATCGGG), and its quinoxaline ((MeCys3, MeCys7]TANDEM) drug duplex complex has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. The 31P chemical shifts of this 10 base pair oligonucleotide follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. While the 31P chemical shifts show sequence-specific variations, they also do not generally follow the Calladine "rules" previously demonstrated. 31P NMR also provides a convenient monitor of the phosphate ester backbone conformational changes upon binding of the drug to the duplex. Although the quinoxaline drug, [MeCys3, MeCys7]TANDEM, is generally expected to bind to duplex DNA by bis-intercalation, only small 31P chemical shift changes are observed upon binding the drug to duplex d(CCCGATCGGG). Additionally, only small perturbations in the 1H NMR and UV spectra are observed upon binding the drug to the decamer, although association of the drug stabilizes the duplex form relative to the other states. These results are consistent with a non-intercalative mode of association of the drug. Modeling and molecular mechanics energy minimization demonstrate that a novel structure in which the two quinoxaline rings of the drug binds in the minor groove of the duplex is possible.     

H NMR study of a complex between the lac repressor headpiece and a 22 base pair symmetric lac operator

A complex between the lac repressor headpiece and a fully symmetric tight-binding 22 bp lac operator was studied by 2D NMR. Several 2D NOE spectra were recorded for the complex in both H2O and 2H2O. Many NOE cross-peaks between the headpiece and DNA could be identified, and changes in the chemical shift of the DNA protons upon complex formation were analyzed. Comparison of these data with those obtained for a complex between the headpiece and a 14 bp half-operator, studied previously [Boelens, R., Scheek, R. M., Lamerichs, R. M. J. N., de Vlieg, J., van Boom, J. H., & Kaptein, R. (1987) in DNA-ligand interactions (Guschlbauer, W., & Saenger, W., Eds.) pp 191-215, Plenum, New York], shows that two headpieces form a specific complex with the 22 bp lac operator in which each headpiece binds in the same way as found for the 14 bp complex. The orientation of the recognition helix in the major groove of DNA in these complexes is opposite with respect to the dyad axis to that found for other repressors.     

Assignment of phosphorus-31 and nonexchangeable proton resonances in a symmetrical 14 base pair lac pseudooperator DNA fragment

The 31P chemical shifts of all 13 phosphates and the chemical shifts of nearly all of the non-exchangeable protons of a symmetrical 14 base pair lac pseudooperator DNA fragment have been assigned by regiospecific labeling with oxygen-17 and two-dimensional NMR techniques. At 22 degrees C, 8 of the 13 phosphorus resonances can distinctly be resolved while the remaining 5 resonances occur in two separate overlapping regions. The 31P chemical shifts of this particular 14 base pair oligonucleotide do not follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence the more upfield the 31P resonance occurs, as shown from other 31P assignment studies. Failure of this general rule is believed to be a result of helical distortions that occur along the oligonucleotide double helix, on the basis of the analysis of Callidine [Callidine, C.R. (1982) J. Mol. Biol. 161, 343-352]. Notable exceptions to the phosphate position relationship are 5'-Py-Pu-3' dinucleotide sequences, which resonate at a lower field strength than expected in agreement with similar results as reported by Ott and Eckstein [Ott, J., & Eckstein, F. (1985) Biochemistry 24, 253]. A reasonable correlation exists between 31P chemical shifts values of the 14-mer and the helical twist sum function of Calladine. The most unusual 31P resonance occurs most upfield in the 31P spectrum, which has been assigned to the second phosphate position (5'-GpT-3') from the 5' end. This unusual chemical shift may be the result of the predicted large helical twist angle that occurs at this position in the 14-mer sequence. Further, it is believed that the large helical twist represents a unique structural feature responsible for optimum binding contact between lac repressor protein and this 14-mer lac pseudooperator segment. Assignments of proton resonances were made from two-dimensional 1H-1H nuclear Overhauser effect (NOESY) connectivities in a sequential manner applicable to right-handed B-DNA, in conjunction with two-dimensional homonuclear and heteronuclear J-correlated spectroscopies (1H-1H COSY and 31P-1H HETCOR). Most nonexchangeable base proton and deoxyribose proton (except for some unresolved H4', H5', and H5" protons) resonances were assigned.     

Cooperative alpha-helix unfolding in a protein-DNA complex from hydrogen-deuterium exchange

We present experimental evidence for a cooperative unfolding transition of an alpha-helix in the lac repressor headpiece bound to a symmetric variant of the lac operator, as inferred from hydrogen-deuterium (H-D) exchange experiments monitored by NMR spectroscopy. In the EX1 limit, observed exchange rates become pH-independent and exclusively sensitive to local structure fluctuations that expose the amide proton HN to exchange. Close to this regime, we measured decay rates of individual backbone HN signals in D2O, and of their mutual HN-HN NOE by time-resolved two-dimensional (2D) NMR experiments. The data revealed correlated exchange at the center of the lac headpiece recognition helix, Val20-Val23, and suggested that the correlation breaks down at Val24, at the C terminus of the helix. A lower degree of correlation was observed for the exchange of Val9 and Ala10 at the center of helix 1, while no correlation was observed for Val38 and Glu39 at the center of helix 3. We conclude that HN exchange in the recognition helix and, to some extent, in helix 1 is a cooperative event involving the unfolding of these helices, whereas the HN exchange in helix 3 is dominated by random local structure fluctuations.     

Alkyl phosphotriester modified oligodeoxyribonucleotides. VI. NMR and UV spectroscopic studies of ethyl phosphotriester (Et) modified Rp-Rp and Sp-Sp duplexes, (d[GGAA(Et)TTCC])2.

1H NMR chemical shift assignments for the title compounds were made for all but a few H5' and H5" signals using two-dimensional nuclear Overhauser effect (2D-NOE) data, which was also used for the first time to assign absolute configuration at phosphorus. The chemical shifts were, in general, similar to those reported [Broido, M.S., et al. (1985) Eur. J. Biochem. 150, 117-128] for the B-like conformation of the unmodified, parent duplex, [d(GGAATTCC)]2. Differences in chemical shifts for corresponding protons were mostly localized to the AA(Et)TT region, and showed some stereochemical dependence. Unambiguous assignment of the phosphotriester 31P signals was achieved in a novel way using selective insensitive nucleus enhancement by polarization transfer (selective INEPT) NMR. The Rp-Rp duplex melted ca. 11 degrees C lower than either the Sp-Sp or parent duplexes, as evidenced by Tm and variable temperature 1H/31P NMR measurements. The 2D-NOE data for the Rp-Rp duplex suggested possible steric interactions between the ethyl group and the H3' of the flanking A residue. At low ionic strength, the Sp-Sp and parent duplexes had similar stability but at high ionic strength the Sp-Sp duplex was less stable.     

Photo-CIDNP study of the interaction between lac repressor headpiece and lac operator DNA

Lac repressor headpiece (HP) and intact lac repressor have been studied using the photo-CIDNP method. At neutral pH histidine 29, tyrosines 7, 12 and 17 and methionine 1 are polarised. His-29 polarizations are weaker and broader in HP59 than in HP51 indicating that the C-terminal octapeptide in HP59 adopts a conformation that allows an interaction with His-29. The photo-CIDNP spectra of intact lac repressor and HP51 are very similar, showing that the same residues are accessible to the photo-excited flavin. An equimolar mixture of HP51 and a 14 base pair lac operator fragment strongly suppresses the photo-CIDNP effect of tyrosines 7 and 17 and abolishes the His-29 polarizations. The results are compared with earlier photo-CIDNP measurements on a complex of headpiece with poly[d(AT)] and with a model derived from a 2D NMR study on a lac headpiece-operator complex.     

Analysis of complex formation of daunomycin with the deoxyhexanucleotide 5'-d(TpApCpGpTpA) in an aqueous solution from NMR-spectroscopy data

Self-association of hexadeoxynucleotide 5'-d(TpApCpGpTpA) and its complexation with antitumor antibiotic daunomycin were studied by one- and two-dimensional homonuclear 1H NMR spectroscopy and heteronuclear 1H-31P NMR spectroscopy in water-salt solution. The concentration and temperature dependences of proton chemical shifts of the hexadeoxynucleotide and the ligand were measured, and equilibrium constants and thermodynamic parameters of corresponding reactions were calculated on this basis using models for the formation of hexadeoxynucleotide duplex and its complex with the antibiotic. The spatial structure of daunomycin-d(TACGTA)2 complex in solution was calculated using X-PLOR software on the basis of 2D NOE spectral data and the limit values of proton chemical shifts of the ligand. Comparative analysis of different intermolecular interactions in sequence-specific binding of the antibiotic to the DNA fragment was carried out.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Duplex and quadruplex DNA binding and photocleavage by trioxatriangulenium ion

The stable trioxatriangulenium ion (TOTA) has previously been shown to bind to and photooxidize duplex DNA, leading to cleavage at G residues, particularly 5'-GG-3' repeats. Telomeric DNA consists of G-rich sequences that may exist in either duplex or G-quadruplex forms. We have employed electrospray ionization mass spectrometry (ESI-MS) to investigate the interactions between TOTA and duplex DNA or G-quadruplex DNA. A variety of duplex decamer oligodeoxynucleotides form complexes with TOTA that can be detected by ESI-MS, and the stoichiometry and fragmentation patterns observed are commensurate with an intercalative binding mode. TOTA also forms complexes with four-stranded and hairpin-dimer G-quadruplex oligodeoxynucleotides that can be detected by ESI-MS. Both the stoichiometry and the fragmentation patterns observed by ESI-MS are different than those observed for G-tetrad end-stacking binding ligands. We have carried out (1)H NMR titrations of a four-stranded G-quadruplex in the presence of TOTA. Addition of up to 1 equiv of TOTA is accompanied by pronounced upfield shifts of the G-tetrad imino proton resonances in the NMR, which is similar to the effect observed for G-tetrad end-stacking ligands. At higher ratios of added TOTA, there is evidence for additional binding modes. Duplex DNA containing either human telomeric repeats (T(2)AG(3))(4) or the Tetrahymena telomeric repeats (T(2)G(4))(4) are readily photooxidized by TOTA, the major sites of oxidation being the central guanine residues in each telomeric repeat. These telomeric repeats were incorporated into duplex/quadruplex chimeras in which the repeats adopt a G-quadruplex structure. Analysis by denaturing polyacrylamide gel electrophoresis reveals significantly less TOTA photocleavage of these quadruplex telomeric repeats when compared to the duplex repeats.     

Two-dimensional 1H and 31P NMR spectra and restrained molecular dynamics structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex

Assignment of the 1H and 31P NMR spectra of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, has been made by two-dimensional 1H-1H and heteronuclear 31P-1H correlated spectroscopy. The downfield 31P resonance previously noted by Patel et al. (1982) has been assigned by both 17O labeling of the phosphate as well as a pure absorption phase constant-time heteronuclear 31P-1H correlated spectrum and has been associated with the phosphate on the 3' side of the extrahelical adenosine. JH3'-P coupling constants for each of the phosphates of the tridecamer were obtained from the 1H-31P J-resolved selective proton-flip 2D spectrum. By use of a modified Karplus relationship the C4-C3'-O3-P torsional angles (epsilon) were obtained. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. The 31P chemical shifts and epsilon torsional angles follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. Because the extrahelical adenosine significantly distorts the deoxyribose phosphate backbone conformation even several bases distant from the extrahelical adenosine, 31P chemical shifts show complex site- and sequence-specific variations. Modeling and NOESY distance-restrained energy minimization and restrained molecular dynamics suggest that the extrahelical adenosine stacks into the duplex. However, a minor conformation is also observed in the 1H NMR, which could be associated with a structure in which the extrahelical adenosine loops out into solution.     

Interaction between the lac operator and the lac repressor headpiece: fluorescence and circular dichroism studies

The interaction between the lac repressor headpiece and a small operator DNA fragment has been examined by fluorescence and circular dichroism (c.d.) measurements. Binding of the headpiece to the DNA fragment induces a strong quenching of the fluorescence of its tyrosine residues. Quantitative analysis of the fluorescence data demonstrates that, in a first step, two headpieces bind very strongly to the DNA fragment then weaker binding occurs. C.d. demonstrates that the binding induces conformational changes of the DNA. The c.d. change produced upon binding of the first two headpieces differs from that induced upon binding of two further headpieces . Binding of the second pair of headpieces is similar to non-specific binding to non-operator DNA. The conformation of the operator DNA in the presence of two headpieces differs drastically from that in presence of lac repressor. Addition of the core to the lac operator does not induce any conformational change of the nucleic acids. These results are discussed with respect to the relative roles of core and headpieces in the lac repressor-lac operator interaction.     

lac repressor mutants with double or triple exchanges in the recognition helix bind specifically to lac operator variants with multiple exchanges.

Several lac repressor mutants have been isolated which repress beta-galactosidase synthesis in Escherichia coli up to 200-fold. They do so by binding specifically to particular symmetrical lac Oc operator variants. The mutations in the lac repressor are localized in two separate parts of the recognition helix comprising (i) residues 1 and 2 which interact with base pairs 4 and 5 of lac operator and (ii) residue 6 which recognizes operator base pair 6. Mutations of residues 1 and 2 may be combined with a mutation of residue 6. The resulting mutant protein binds specifically to an operator variant with three symmetric exchanges in base pairs 4, 5 and 6.     

Structure of the complex between lac repressor headpiece and operator DNA from measurements of the orientation relaxation and the electric dichroism.

The complex between lac repressor headpiece and short rodlike DNA fragments containing the lac operator sequence is characterised by measurements of the rotation diffusion. Using the method of electric dichroism we measure the rotation relaxation and determine changes in the length of the DNA upon ligand binding with high accuracy. According to these measurements any change in the length of the operator DNA upon binding of the first two headpiece molecules remains below 1A; the electric dichroism also remains virtually unchanged. At high degrees of (unspecific) binding we observe an increase in the rotation relaxation time, which is attributed to an increase of the apparent mean radius of the complex. As a control of our procedure for the determination of length changes we use the intercalation of ethidium bromide and arrive at an increase of the DNA length per bound ethidium of 3.2A (at 3.4A rise per base pair). The results obtained for the headpiece operator complex are not consistent with models assuming large changes of the DNA structure or intercalation of tyrosine residues.     

Effect of distortions in the deoxyribose phosphate backbone conformation of duplex oligodeoxyribonucleotide dodecamers containing GT, GG, GA, AC, and GU base-pair mismatches on 31P NMR spectra

We have previously suggested that variations in the 31P chemical shifts of individual phosphates in duplex oligonucleotides are attributable to torsional angle changes in the deoxyribose phosphate backbone. This hypothesis is not directly supported by analysis of the 1H/31P two-dimensional J-resolved spectra of a number of mismatch dodecamer oligonucleotide duplexes including the following sequences: d-(CGTGAATTCGCG), d(CGUGAATTCGCG), d(CGGGAATTCGCG), d(CGAGAATTCGCG), and d(CGCGAATTCACG). The 31P NMR signals of the dodecamer mismatch duplexes were assigned by 2D 1H/31P pure absorption phase constant time (PAC) heteronuclear correlation spectra. From the assigned H3' and H4' signals, the 31P signals of the base-pair mismatch dodecamers were identified. JH3'-P coupling constants for each of the phosphates of the dodecamers were obtained from 1H/31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. JH3'-P coupling constants were measured for many of the oligonucleotides as a function of temperature. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. This correlation can be further extended to the C3'-O3'-P-O5' torsional angle (zeta) by using a linear relationship between epsilon and zeta obtained from crystal structure studies. The 31P chemical shifts follow the general observation that the more internally the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. In addition, 31P chemical shifts show sequence- and site-specific variations. Analysis of the backbone torsional angle variations from the coupling constant analysis has provided additional information regarding the origin of these variations in 31P chemical shifts.