DNA, The versatile vector of life: two-dimensional NMR studies

DNA, The versatile vector of life: two-dimensional NMR studies

109 Journal ofMolecuhr Structure, 141(1986) 109-126 Elsevier Science publishers B.V., Amsterdam - printed in The Netherlands DNA, THE VERSATILE VECT...

1000KB Sizes 1 Downloads 158 Views

109

Journal ofMolecuhr Structure, 141(1986) 109-126 Elsevier Science publishers B.V., Amsterdam - printed in The Netherlands

DNA, THE VERSATILE VECTOR OF LIFE:

TWO-DIMENSIONALNMR STUDIES

C. ALTONA Gorlaeus

Laboratories,

State

University,

P.O.Box

9502,

2300 RA Leiden

(The Netherlands)

ABSTRACT Deoxyribonucleic acid (DNA), until recently regarded as a relatively stiff in reality displays a multitude of surprising and regularly built double helix, local structural variations. According to X-ray crystallographic findings, three different families of DNA structures exist: the right-handed duplexes A DNA and Modern 2D NMFt techniques now alB DNA as well as the left-handed duplex Z DNA. low for unequivocal assignment of base proton resonances and at least Hl’, H2’, Nuclear Overhauser Enhancement (NOE) intenH2” signals in intact DNA duplexes. sities then lead to a quick determination of the overall structure (B or Z DNA, the A form thus far has not been detected in aqueous solution). However, the study of finer structural details requires the determination and interpretation Examples of local structural variations of the of vicinal coupling constants. sugar ring in single-helical as well as in double-helical DNAs are given.

INTRODUCTION A DNA molecule backbone;

any of

tached

the

to

the

sugars

(or

In the

course

translated access sary

into

only for

of

to

background

its

a protein a small

of

own function.

is

off

the

genome at the cular to

biology:

recognize

proper

of

time.

what are faithfully

0022-2660166/$03.50

the

the

the

time All

refs

a cell

the

a specific

above location

two complementary

copied.

requires

call

that

the the

the

cell

the

on DNA in order

-

of

e.g., DNA

on the

in mole-

allow

block

neces-

of

place

problem that

to

has

part

expression

proper

a central level

0 1966 Elaevier Science Publishers B.V.

specialized

and finally

upon synthesis

at

illustrates at atomic

In the

namely,

means that

and must be activated

features

may pair

a double-helical

must be read

information, would

by

1-5.

faithfully

This

in-

characterized

to

contains

information

total

needed.

rise

at-

pyrimidines

equilibrium

is

two DNA molecules

Specialization

of

is

see

DNA is

coded

Moreover,

when it

most

the

which

C, giving

DNA polymer

are

A and G and the

stack

In addition,

sequence.

enzyme only

purines

diester

thereof)

in a conformational

information

fraction

a specific blocked

exists

T and G with

reproduction parts

the

a duplexed

sugar-phosphate

some modifications

single-helical

overlap.

Thus,

specific

cell,

solution)

A with

structure.

a deoxyribose (and

position,

and the

bonds:

For general

strands.

(in

of

bases

Cl’

base-base

hydrogen

duplex)

living

at the

random coil

nearest-neighbour through

1) consists different

Such a chain

C and T. volving

(Fig. four

it

a protein or

to

110

Fig, 1. Arbitrary sequence of a self-complementary deoxyribotetranucleoside triphosphate: d(TGCA). The standard numbering of base and sugar entities is shown. The residues themselves are numbered (l), (2) ... starting with the St-OH terminal (T(1) in this particular example). activate it. Most of these features are poorly understood and much further work remains necessary. Many questions concerning the structure and dynamics of protein-DNA complexes and of protein-free DNA are presently under investigation in many laboratories. Two powerful techniques, besides others, are currently applied in studies on DNA and DNA complexes:

(a) X-ray analysis of single crystals;

(b) nuclear mag-

netic resonance (NHR). The purpose of the present paper is twofold. First, a brief review of modern aspects of the conformational analysis of the deoxyribose ring in DNA constituents -

both single-helical as well as double-helical -

will be given, together with selected NMR techniques that have proven their worth. Emphasis will be laid on the use and interpretation of NMR coupling constants. Secondly, our results concerning B DNA will be contrasted to findings by other researchers.

CONFORMATIONAL ANALYSIS OF DEOXYRIBOSES In biochemical literature the two distinct B-D-deoxyribose conformations are commonly denoted as C3'-endo (3E) and C2'-endo (zE). Note that these terms tend to give the false impression that these sugar conformations are rigid, whereas in reality the pseudorotational movement

161 of the furanose ring constitutes

one of the major degrees of freedom of the backbone of B DNA. For this reason we prefer to distinguish between the N-type genus of forms (P = 0' f 90') and the S-type genus (P = 180' * 900), Fig. 2. The endocyclic torsion angles +j in an

111

H2

BASE

N-conformer

S-conformer

Fig. 2. N- and S-type conformations of furanose rings. Idealized forms are shown, whereas in fact N and S stand for large ranges of possible conformations [6,7] . For the sake of clarity, the HZ" proton has been omitted. The classical A-DNA form is characterized by an N-N-N-N sequence of sugars, the classical B-DNA by an S-S-S-S sequence, whereas the Z-DNA structure displays alternation: S-N-S-N .

-

approximately -

rameter equation

oj

=

equilateral S-membered ring are mutually related by a 2-pa-

[6] , es.(l).

$m cos[P * 4n(j-2) / 51

j = 0, 1, 2, 3, 4

in which 4, represents the maximum amplitude of pucker, P the phase angle and the endocyclic torsions are numbered clockwise, starting with C4'-04'-Cl'-C2' (j = 0). Knowledge of two or three endocyclic torsions, e.g., from X-ray crystallography or from proton NMR coupling constants, enables one to uniquely deduce the pseudorotation parameters P and $,, from which not only the remaining torsions but also the important backbone angle 6(C5'-C4'-C3'-03')can be calculated, eq. (2).

6

=

120.6 +

1.1 $m cos(P + 145.2)

Eq. (2) is based upon a corrected pseudorotation equation

(2) [S]

. The constants

of eq. (2) were determined from a least-squares fit to our crystal structure data set

[9]

.

Given the fact that in aqueous solution of nucleic acids the N- and S-type conformers most often occur side-by-side in fast equilibrium, the NMR spectroscopist is faced with the simultaneous determination of four geometrical parameters (PN, $,, Ps, 4~) and one equilibrium constant K. The Leiden group has spent considerable effort in order to provide a dependable and routine solution to this problem. Two keys to success were developed. First, de Leeuw et at. [9] derived a set of empirical equations which relate proton-proton torsion angles to the endocyclic furanose torsions. Eq. (3) applies to 8-D-deoxyribose:

112 =

121.4O

+

1.03 $m cos(P

- 144O)

Ua)

$1’2”

=

0.9O

+

1.02 +m cos(P

- 144O)

WI

$2,3,

=

2.40

+

1.06@mcosP

(3c)

=

122.9O

+

1.06$mcosP

(3d)

+

1.09 0, cos(P

%,2,

‘2”3

= -124.0”

$3’4’

Secondly,

Haasnoot et ~2.

equation

by the introduction

and orientation [4,

iterative

i.e.

acids

not be treated

couplings

curacy of the procedure

is greatly

displays

with temperature,

fact,

shift

the geometry of the sugars,

accuracy

here.

acid

fragments

that is at least

Karplus

of these findings it to say that an

yields

the pseudorotation

sugars,

and the N/S

along the bonds Cl’-C2’

,

are shown in Table I. The ac-

enhanced in case the N/S equilibrium as is usually often

constant

found for oligomers.

and thus the important in solution‘now

to

has been documented in

Suffice

set of 35HH couplings

. Some representative

of small nucleic

application

and prolines

geometry of the N- and S-type

given a complete

a strong

three-parameter

for the electronegativity

computer program PSEUROT [12]

the detailed

, C3’-C4’

The successful

of nucleic

and will

least-squares

molar ratio, C2’-C3’

analysis

11, 121

parameters,

extended the classical of terms which account

of substituents.

the conformational estenso

[lo]

+ 144O)

In

backbone angle 6 (Eq.2),

can be determined

equal to that of a high-resolution

with an

X-ray crystal

struc-

ture analysis. The determination the determination COSY spectra. ed

-

of the sugar geometry in intact

of signal

These signal

from the foregoing

pseudorotational each residue, prerequisite haviour

widths of Hl’ width contain equations

-

for each conformer

before

valid

from at least conclusions

can be drawn. In addition,

31-phosphorus

, H2’ and H2” from either sums of couplings

pathway, Table I and Fig. information

decoupled

spectra

three

concerning

analysis

(N/S) ratio

formational

species

with two conformational

P(S)

=

(31.7

-

(ii)

and

observables

(PS in B DNA). Clearly,

the population

ratio

C2”) / 10.7

(Hl’

or

for

, H2’ , H2”) is

the sugar conformational

two main degrees

be-

width in This

of freedom:

the phase angle P of the major conself-consistency

parameters

is best

along the

that,

should be attempted whenever possible.

(i)

upon

1D or from

encountered

of the H3’ resonance

that one has to determine

practice,

resonances

rests

which can be predict-

3. It should be stressed

stems from the fact the population

double helices

approximated p(N)

of three or four

lends credence

=

from

22”,

(22” - 21.0)

to the results. eq.

In

(4):

/ 10.7

(4)

113

TABLE I Calculateda

coupling

formations

constants

(Hz) of deoxyribose

running from N-type via Ol’-endo

principle

obtainable

from COSY spectra,

rings

to S-typeb.

are also

in DNAfor

sugar con-

sums of couplings,

in

presented.

Deoxy N

21T

ZT ZT !T

2E

3E

qE

OE

P

6c

1’2’

1’2”

2’3’

2”3’

3’4’

324 342 0 18 36 54 72 90

108 97 89 84 82 84 90 99

1.3 1.3 1.5 2.1 3.3 5.0 6.8 8.4

6.9 7.1 7.7 8.4 8.9 9.0 8.4 7.5

7.5 7.3 7.2 7.5 8.1 9.0 9.6 9.6

8.7 9.5 9.7 9.5 8.7 7.3 5.5 3.6

5.3 6.7 7.7 8.3 8.4 8.3 7.7 6.7

108 126 144 162 180 198 216

110 121 133 144 152 158 159

9.5 10.0 10.2 10.0 9.5 8.4 6.8

6.5 5.9 5.6 5.7 5.8 6.2 6.6

8.9 7.8 6.6 5.8 5.5 5.6 5.9

2.1 1.3 1.2 1.2 1.3 1.2 1.2

5.3 3.7 2.3 1.4 1.0 0.9 0.9

d

C2fe

12”f

E3’g

8.2 8.4 9.2 10.5 12.2 14.0 15.2 15.9

22.8 22.6 22.7 23.6 25.4 28.0 30.4 32.0

29.6 30.6 31.4 32.0 31.6 30.3 27.9 25.1

21.5 23.5 24.6 25.3 25.2 24.6 22.8 19.9

16.0 15.9 15.8 15.7 15.3 14.6 13.4

32.4 31.8 30.8 29.8 29.0 28.0 26.7

22.6 21.2 20.8 20.9 21.1 21.4 21.8

16.3 12.8 10.1 a.4 7.8 7.7 7.7

Zl’

Deoxy S !T TT

lE 2E

ZT 3B ZT

aThe calculations were carried out by use of the extended Karplus equation [lo] and eq. (3); the initial values were corrected for the Barfield through-space btransmission effect [13]. A mean puckering amplitude Ornof 35” was assumed[4]. See refs 6,7 for the definition of the phase angle of pseudorotation P and the relation between P and the various twist (T) and envelope (E) conformations. iThe backbone angle 6 (C5’-C4’-C3’-03’) was calculated from eq. (2). El’ = 51121 + Jl’2” , the line width of the Hl’ signal. eZ2’ = 51121 + J2~3t + ,J212” . An average value of 14:DHz was used for J2’211 . 52” = 51~2~~+ J2”3 + 5212” ; see footnote e. gz3 ’ = J213t + ~J2~~31+ 53’41 , the line width of the H3’ signal in 31-phosphorus decoupled spectra. Eq. (4) is valid Table VIII)

521211 coupling is known cation

for PN O-18’,

a similar

PS 126-l&?O’, and brn 35”.

equation

was presented

has been subtracted

for a given residue,

of X2’ gives

from eq.

Recently, ed [13] -

attempts sight

these distances was claimed

a,

Once the S/N ratio of the lo-

pathway. Alternatively,

by means of program PSEUROT[12]

. Once

species

is known, backbone angle 6

sugar puckerings

from NOEShave been publish-

(2). to determine

. Since NOE intensities

at first

4 (footnote

a good impression

on the pseudorotational

can be done automatically

the phase angle of the major conformational follows

In ref.

however, the geminal

from the sums of couplings.

analysis

of the sugar conformation

the calculation

in which,

-

(iteratively)

[13]

depend upon l/r6

be tempted to translate to derive

that distances

(r = H . . . H distance)

the DNA sugar conformation.

in an intact

one could

NOESinto H . . . H distances duplex could

be fitted

and use

In fact,

it

with an

114 r.m.s.

accuracy

of 0.2 A.

tered at Cl’-exo

(‘E,

P 126”).

cism,

however.

falls

such as different

ization

according

and local

protons

[14] ).

Two pitfalls

B DNAduplexes

second,

usually

analyses exhibit

remain constant

within

and that concerns

H . . . H distances

to Cl’-exo

the stated

the Hl’

(&de

(126”),

freedom

insensitive

H . . . H distances

going from the center

of the

occurs

distances,

to N-type sugar.

cannot be detected

the outcome of the iterative

all

toas

generally

One exception

(c 0.35 A). Moreover,

except HZ” . . . H4’ , do not change much on going from S-type it might influence

completely

deoxyribose

are quite

limit.

pit-

then loom up:

the H . . . H distances

[13] 0.2 A error

. . . H4’ distance

origin infra),

conformational

Fig.. 4 shows the endocyclic

This means that a 10 - 20% N-type population although

critiof magnet-

can be avoided

residual

of the phase angle P. It can be seen that,

normal Ps range (150 - 160”)

motions and relay

protons

of conformational

constant

most intrasugar

ward changes of sugar pucker.

molecular

or via solvent

to our coupling

sugars in intact

function

thus found cen-

Such an approach remains open to serious

overall

U&J neighbouring

(pN 6 20%):

of the sugar puckerings

Let us, for the sake of argument, grant that major technical

(which is doubtful first,

The majority

by a NOEanalysis,

structure

calculation.

Fig. 3. Predicted sums of couplings (see footnotes d, e, f of Table I) along the pseudorotational pathway of the deoxyribose ring. In this example Ornwas kept constant at 35”.

115

u.00

3.75

2

p.50

5

5.25

P

E3.00 '" 0 I t 2.75 x

2.50

2.25

Fig. 4. Endocyclic H . . . H distances in deoxyribose as function of the phase angle of pseudorotation P at constant puckering amplitude (4, 35’). (Courtesy of Dr J.-R. Mellema). For example, distance)

the presence

of a small percentage

would increase

terpreted

-

the observed

would result

H2”/H4’

in a calculated

by a NQE study

distance

of 0.29 - 0.31 A is found for all

Fig.

whereas it is stated

4),

(S-type

[15] of double-stranded

sugars).

H2”...H4’

improperly

in-

d(CGTACG). An apparent H2”...H4’

residues

[15] that this

The latter

- if

PS that is too low. A good example is

provided

structure

of N-type sugar (short NOEand this

(i.e.,

in the N-type range,

duplex adopts a conventional

implies

a H2”...H4’

distance

B-type

of > 3.8 h .

ASSIGNMENT OF RESONANCES The principles recently

by Bax [16].

information parameters interest

of two-dimensional In most

-but

(2D) NMRspectroscopy not all

is added, compared to classical into

two different

for NMRof larger

frequency

molecules

and Hl’,

in a duplex decamer [la],

dimensions powerful

H2’, H2” resonances

current

types of 2D experiments

1D NMR, but the separation

assignment

no new of NMR

F1 and F, is of particular

with many overlapping

the 1D NOEmethod was shown to be quite ments of base proton

-

have been summarized

resonances.

and actually

in a duplex octamer

strategies

Although

led to full

for nucleic

assign-

[17] and

acids

favour

116

0

INTRARESIDUE SCALAR

COUPLINGS

KEY

NOES

(COSY)

ACTORS

(NOESYI

H5/CH3-H6 Hl’l

Fig. 5. Intraresidue coupling networks which are revealed by COSY and intraresiwhich give rise to NOESYcross peaks (for the sake of due short H... H distances clarity some NOESYconnectivities have been omitted). the simultaneous 19-241. lar

study of both COSY and NOESYspectra

In brief,

cross

(J) couplings

the existence nectivities

peaks in a COSY contour

between protons,

The normal COSY experiment strand

or duplex)

array.

In practice,

volving nicity

of e.g.

peaks.

often

Problems

ly introduced [20].

Fig.

experiment,

causes (ii)

and (iii)

extension

20

H4’ resonances

(iii)

HOMONUCLERR

7 gives

5.

in single

cross

accidental

peaks inisochro-

the complete

as-

in the H4’,

the envelope

HS’,

of diagonal

by the employment of a recent-

termed homonuclear

double-RELAY

schemes of a normal COSY and of the double-relayed

dubbed DAYCOSYby the present Fig.

con-

sugar in the

of signals

into

can be circumvented

(i)

(ii)

prohibits

overcrowding

of the COSY experiment,

6 shows the pulse

network.

the noise;

of sca-

reflect

are shown in Fig.

to each individual

peaks to disappear

mixing time unambiguous subspectra plete)

into

[S,

The intraresidue

may be encountered:

disappear

involved; cross

the existence

(n = number of residues

that belong

difficulties

often

(< Q, 5 A).

the n sets

two or more H3’ and/or

signment of the networks H5” region

yields

several

distances

betray

thereof

peaks in a NOESYplot

from COSY and NOESYspectra

of seven signals

small couplings

plot

whereas cross

of short proton-proton that one may expect

or variants

author.

are obtained

By a judicious for

each complete

choice

of the

(or almost com-

an example of a DAYCOSYspectrum for the simple

CORRELRTION

SPECTROScOPY

Fig. 6. Pulse scheme of a normal COSY (left) ment [ZO].

20

HOMONUCLEAR COHERENCE

DOUBLE-RELAYED SPECTROSCOPY

and of a DAYCOSY(right) experi-

117

0

wm

Fig. 7. High-field region of the 300 MHz DAYCOSY spectrum of d(TA). Extra cross peaks [compared to a standard COSY experiment) are encircled in the lower righthand half of the contour plot only. Chemical shifts are referenced to the methyl peak of tetramethylammonium chloride (TMA, 3.18 ppm upfield from DSS or TMS). dimer d(TA).

The “extra”

cross

peaks (absent

is seen that HI’ to H4’ can be correlated region

offers

no particular

A similar

problems.

Drobny [ZO] on the trimer AZ’-5*A2’-5’A. tell

us which set of signals

e.g.

refs

dure rests

belongs

to which residue. [19].

was obtained

22-24)

upon the fact

located

on residue own residue

and will

that a base proton

n in a right-handed (Fig.

5) and Hl’,

of the 5’ terminal

spectroscopist

first

not be treated

Therefore,

stack,

here in

(H8 of purines

of the strand).

does not

the information

of NOESYin the matter

or H6 of pyrimidines),

n-l

It is for this

(See

The proce-

extenso.

is sandwiched between Hl’,

HZ’, HZ” of residue

turns his attention

by Bax and

experiment

The application

It

that the HS’, H5”

of A- and B-DNA fragments has become an almost routine 5, 13-15,

of its

direction

result

and also

A COSY or DAYCOSY

must be supplemented by a NOESYspectrum assignment

in a normal COSY) are encircled.

at a glance

“above”

(i.e.

HZ’, HZ” in the

reason that the NMR

to the NOEconnectivities

between the

118 residue n-l

5'

-

n -

n+l

3'

H8/H6 ZI=L. Hl' - HZ' - HZ" t----,

H8/H6

--Zc__-W--, Hl' - HZ' - HZ" a

-

_??-Hl' - HZ' - HZ" 1_ Scheme 1. NOESY cross peaks most used in assignment of A- and B-DNA lH NMR spectra. H8/H6 protons and the Hl', HZ', HZ" regions. The procedure is outlined in Scheme 1. It shows that, except at the 5' terminal, each H8 or H6 should display two connectivities to Hl' resonances, one intra-, the other interresidue. An independent check is possible through the H8/H6 - HZ'/H2" and Hl' - HZ'/HZ" NOESY cross peaks. It should be mentioned here that, loosely speaking, the NOE is a function of l/r6 (r = H...H distance), of the motional correlation time rc' and of the mixing time used. Near the condition wrc = 1.1, where w is the angular Larmor frequency, the classical (i.e. laboratory frame) NOE vanishes. Thus, for a slowly tumbling nucleic acid strand (.ccg O.S~lp-~s) it is often difficult to observe NOES on NMR spectrometers operating in the 200-300 MHz range. An interesting way out of this difficulty appears to be offered by the ROESY experiment, recently advocated by Bax and Davis [Zl]. The rotating-frame NOE (under spinlocked conditions) is always positive and monotonically increases for increasing values of 'c C’

CONFORMATIONAL ANALYSIS OF DNA FRAGMENTS Single helices A single-helical DNA fragment is stabilized solely by "vertical" base-base overlap interactions (stacking) [l-3]. In contrast, the DNA duplex, besides stacking, has an additional stabilization through the formation of hydrogen bonds between bases on opposite strands, as well as the possibility of interstrand steric and/or charge-charge interactions. In order to understand the structure and dynamics of DNA, it is therefore important to study the behaviour of single helices of varied sequence in conjunction with the behaviour of the corresponding duplexes. We have set out to answer the following questions: (i) what conformational changes, if any, occur in going from the single helix to the double helix? (ii) do single-helical and duplex regions display conformational purity of the sugar rings or is N/S ring flip still possible:

119 TABLE II Pseudorotation

parameters

and predicted

6 deviations

mers,

see

[6,7],

backbone angle 6 (CS’-C4’-C3’-03’)

from mean (141’)

of various

single-helical

and observed DNAoligo-

text.

Compound

Residue

d (TAAT)1

T(1) A(2) A(3) T(4)

160 167 153 143

A(2) T(3) A(4)

T(1)

pS

Temp. range (“C)

6S(o)

A60bs

*‘pred

35 36 40 32

143 147 141 132

+2 +6 0 -9

0 +5 0 -5

8-43 I, II II

158 160 149 166

35 36 35 32

142 143 136 144

+l +2 -5 +3

0 +5 -5 0

a-44 ,I II ,t

A(1) T(2) A(3) T(4)

173 150 170 148

36 36 36 31

150 138 149 134

+9 -3 +8 -7

+5 -5 +5 -5

8-35 II

d (‘ITA) 3

T(1) T(2) A(3)

155 155 162

34 34 34

140 140 143

-1 -1 +2

0 0 0

5-65 II II

d&W4

A(1) A(2) A(3)

170 157 153

36 36 34

149 142 138

+8 +l -3

+5 0 -5

11-25 ,1 11

d(AA)4

A(1) A(2)

167 153

36 35

147 139

+6 -3

+5 -5

o-25 II

d (TGTG)5

T(1) G(2) T(3) G(4)

150 158 149 166

37 37 35 33

138 142 137 144

-3 +l -4 +3

0 +5 -5 0

27 11 ,I 11

d (CAACTT)6

C(1) A(2) A(3) C(4) T(5) T(6)

158 171 153 162 159 -

34 36 38 35 35 -

141 149 140 144 142 -

0 +8 -1 +3 +l -

0 +5 0 -5 0 0

15-35 II 11 II 11 11

d (CGT)7

C(l) G(2) T(3)

154 163 146

35 37 33

140 145 134

-1 +4 -7

0 +5 -5

2-52 11 II

d(TCG)7

T(l1 C(2) G(3)

158 157 160

33 35 34

141 141 142

0 0 +l

0 0 0

12-52 11 II

d (TATA)1

d (ATAT]2

8 35

l) J.-R.Mellema, A.K.Jellema, C.A.G.Haasnoot , .J.H. van Boom, and C.Altona, Eur. J.Biochem. 141 (1984) 165-175. 2, Refs 25.26. 3, C.A.G.Haasnoot, J.-R.Mellema, and C.Altoncunpublished results. 4, C.S.M.Olsthoorn, L.J.Bostelaar, J.H. van Boom, and C.Altona, Eur. J. Biochem. 112 (1980) 95-110. 5, L.J.Rinkel, J.-R.Mellema, and C.Altona, to be pubzhed. 6, L.P.M.Orbons and C.Altona, to be published. 7, J.-R. Mellema, R. van der Noerd, G.A. van der Marel, J.H. van Boom, and C.Altona, Nucleic Acids Res. -12 (1984) 5061-5078.

120 Crystallographic

information

concerning

single-helical

and knowledge must be gained from NMRspectroscopy. ber of such fragments have been studied ings will

be summarized here.

and 6S of deoxyribose vicinal

coupling

sugar ring

constants,

by PS ‘L 150’. character

Closer

(purine

of roughly

sugar conformational particular

occur

variation

Cg proposed

simple additivity sequence.

error

of such a step, [27],

especially

Elellema et al.

the values

values

negative

concerned

Analogous

[25,

appear

values to the

261 devised

6 variation

a

with DNA

the values

+1, -1,

0, 0. The unit of sum function

Fig.

g6

8 shows the A60bs (= 6obs-1410)

61 residues.

for the observed

already

noted

concomitant

results

in an in-plane

similar

sliding

effect,

in the single

in a B-DNA single-helix

An overall

with an increase sliding

leading helix.

r.m.s.

Fig.

deviation expe-

to 133” slides

less

aligned

with the adenine C6-N6 dipole.

tion

of electrostatic of 6

has another

the purine

interesting

a change in 6 leads

of the Cl’-N

vectors

this

where the

and becomes more or

suggests

a.B,Y,e,<)

to an appreciable residues,

a dA-dT step

that optimiza-

role.

conformational

of successive

is

6 angle from the

a position

6-membered ring This result

for

6

that a

interaction,

principle

the T base into

(angles

angle

base-paired

It is proposed

stacking

of the thymidine

may play a certain

pathway of the phosphate backbone fixed,

of 6 of the opposite

9 illustrates Reduction

by computer

of a pyrimidine

of the base pair.

is positioned,below interaction

is suggested

to optimum base-base

conformation.

value in the monomer (152’)

6 variation

[27] that a decrease

C4-04 dipole

more or less

changes of

when the residue

us to correlate

angle difference.

purine,

entation

II shows that positive

in a dR-dR or dR-dY step;

incorporating

explanation

Dickerson

Variation

the important

on local

which is of the same order of magnitude as the estimated

in a B-DNA duplex,

general

into

Since the

in 6.

A possible modelling.

operative

to the sugar in ques-

now be centered

in Table

by Dickerson

at 5O of torsion is found,

in

to the

in the sequence.

of 6 is seen in dY-dR or dY-dY steps.

are assigned

of 16 sequences,

rimental

variability

is not related

The sugars at each dR-dR and dR-dY step are assigned

is fixed of 2.6’

bases

will

scheme I& which allows

dY-dR and dY-dY steps values

pattern

P and 4, can be translated

R residue

in the second residue

sum function

of

angle.

with the first

is dY. No large

PS, 0,

analysis

in the case of d(ATAT) the

Y) of the base attached

Z), the discussion

torsion

parameters

whereas the dT sugars are characterized

of two neighbouring

parameters

years a large num-

sequence-dependent

that this

Study of the column headed Mobs to correlate

A clear

170°,

reveals

R or pyrimidine

but to the character

the geometrical

261. For example,

is scanty

Some of the main find-

by means of pseudorotation

supra. [25,

scrutiny

backbone angle 6 (eq. this

vi&

geometry is evident

dA sugars have PS values

tion,

obtained

In recent

in our laboratory.

Table II lists

sugars,

DNA fragments

consequence.

If the

is assumed to remain change in relative

viewed in a direction

ori-

121

-+-c-rT--R--T--R 0

0 .I

-I 0

0

-1

0

-.

b*

0 .,

G-G 1, -1 _f

&.*I

.I

-I

1

&.

-I 0

0

-1

0

. 0 f,

0

Fig. 8. Correlation of the variation of the backbone angle 6 observed for a series of single-helical DNAfragments in aqueous solution with the expected variation according to the sum function C6. Experimental values are connected by full lines, predicted ones by broken lines. The unit of change in Lg corresponds to So of 6. (From ref. 26, see also ref. 25 and Table II). perpendicular ly identical

to the base planes. to the local

From the model depicted

twist in Fig.

The relative

orientation

so defined

is virtual-

angle tl as defined by Dickerson and Drew [28]. 9 it is seen that a change from 6 152” to 6 133”

122

Fig. 9. Schematic illustration of the effect of 6 variation on base-base overlap in a regular B-type single-helical dA-dT step, viewed in the 5’+3’ direction. The broken outline shows the expected position of the dT base for 6 152’, i.e. the approximate 6 value of the thymidine monomer. A decrease of the dT 6 angle to 133’ (Table II) causes an appreciable sliding of the dT base with respect to the overlying dA base (solid line), concomitant with a decrease of the local twist angle tl from about 42’ to about 25’: allows

the lower base to slide

about 25’,

i.e.,

helical

dependent adjustments stitutes

in such a way that tl

“untwisting”.

These considerations

of the deoxyribose

a mechanism for

fine

decreases

suggest

expressed

pucker,

tuning of the base-base

from about 42’

to

that sequence

as 6 variation,

con-

overlap.

Double helices DNAs can adopt several mental conditions DNAduplexes E,T),

double-helical

such as counterions,

are classified

the B family

this

that several

deoxyoligonucleotides

clearly

mode is S (2’-endo).

indicate

At this

point

sequence d(GGCCGGCC)[29,30]. in solution

indicated contrast

a major contribution to present

Thus far,

lines

the sole

is read off

Similarly, of B-like

representative

These duplexes

display

interesting

[33] to steric

the duplex formed by

It is well

to remember that

clashes

N- or S-type

RNA-DNAhybrid

forms in the deoxyribose

studies

[18]

strands,

in

[3]. of the B-type duplex in the single

is the dodecamer sequence d(CGCGAATTCGCG) and its by Calladine

as

compounds in aqueous

in the A form [32] whereas solution

of thought

to note

and characterized

from two criteria:

the covalent

(B,C,D,

whereas for

is of interest

This concerns

the self-complementary

[31].

the B-family

of one of these

A- or B-type character sugar and NOE intensities

it

have been crystallized

B-type character.

r(GCG)d(TATACGC) crystallizes

sequence and composition.

A-type,

mode for A DNA.is N (3’-endo),

DNA, whereas our NMRinvestigations

solution

depending upon environ-

water activity,

in three broad groups:

and Z DNAs. The sugar puckering

A-type

conformations,

local

C(9) brominated

changes of geometry

between adjacent

purines

crystal

derivative.

[27,28,34] in opposite

ascribed strands.

123

d(C A C A T G T Gl

d(C A C Al d(T G T G)

Fig. 10. Populations of S conformer (%, hatched areas) in two single-helical tetramers - d (CACA) and d (TGTG) - and in the composite octamer duplex d(CACATGTG), ref. 35. Note that the sugars in the core of the duplex appear to retain some conformational freedom. Obviously,

such clashes

can occur

posed to be exceptionally mations are concerned, (unimodal) (P 90°)

variation

tend to favour

and coworkers

of sugar puckers,

In the crystal

about Cl’-exo

(P 126’)

favour

a similar

picture,

prefer

albeit

d(CACATGTG), together

high 6 values.

Moreover,

with those

in the two constituent

single

[35] by means of line

The results

(Fig.

lo),

it

of the duplex existence

still

distribution). studies hybrid

retain

flip

duplex octamer d(GGm5Cm5C)2[17],

alternating

already

influence

pared to that observed

the geometrical a considerable offered

principle

parameters

since

a certain

found.

the marked decrease

RNA-DNA

amount of residual

that a limited Second (Fig. 6 variation,

fragments.

by the Calladine-Dickerson

[27,33],

(bimodal

our earlier

amount

at Zft resolution,

amount of additional single-helical

confor-

Thus, the

and on the covalent

in an X-ray analysis

in the constituent

to the rationale

anticomplementarity

freedom.

It must be mentioned in passing

duplex induces

on

the sugars in the core

duplex is indicated

revealed

12”

10 and 11.

towards the S-type

does not come as a surprise,

might remain undetected

though it could

contrast

(S Z N) in the intact

duplex mentioned above [IS],

freedom of the sugar rings. of ring

driven

tetramers,

some measure of conformational flip

This conclusion

on B-DNA

although

single-helical

fragments X2’,

in the section

are summarized in Figures

is noted that,

of dynamic ring

helical

widths Cl’,

with the aid of the method outlined

First,

[35]

dY-dR octamer duplex

conformational

mer compared to the constituent

(omitting

pyrimidines

one.

taken from COSY spectra, analysis.

to 04’-endo

Our NMRstudies

sugars in the alternating

d(CACA) and d(TGTG), were analyzed

(P 162’)

with a mean 6 = 125’

a more complicated

The puckers of the individual

a continuous

the sugar conformations

and a 6 range running from 155’ to 80’.

low 6 and purines

and are sup-

As far as the sugar confor-

[27,28,34]

running from C2’-endo

6 variation.

distributed

the 3’ terminals),

in the minor groove.

Dickerson

with a concomitant

are apparently

reveal

strong

only in dR-dY and in dY-dR steps

rules,

alll),

the com-

However, in i.e.,

the

in 6 of the pyrimidine

124

110-L

CACATGTG

Fig. 11. The sequence-dependent variation of backbone angle 6 of the major Stype conformation in two single-helical tetramers d(CACA) and d(TGTG), indicated deduced for the B-type double helix formed by +, compared to the 6 variation by the octamer d(CACATGTG), indicated by o. residues

dC and dT is not accompanied by a similar

residues

dA and dG (perhaps with the exception

the &-range seen in aqueous solution smaller

than the b-range

(155O to 80°, (134’) this

omitting

therefore point

of local

it

strand steric generalizations

(157” to 115’)

the 3’ terminals).

can be concluded

Moreover,

of the dodecamer

The value of Gmean in our octamer

than that found for the dodecamer (125’). that the NMRresults

of geometry in a B-DNA duplex,

effects.

in 6 of the purine

appears to be decidedly

deduced from X-ray crystallography

remains greater

variations

increase

of the dG(6) sugar).

More solution

studies

decidedly possibly

favour

At

the concept

mediated tria inter-

appear necessary,

however,

before

can be made.

ACKNOWLEDGEMENTS This research was supported by the Netherlands Foundation for Chemical Research (S.O.N.) with financial aid from the Netherlands Organization for the Advancement of Pure Research (Z.W.O.). Dr J.-R.Mellema, Mr L.P.M.Orbons and Mr L.J.Rinkel kindly supplied some of the illustrations used in this paper. REFERENCES P.O.P. Ts’o, “Basic Principles in Nucleic Acid Chemistry”, vols 1-3, Academic Press, New York (1974). C.R. Cantor and P.R. Schimmel, “Biophysical Chemistry”, vols 1-3, Freeman, San Francisco (1980) . of Nucleic Acid Structure”, Springer, New York (1984). W. Saenger, “Principles C. Altona, Reel. Trav. Chim. Pays-Bas 101 (1982) 413-433. C. Altona, in Stud. Org. Chem. Vol. 20 “Natural Products Chemistry”, Elsevier, Amsterdam (1985) 285-304. C. Altona and M. Sundaralingam, J. Am. Chem. Sot., 94 (1972) 8205-8212; idem ibid. 95 (1973) 2333-2344. IUPAC-IUB Joint Commission on Biochemical Nomenclature, Eur.. J. Biochem., 131 (1982) 9-1s. F.A.A.M. de Leeuw, P.N. van Kampen, C. Altona, E. Diez, and A.L. Esteban, J. Mol. Struct ., 125 (1984) 67-88. H.P.M. de Leeuw, C.A.G. Haasnoot, and C. Altona, Isr. J. Chem.,20 (1980) 108-126.

125

10 C.A.G. Haasnoot, F.A.A.M. de Leeuw, and C. Altona, Tetrahedron, 36 (1980) 2783-2792. 11 C.A.G. Haasnoot, F.A.A.M. de Leeuw, H.P.M. de Leeuw, and C. Altona, Biopolymers, 20 (1981) 1211-1245. 12 F.A.A.M. de Leeuw and C. Altona, J. Comput. Chem., 4 (1983) 428-437; iciem, Quant. Chem. Progr. Exch. No 463 (1983); idem, J. Chem. Sot. Perkin II, (1982) 375-384. 13 G.M. Clore and A.M. Gronenborn, EMBO J., 4 (1985) 829-835; i&m, FEBS Lett., 179 (1985) 187-198 and references therein. 14 .J.W.Keepers and T.L. .James,J. Magn. Reson. 57 (1984) 404-426: M.S. Broido, T.L. James, G. Zon and J.W.Keepers, Eur. J. Biochem., 150 (1985) 117-128. 15 A.M. Gronenborn, G.M. Clore and B.J. Kimber, Biochem. J., 221 (1984) 723-736. 16 A. Bax, "Two-Dimensional Nuclear Magnetic Resonance in Liquids", Reidel, Dordrecht, Holland, 1982. 17 M.R. Sanderson, J.-R. Mellema, G.A. van der Marel, G. Wille, J.H. van Boom, and C. Altona, Nucleic Acids Res. 11 (1983) 3333-3346. 18 J.-R. Mellema, C.A.G. Haasnoot, G.A. van der Marel, G. Wille, C.A.A. van Boeckel, J.H. van Boom and C. Altona, Nucleic Acids Res. 11 (1983) 5717-5738. 19 COSY = Correlated Spectroscopy, A. Bax, R. Freeman and G.A. Morris, J. Magn. Reson.,42 (1981) 164-168; NOESY q Nuclear Overhauser Enhancement Spectroscopy, J. Jeener, B.H. Beier, P. Bachmann and R.R. Ernst, .J.Chem. Phys., 71 (1979) 4546-4553. 20 DAYCOSY = double-RELAY experiment, A. Bax and G. Drobny, J. Magn. Reson., 61 (1985) 306-320; G. Eich, G. Bodenhausen and R.R. Ernst, J. Am. Chem. Sot., 104 (1982) 3731-3733. 21 ROESY = Rotating-frame Overhauser Enhancement Spectroscopy, A. Bax and D.G. Davis, J. Magn. Reson., 63 (1985) 207-213; A.A. Bothner-By, R.L. Stephens, J. Lee, C.D. Warren and R.W. Jeanloz, J. Am. Chem. Sot., 106 (1984) 811-813. 22 J. Feigon, J.M. Wright, W. Leupin, W.A. Denny and D.R. Kearns, J. Am. Chem. sot., 104 (1982) 5540-5541. 23 D.G. Reid, S.A. Salisbury, S. Bellard, Z. Shakked and D.H. Williams, Biochemistry, 22 (1983) 2012-2025. 24 R.M. Scheek, N. Russo, R. Boelens, R. Kaptein and .J.H.van Boom, J. Am. Chem. sot., 105 (1983) 2914-2916. 25 J.-R. Mellema, J.M.L. Pieters, G.A. van der Marel, J.H. van Boom, C.A.G. Haasnoot and C. Altona, Eur. .J.Biochem. 143 (1984) 285-301. 26 J.-R. Mellema, Appendix in Thesis, Leiden, 1984. 27 R.E. Dickerson, J. Mol. Biol. 166 (1983) 419-441. 28 R.E. Dickerson and H.R. Drew, J. Mol. Biol. 149 (1981) 761-768. 29 A.H.-J. Wang, S. Fujii, J.H. van Boom and A. Rich, Proc. Natl. Acad. Aci. U.S.A., 79 (1982) 3968-3972. 30 L.J. Rinkel, M.R. Sanderson and C. Altona, to be published. 31 C.A.G. Haasnoot, H.P. Westerink, G.A. van der Marel, and J.H. van Boom, J. Biomol. Struct. Dyn., 2 (1984) 3445-3460. 32 A.H.-J. Wang, S. Fujii, J.H. van Boom, G.A. van der Marel, C.A.A. van Boeckel and A. Rich, Nature, 299 (1982) 601-604. 33 C.R. Calladine, J. Mol. Biol. 161 (1982) 343-352. 34 A.V. Fratini, M.L. Kopka, H.R. Drew and R.E. Dickerson, J. Biol. Chem., 257 (1982) 14686-14707. 35 L..J.Rinkel, J.-R. Mellema, L.P.H. Orbons, G.A. van der Marel, J.H. van Boom and C. Altona, Abstracts 4th Conversation in Biomolecular Stereodynamics, R.H. Sarma @d.)State University of New York at Albany, 1985, pp 126-127; L.J. Rinkel et al. to be published.