Structural and Motional Changes Induced in apo-S100A1 Protein by

Systematic, structural studies of various mixed disulfides of S100A1 in the apo and holo ... S100A1 protein, intensively studied since its discovery i...
0 downloads 0 Views 512KB Size
640

Biochemistry 2008, 47, 640-650

Structural and Motional Changes Induced in apo-S100A1 Protein by the Disulfide Formation between Its Cys 85 Residue and β-Mercaptoethanol†,‡ Igor Zhukov,§,| Andrzej Ejchart,§ and Andrzej Bierzyn´ski*,§ Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawin´ skiego 5A, 02-106 Warsaw, Poland, and SloVenian NMR Centre, National Institute of Chemistry, HajdrihoVa 19, P.O.B. 660, SI-1001 Ljubljana, SloVenia ReceiVed September 4, 2007; ReVised Manuscript ReceiVed October 26, 2007

ABSTRACT: Recently, we have shown (Goch, G., Vdovenko, S., Kozłowska, H., and Bierzyn´ski, A. (2005) FEBS J. 272, 2557-2565) that the chemical modification of Cys 85 residue of S100A1 protein by disulfide bond formation with small thiols such as glutathione, cysteine, or β-mercaptoethanol (βME) leads to a dramatic increase of the protein affinity for calcium. Therefore, the biological function of S100A1 as a calcium signal transmitter is probably regulated by the redox potential within the cell. Systematic, structural studies of various mixed disulfides of S100A1 in the apo and holo states are necessary to elucidate the mechanism of this phenomenon. Using NMR methods we have determined the structure of apo-S100A1-βME and, on the basis of 15N nuclear magnetic relaxation data, we have characterized the structural dynamics of both: modified and unmodified molecules of apo-S100A1. The following effects of βME modification have been observed: (1) Helices IV and IV′ of two protein subunits are elongated by five residues (85-89). (2) Conformation of the calcium binding N-terminal loops is dramatically changed, and structural flexibility of the N-loops as well as C-loops markedly increases. (3) The angle between helices I and IV increases by ∼20° and between helices IV and IV′ decreases by ∼35°. All these observations lead to the conclusion that βME modification of apo-S100A1 makes its structure more similar to that of holo-S100A1, so that it becomes much better adjusted for calcium coordination.

Calcium ion is one of the most important messengers regulating numerous vital biological processes. A crucial role in the calcium signal transduction is played by EF-hand proteins that coordinate calcium and change their conformation exposing hydrophobic patches to which a large variety of target proteins bind. S100 is a subfamily of EF-hand proteins found only in mammals and avians, implicated in cell growth, motility, and cell differentiation, by a variety of processes such as regulation of some enzyme as well as transcription factor activities, Ca2+ homeostasis, and regulation of cytoskeleton dynamics. Some of them have extracellular cytokine-like activities (1, 2). Recently, S100 proteins have received increasing attention due to their close association with several human diseases including cardiomyopathy, neurodegenerative disorders, and cancer (2). S100A1 protein, intensively studied since its discovery in 1965 (3), is, in every respect, a typical representative of the S100 family. The structure of S100A1 in the apo state has been determined by Rustandi et al. (4) using NMR methods. Recently, its structure in the holo state has also been solved (5). The protein molecule is a symmetrical homodimer, † This work was supported by Polish Committee for Scientific Research Grant 6 P04 009 16. ‡ The coordinates have been deposited in the Protein Data Bank (code 2JPT). * To whom correspondence should be addressed. Telephone: (48 22) 592-23-71. Fax: (48 22) 823-71-94. E-mail: [email protected]. § Polish Academy of Sciences. | National Institute of Chemistry.

strongly stabilized by hydrophobic interactions between its two R subunits. Each subunit contains two EF-hand motives that coordinate calcium ions in a noncooperative way with the binding constants, at physiological pH (7.2) and NaCl concentration (100 mM), of the order of 103 M-1 for C-terminal motives and 102 M-1 for N-terminal ones (6). These values are too low for a calcium-signaling protein, because the intracellular calcium concentration of Ca2+ ions never exceeds ≈1 µM, even during the calcium stress (7). We have found, though, that the formation of a mixed disulfide between glutathione, cysteine, or β-mercaptoethanol (βME)1 and the single cysteine residue present in position 85 in both protein subunits leads to a dramatic increase of S100A1 affinity for calcium (6). This observation suggests that the calcium-dependent biological activity of the protein can be regulated by thionylation, most probably glutathionylation, because of high concentration of glutathione in the cell (∼10 mM). Therefore, S100A1 can play the role of a linker between the two most important signaling pathways, i.e., the redox- and calcium-dependent ones. What is the structural mechanism responsible for this phenomenon? To answer this question it is necessary to compare the structures of S100A1 and its mixed disulfides 1 Abbreviations: βME, β-mercaptoethanol; DSS-d4, 3-trimethylsilyl2,2,3,3-tetradeuteropropionic acid sodium salt; EDTA, ethylenediaminetetraacetic acid; HSQC, heteronuclear single quantum coherence; HPLC, high performance liquid chromatography; NMR, nuclear magnetic resonance; NOE, nuclear Overhauser effect; S100A1-βME, mixed disulfide of apo-S100A1 protein with β-mercaptoethanol; TRISd11, perdeuterated 2-amino-2-(hydroxymethyl)1,3-propanediol.

10.1021/bi701762v CCC: $40.75 © 2008 American Chemical Society Published on Web 12/19/2007

Changes Induced in apo-S100A1 by Its Cys 85 Modification in the apo and holo states. Valuable information can also be obtained from comparative studies of the structural dynamics of these molecules. In particular, comparative studies of the protein modified with glutathione and βME, the substituents so different chemically, yet similarly enhancing the protein calcium-binding ability, is interesting, making it possible to discern between relevant and nonrelevant conformational changes induced by a substituent. In this paper we present the NMR structure of apo-S100A1 modified by βME and compare it with the structure of the unmodified protein determined by Rustandi et al. (4). Structural flexibilities of both molecules, as determined by nuclear magnetic relaxation methods, are also compared. MATERIALS AND METHODS Sample Preparation. 15N-labeled and 13C,15N-double labeled S100A1 protein was expressed as previously described (8). Shortly: The synthetic gene coding for bovine S100R subunit with the amino acid sequence determined by Isobe and Okuyama (9) was cloned into a derivative of pAED4 plasmid and expressed in Escherichia coli utilizing the T7 expression system. The bacteria were grown in M9 medium containing (15NH4)2SO4 as the sole nitrogen source and either unlabeled or 13C-labeled glucose as the sole carbon source. The expression products were isolated using a phenylsepharose column, purified by reverse-phase HPLC on a semipreparative Vydac C18 column, and identified by the electrospray ionization mass spectrometry using a Macromass Q-Tof spectrometer (6). Because 1 mM concentration of β-mercaptoethanol was maintained during E. coli cell sonification and isolation of the recombinant protein, its disulfide with β-mercaptoethanol (S100A1-βME) predominates in preparations and can be easily separated by HPLC from the reduced S100A1 molecules. Two forms of the protein, (a) with the sequence strictly corresponding to its gene sequence and (b) with the additional initiator methionine at the N-terminus, come from HPLC as partly overlapping peaks. NMR measurements indicated that structural differences between both forms are small and localized to the close vicinity of the N-terminal Met residue (10). Therefore, a mixture of (a) and (b) species was used in all experiments. 600 µL NMR samples contained 1.1 mM protein solution (monomer concentration) in 90%/10% H2O/D2O, 50 mM TRIS-d11, 2.9 mM EDTA, 0.1 mM NaN3, and 96 mM KClO4 with pH adjusted to 6.3 (uncorrected value). For 13C,15Ndouble labeled protein the sample volume was limited to 250 µL and NMR measurements were performed using a Shigemi tube. NMR Spectroscopy. NMR measurements were performed on Varian Inova 400 MHz (9.4 T), Varian Unity+ 500 MHz (11.7 T), and Varian Inova 750 MHz (17.6 T) spectrometers at 35 °C carefully adjusted with ethylene glycol reference sample. All spectrometers were equipped at least with three channels, a z-axis pulsed field gradient unit, and a triple {13C,15N}1H probe with gradient coil(s) and the inverse detection. All heteronuclear NMR data were acquired using the States-TPPI quadrature detection (11) followed by the sensitivity enhanced detection introduced by Kay et al. (12). 1.4 s recycling delay was used, if not stated otherwise. All chemical shifts in heteronuclear NMR spectra were reported

Biochemistry, Vol. 47, No. 2, 2008 641 Table 1: Parameters of NMR Experiments Used for the Structure Determination of apo-S100A1-βME Proteina

experiment 2D [1H-15N]HSQC 3D HNHA (17.6 T) 3D 15N-edited NOESY-HSQC (17.6 T), 150 ms mixing time 3D HNCO, 3D C(CO)NH, 3D (HCA)CO(CA)NH, 3D H(CCO)NH 3D HN(CO)CA, 3D HNCA, 3D CBCA(CO)NH, 3D HNCACB 3D 13C-edited NOESY-HSCQ, (17.6 T), 100 ms mixing time

dim

nucl

time pts

acq time (ms)

t1 t2 t1 t2 t3 t1

15

N H 1 H 15 N 1 H 1 H

64 512 64 16 512 320

49.2 83.9 12.8 16.0 83.9 35.5

t2 t3 t1 t2

15

N H 13 C 15 N

24 720 32 16

19.3 80.0 16.8 12.3

t3 t1 t2 t3 t1

1

H C 15 N 1 H 1 H

384 32 16 384 320

62.9 9.4 12.3 62.9 35.5

t2 t3

13

24 720

19.3 80.0

1

1

13

C H

1

a If not indicated otherwise, experiments were carried out at 11.7 T. Zero filling in all dimensions was applied prior to the Fourier transformation.

with respect to external DSS-d4. Chemical shifts of 13C and N signals were referenced indirectly using the following ratios of the zero-point frequencies: f(13C)/f(1H) ) 0.251449530 and f(15N)/f(1H) ) 0.101329118 (13). Experimental data were processed using the NMRPipe software package (14). The processed spectra were analyzed with both the XEASY (15) and SPARKY (16) software. The sequence-specific assignment was performed on the uniformly 13C,15N-double labeled sample of apo-S100A1βME. The 3D heteronuclear HNCACB (17), CBCA(CO)NH (18), HNCA (19), HN(CO)CA (20), HNCO (21), (HCA)CO(CA)NH (22) spectra were used to obtain assignments of backbone 1H, 13C, and 15N resonances using standard methods (21). Aliphatic side chain 1H and 13C resonances were assigned from the analysis of the heteronuclear 3D C(CO)NH and H(CCO)NH (23), and 2D [1H13 C] HSQC spectra. Aromatic side chain resonances were assigned from the 13C-edited NOESY-HSQC spectrum recorded with the offset, spectral widths, and 13C-1H coupling constants tuned to aromatic carbons. 74 backbone scalar coupling constants 3J(HN,HR) were determined using the HNHA experiment. The detailed information about all experiments is collected in Table 1. 1H, 13C, and 15N chemical shifts were deposited in the BMRB database (http://www.bmrb.wisc.edu) under the accession number 4982. Assignments of HN, N, C′, CR, and Cβ chemical shifts in apo-S100A1 were performed using the same techniques as in the case of apo-S100A1-βME protein. Relaxation experiments were performed at two magnetic fields, 9.4 and 11.7 T, on 15N-labeled samples of either apoS100A1 or apo-S100A1-βME protein. Longitudinal (R1) and transverse (R2) relaxation rates were measured using a sensitivity enhanced [15N-1H] HSQC pulse sequence (12) with the option of either R1 or R2 measurements of 15N nuclei (24). Each spectrum was acquired with 512 (t2) × 128 (t1) 15

642 Biochemistry, Vol. 47, No. 2, 2008 complex data points and the sweep widths of 6000 Hz in 1H and 1300 Hz in 15N dimensions. Zero filling was performed prior to the Fourier transformation. In the case of R1 eleven relaxation delays of 10, 30, 90, 130, 210, 290, 370, 490, 610, 730, and 850 ms were used. The R2 relaxation rate measurements were performed with the Carr-PurcellMeiboom-Gill (CPMG) pulse train. Refocusing time of 650 µs was used during the relaxation delays which were set to 10, 30, 50, 70, 90, 130, 170, 210, and 250 ms. The acquisition parameters for R1 and R2 measurements were identical with the exception of the delay between π (1H) pulses used for the cross-correlation effect suppression (25). 5 and 10 ms delays were used in R1 and R2 measurements, respectively. Delay between the scans of 1.8 s was employed in both experiments. {1H}-15N NOEs were measured with the pulse sequence included in the ProteinPack Varian Inc. (Palo Alto, CA) software. To increase its accuracy the dynamic variant of the experiment was chosen (26). The 1H saturation time was set to 0.1, 0.15, 0.3, 0.7, 1.3, and 2.1 s. 32 scans per t1 increment and the recycling delay of 3.2 s were employed. In all relaxation measurements the 15N decoupling during the acquisition was applied using the GARP sequence with γB2/2π ) 3200 Hz (27). All spectra were processed using the NMRPipe (14) and analyzed with the XEASY (15) programs. Structure Calculations. The cross-peak assignment in NOESY spectra of S100 proteins is seriously complicated by poor resonance dispersion owing to a high content of the R-helical secondary structure and by the ambiguity arising from the symmetrical, dimeric structure of their molecules. Nevertheless, the cross-peak dispersion in the 1H/15N-HSQC spectrum of S100A1-βME protein was fairly good (Figure 1), much better than in the case of apo-S100A1. Therefore, the 15N-edited NOESY-HSQC and 13C-edited NOESYHSQC spectra of S100A1-βME performed at high magnetic field strength (17.6 T) were informative enough to assign unequivocally intersubunit constraints among protons situated on the subunit interface of the protein (helices I, I′, IV, and IV′), so that measurements of 4D spectra were not necessary, although they had been indispensable for structure determination of apo- as well as holo-S100A1 protein (4, 5). Strong (1.8-2.9 Å), medium (1.8-3.5 Å), weak (1.84.5 Å), and very weak (1.8-5.5 Å) intrasubunit as well as intersubunit distance constraints were assigned from crosspeaks in the heteronuclear 3D 15N-edited NOESY-HSQC (150 ms mixing time) and 13C-edited NOESY-HSQC (100 ms mixing time) spectra. The pseudo-atom corrections were applied to the upper limits of constraints for methyl, methylene, and aromatic ring protons (28). +0.3 Å correction was added to upper limits of distance constraints involving methyl groups to account for a higher intensity of methyl resonances (29). In order to distinguish between intra- and intersubunit NOE-derived distance constraints, the protocol introduced by Nilges (30) was used, resulting in identification of 286 intersubunit constraints. Backbone dihedral constraints (-80° < φ < -40° and -55° < ψ < -25° for R-helices and -160° < φ < -80° and 100° < ψ < 180° for β-strands) were included in the structure calculation on the basis of 3 J(HN,HR) scalar coupling constants (J < 5 Hz for R-helices and J > 8 Hz for β-strands) (28). Additional constraints for φ and ψ backbone torsion angles were evaluated using the

Zhukov et al.

FIGURE 1: Amide correlations in the [15N-1H] HSQC spectrum of apo-S100A1-βME protein at 11.7 T. Pairs of side chain primary amide correlations in asparagines and glutamines are connected by horizontal solid lines.

TALOS software package (31). The predictions for 70 residues out of 93 were classified as “good” and used in the structure calculations after the fivefold increase of their deviations. 230 constraints for χ1 and χ2 torsion angles were evaluated by the GLOMSA and HABAS (gridsearch procedure) programs (32) included in the CYANA 1.0.6 software (33). The distance restraints (rHN-O ) 1.5-2.8 Å and rN-O ) 2.4-3.5 Å) for 72 backbone hydrogen bonds in regions of regular secondary structures (28) were included in the final stage of structure calculations. This procedure was the same as that applied previously by Wright et al. during the structure calculations of holo-S100A1 (5), which makes more trustworthy the comparisons between the structure solved by us and that of holo-S100A1 published by Wright et al. (PDB accession no. 1ZFS). Determination of the spatial structure of apo-S100A1-βME was done using the XPLOR-NIH package based on the X-PLOR software (34), extended with the NMR-specific modules (35). Calculations started from a monomer template structure with randomly chosen φ and ψ torsion angles. In the first step, the structure was duplicated providing two identical, completely overlapping monomers. The standard calculation protocol included simulated annealing (12 ps at 3000 K) followed by slow cooling to 100 K combined with the final refinement procedure. The potential energy term for a noncrystallographic symmetry (NCS) was included in the simulated annealing procedure as described by Nilges (30).

Changes Induced in apo-S100A1 by Its Cys 85 Modification The β-mercaptoethanol molecule was created using the standard molecular geometry and introduced as a hetero compound to the X-PLOR topology and parameter files (34). The generation of topology for apo-S100A1-βME protein was done by application of a specific patch to create the S-S bond between βME and Cys 85. Analysis of 15N Relaxation Data. All relaxation parameters (R1, R1, NOE) were determined by a two-parameter nonlinear least-squares fit of cross-peak heights to a single-exponential decay. Two relaxation mechanisms have to be taken into account for the nuclear spin relaxation of amide nitrogens in proteins: the chemical shift anisotropy of a nitrogen nucleus and the dipole-dipole interaction between a nitrogen and the hydrogen directly bound to it. Equations describing relaxation parameters in terms of spectral density functions are given by (36)

1 R1 ) D2[J(ωH - ωN) + 3J(ωN) + 4 1 6J(ωH + ωN)] + C2J(ωN) 3 1 R2 ) D2[4J(0) + J(ωH - ωN) + 3J(ωN) + 8 1 6J(ωH) + 6J(ωH + ωN)] + C2[4J(0) + 3J(ωN)] + Rex 18 D2 NOE ) 1 + (γH/γN) [6J(ωH + ωN) - J(ωH - ωN)] 4R1 where D ) [µ0/(4π)][(γNγHp)/〈rNH3〉], C ) ωN∆σ, and other symbols have their usual meaning. The values of rNH ) 0.105 nm and ∆σ ) -160 ppm (36) were used in calculations. The additional term Rex takes into account the conformational exchange contribution to R2 resulting from processes in the micro- to millisecond time scale often referred to as chemical exchange effects (36, 37). Such processes, slower than the molecular tumbling, but fast enough to average chemical shifts, can influence transverse relaxation rates determined using the CPMG method. Rex contribution to the transverse relaxation rate is proportional to the square of the chemical shift difference between exchanging states, ∆δ, and to ωN, the Larmor frequency. It should be pointed out that the conformational exchange mechanism can affect the apparent transverse relaxation rate only if ∆δ * 0. In the case of an isotropic molecular tumbling the model-free approach spectral density function takes the form (38)

2 J(ω) ) [S2τR/[1 + (ωτR)2] + (1 - S2)τ/[1 + (ωτ)2]] 5 with τ-1 ) τR-1 + τint-1. The isotropic overall motion is described by the correlation time τR and internal motion(s) by the generalized order parameter S, which is a measure of spatial restriction of the motion and of the effective correlation time τint describing the rate of this motion. In the case of an anisotropic overall motion the spectral density function becomes more complex. Combining the model-free approach with an axially anisotropic overall

Biochemistry, Vol. 47, No. 2, 2008 643 Table 2: Diffusion Tensor Parameters Obtained for apo-S100A1-βME and apo-S100A1 Proteins from a Model-Free Analysis Assuming the Axial Symmetry of the Diffusion Tensor protein

apo-S100A1

apo-S100A1-βME

ratios of the principal values of moment of inertia expected ratio of diffusion coefficients experimental D|

0.722:0.823:1.000

0.689:0.671:1.000

0.965

0.950

(2.04 ( 0.17) × 107 s-1 (2.31 ( 0.14) × 107 s-1 0.88 ( 0.09 7.51 ( 0.17 ns

(2.06 ( 0.11) × 107 s-1 (2.14 ( 0.09) × 107 s-1 0.96 ( 0.07 7.89 ( 0.12 ns

experimental D⊥ experimental D|/D⊥ average τRa a

Calculated using the equation τR ) (2D| + 4D⊥)-1.

tumbling (39), the spectral density function takes the form (40)

J(ω) ) ∑kAk(R) JMFA(S,τR,k,τk;ω) where R is the angle between the N-H vector and the unique axis of the rotational diffusion tensor. Three coefficients Ak(R) have the form

A1 ) 0.75(sin4 R), A2 ) 3(sin2 R)(cos2 R), A3 ) 0.25(3cos2 R - 1)2 and corresponding correlation times τR,k are defined as

τR,1 ) 1/(4D| + 2D⊥), τR,2 ) 1/(D| + 5D⊥),

τR,3 ) 1/(6D⊥)

where D| and D⊥ are the parallel and perpendicular components of the rotational diffusion tensor. Atomic coordinates of the NMR-derived structure of S100A1-βME were used to calculate the inertia tensor. The ratios of its principal values are given in Table 2 with the z-axis collinear with the symmetry axis of the homodimeric protein. For S100A1 protein the solution structure of rat apoS100A1 was used (4) (PDB accession no. 1K2H) resulting in similar ratios of principal values of the inertia tensor (Table 2). Both structures were also used to calculate the geometric coefficients Ak(R). Because of the symmetry requirement it was assumed that rotational diffusion of both proteins is well represented by axially symmetric oblate tensors with their unique axis collinear with the protein symmetry axis. The least-squares procedure used to optimize model parameters consisted of minimization through a grid-search of the target function χ2 defined as follows:

χ2 ) ∑i∑j∑k(Pijk,ex - Pijk,calc)2/σijk2 where the sums are over the two magnetic fields used, the number of amino acid residues for which data are available (N), and the three relaxation parameters. The subscript ex refers to the experimentally obtained relaxation parameters and the subscript calc to the respective relaxation parameters

644 Biochemistry, Vol. 47, No. 2, 2008

Zhukov et al.

FIGURE 2: Stereoview of ten lowest total energy structures of apo-S100A1-βME protein.

calculated from the model. σijk values are the standard deviations of the experimentally obtained parameters. At this stage of the analysis the relaxation data for amino acid residues for which only the data at 11.7 T were available (Gly 23, Lys 27, Glu 40, Glu 73, Val 78) or the NOE values were distinctly small (Ser 2 at the N-terminus and Trp 90 through Ser 93 at the C-terminus) were excluded from the calculation. The minimization procedure delivered two global parameters D| and D⊥ and N sets of local parameters S2j, τj, and Rex,j. Uncertainties of the calculated parameter values were estimated as standard deviations from 200 Monte Carlo simulations. Using the D| and D⊥ values determined in this way, the local model-free parameters for 5 residues, previously omitted due to the lack of data obtained at two magnetic fields, were also calculated. On the other hand, the highly mobile terminal residues, Ser 2 and Trp 90 to Ser 93, cannot be regarded as a part of a rigid rotor, and description of their motions in terms of a single rotational diffusion tensor is not appropriate. In such a case the motion of an N-H vector in the ith residue should be described by the model-free approach spectral density functions in which the global parameter τR is substituted for a local one, τR,i (41, 42). RESULTS Structure Determination of apo-S100A1-βME. Initial calculations based on distance constraints derived from the 15 N-edited NOESY-HSQC spectrum provided a medium resolution (∼3.5 Å) three-dimensional structure of the apoS100A1-βME protein. This initial structure was used in the further procedure of intersubunit/intrasubunit cross-peak discrimination in the 3D 13C-edited NOESY-HSQC spectrum. Finally, 1098 intraresidual, 282 sequential, 398 medium range, 248 long range, and 286 intersubunit interproton distance constraints were included in the structure calculations. The overlay of the 10 lowest energy structures of the protein is shown in Figure 2, and the structure calculation statistics is given in Table 3. None of these structures contain NOE violations exceeding 0.3 Å and dihedral angle violations larger than 5°. The analysis made using the PROCHECK-NMR program (43) indicates that the average structure is well-defined, corresponding to the X-ray resolution of 1.6 Å.

Table 3: NMR Constraints and Statistics of 10 NMR-Derived apo-S100A1-βME Structures with the Lowest Total Energy rmsd from distance constraintsa total (2384) intraresidual (|i-j| ) 0) (1098) sequential (|i - j| ) 1) (282) medium range (|i - j| e 5) (398) long range (|i - j| > 5) (248) intersubunits (286) hydrogen bonds (72) rmsd from dihedral angles constraintsa total (558) backbone (φ/ψ) (328) side chain (χ1/χ2) (230) rmsd from idealized covalent geometry bonds (Å) angles (deg) impropers (deg) quality of the structure % of residues in the most favorable region % residues in additional allowed regions % residues in generously allowed regions rmsd to the mean structure ordered backbone atoms (3.90) ordered heavy atoms (3.90) all backbone all heavy atoms equivalent X-ray resolution



best

0.0133 ( 0.0005 0.0078 ( 0.0017 0.0182 ( 0.0019 0.0171 ( 0.0011

0.0126 0.0066 0.0208 0.0163

0.0150 ( 0.0022 0.0152 ( 0.0031 0.0130 ( 0.0026

0.0131 0.0115 0.0131

0.360 ( 0.040 0.462 ( 0.057 0.101 ( 0.042

0.348 0.449 0.092

0.00220 ( 0.00008 0.401 ( 0.005 0.220 ( 0.005

0.00216 0.395 0.221

88.2 ( 2.1

86.0

9.4 ( 3.0

12.8

2.5 ( 1.5

1.2

0.71 ( 0.25 Å

0.62 Å

1.11 ( 0.32 Å 0.92 ( 0.31 Å 1.35 ( 0.38 Å 1.6 Å

0.99 Å 0.93 Å 1.34 Å

a None of the 10 structures has a distance violation greater than 0.3 Å or a dihedral angle violation greater than 5°.

Location of the R-helical structure elements within the protein subunits have been unequivocally determined by the characteristic NOE patterns, (i f i + 3) and (i f i + 4), obtained from the 3D NOESY-HSQC spectra as well as by φ and ψ backbone torsion angles typical for the R-helical conformation. Identification of the R-helical elements is also supported by many medium- and long-range contacts. So, the four helices present in each of the subunits are the best defined regions in the whole protein structure (Table 4). As in other S100 protein structures solved so far, the conformation of the linker between helices II and III and, in

Changes Induced in apo-S100A1 by Its Cys 85 Modification

Biochemistry, Vol. 47, No. 2, 2008 645

Table 4: Rmsd Values of Secondary Structure Elements in 10 NMR-Derived Structures of apo-S100A1-βME Protein with the Lowest Total Energy residues

structure elements

backbone atoms (Å)

all heavy atoms (Å)

3-17 30-40 52-63 72-89

helix I helix II helix III helix IV

0.23 ( 0.10 0.14 ( 0.06 0.19 ( 0.07 0.40 ( 0.17

0.50 ( 0.15 0.72 ( 0.23 0.64 ( 0.20 0.81 ( 0.26

particular, the conformations of the Ca2+-binding loops can be determined much less precisely because of few long-range contacts observed for these protein segments. Nevertheless, from the NOE spectra, and from the values of backbone torsion angles, it is evident that residues Lys 27, Leu 28, and Ser 29 in the N-loop and Glu 68, Val 69, and Asp 70 in the C-loop form an antiparallel β-sheet stabilized by a pair of hydrogen bonds between the amide and carbonyl groups of Leu 28 and Val 69. Analogous structures (44), common for EF-hand proteins, have been found in all S100 proteins studied so far. The three-dimensional structure of the protein can be characterized as an “X-type” bundle with a large intersubunit interface formed by two pairs of antiparallel helices I, I′ and IV, IV′, similarly as in other S100 proteins (Figure 3). Protons of residues situated at the subunit interface display numerous intersubunit NOEs. For example, protons in residues of helix I (Leu 4, Met 8, Leu 11, Ile 12, Val 14) give many contacts to protons in corresponding residues of helix I′ and protons in helix IV (Val 75, Val 76, Val 78, Ala 79, Thr 82) give contacts to residues of helix IV′. A number of cross-peaks have been assigned as arising from contacts between residues located at the N-terminus of helix I and the C-terminus of helix II′. For instance, H protons of Met 8 in helix I give cross-peaks with protons situated not only in helix IV′ (Hδ of Leu 81 and Hγ of Val 78) but also in the C-terminus of helix II′ (Hδ of Leu 37) and in the hinge region (Hδ of Leu 41 and H of Phe 44). Chemical Shift Changes Induced by βME Modification. It is well-known that chemical shifts depend on local protein structure (38, 45, 46). Therefore, chemical shift differences between S100A1 and S100A1-βME molecules can be treated as indicators of these regions of the protein structure that are the most affected by the Cys 85 modification. It has been shown (47) that a weighted combination of HN, C′, N, and CR chemical shift differences, ∆δtot, is a very effective measure of structural changes in proteins. ∆δtot values for the S100A1/S100A1-βME pair calculated using eq 1 given in ref (47) are shown in Figure 4. There are three parts of the protein where they substantially exceed the measurement error limits: Met 8 to Lys 25 segment comprising a large part of helix I and the N-terminal half of the N-loop, Glu 40 to Val 54 segment comprising the linker and the N-terminus of helix III, and Val 83 to Trp 90 segment of helix IV in the vicinity of the modification site (Cys 85). These data can be correlated with global conformational changes induced by the protein modification presented in Figure 5. Rotational Diffusion. R1, R2 and NOE data have been obtained for all backbone amide nitrogen nuclei of apoS100A1-βME protein at 11.7 T and for 86 ones at 9.4 T. For apo-S100A1 the relaxation data for only 72 out of 92 backbone amide nitrogen nuclei have been determined at 11.7

FIGURE 3: Ribbon drawing of the three-dimensional structure of apo-S100A1-βME protein. In both subunits the consecutive helices are labeled with Roman numbers and indicated by different colors. Helices I, I′, IV, and IV′, forming the intersubunit interface, are shown in light colors. The structure is oriented in the same way as those presented in Figures 2 and 5.

T and 67 ones at 9.4 T due to strong overlap of many resonances in the 1H/15N correlation spectrum, predominantly coming from amino acid residues located within helix I, the N-loop, and the C-terminus (see Supporting Information). Ratios of principal values of moments of inertia can be used to calculate the predicted ratios of diffusion coefficients D|/D⊥ (48). For S100A1-βME and S100A1 these ratios equal to 0.950 and 0.965, respectively (see Table 2). This prediction was confirmed by the analysis of experimental data using the anisotropic diffusion model. Experimentally determined diffusion coefficients yield D|/D⊥ ratios very close to the predicted ones (see Table 2). Overall tumbling of both proteins is slightly anisotropic and well represented by the oblate tensors. The averaged isotropic rotational correlation times, τR, are typical for proteins of this size and correspond well to the result obtained for an another protein from the S100 family, i.e., S100B: 7.87 ns (49). Protein Backbone Dynamics. There is no significant differentiation between local dynamics of various structural elements of apo-S100A1. The R-helices as well as the calciumbinding loops and the linker segments display similar average values of S2 parameter, which is a measure of fast local motions in the picosecond time scale (Table 5). On the other hand, local dynamics in apo-S100A1-βME is more differentiated. Although the dynamics of helices III and IV and of the β-sheet structure is not affected by the protein modification, the restriction of N-H vector motions in helices I and II as well as in the linker seems to increase slightly. The most pronounced effects concern both calcium-binding loops N and C which become the “hot spots” of the protein, displaying the greatest motional freedom of their N-H vecors (see Table 5 and Figure 6). Another hot point appears at the linker-helix III junction (Val 51, Asp 52, and Ala 53).

646 Biochemistry, Vol. 47, No. 2, 2008

Zhukov et al.

FIGURE 4: The combined chemical shift difference ∆δtot between apo-S100A1 (unpublished results) and apo-S100A1-βME (BMRB database, acc. no. 4982) calculated according to eq 1 of ref 47. Average value of ∆δtot is equal to 0.575, and its standard deviation (0.423) is marked by the horizontal black line. Black bars correspond to those residues that are in the R-helical conformation in both apo-S100A1 and apoS100A-βME. ∆δtot values for residues 86 and 87 are missing because chemical shifts for C′ of Asn 86 and for HN and backbone N of Asn 87 could not be determined in apo-S100A1. Table 5: Weighted Mean Values of the Order Parameter S2 of Structural Elements of Modified and Unmodified apo-S100A1 Proteina structure element (a.a. residue)/no. of residues used in calcn R-helix I (3-17)/11 R-helix II (30-40)/9 R-helix III (52-63)/11 R-helix IV (72-84)b/12 β-sheet (68-70)c/3 central part of N-loop (20-26)/6 central part of C-loop (64-67)/4 part of linker (41-50)/9 residues (3-89)/70

S100A1-βME

S100A1

∆S2

0.95 (0.01) 0.92 (0.01) 0.94 (0.01) 0.95 (0.01) 0.93 (0.02) 0.79 (0.03)

0.88 (0.02) 0.86 (0.02) 0.93 (0.01) 0.93 (0.01) 0.89 (0.02) 0.86 (0.03)

+0.07 ( 0.02 +0.06 ( 0.02

0.76 (0.03)

0.89 (0.02) -0.13 ( 0.04

0.90 (0.01) 0.93 (0.01)

0.85 (0.01) +0.05 ( 0.02 0.91 (0.01)

+0.04 ( 0.03 -0.07 ( 0.04

a

FIGURE 5: Overlay of ensembles of NMR structures of apo-S100A1 (blue) taken from the Protein Data Bank (accession code no. 1K2H) and of apo-S100A1-βME determined in this work (PDB accession code no. 2JPT). In the apo-S100A1-βME structures the residues for which ∆δtot exceeds 0.423 are marked in red and all others in green, including two residues at the C-terminus: Asn 86 and Asn 87 for which the ∆δtot values have not been measured (please compare with Figure 4). The structures are oriented in the same way as those presented in Figures 2 and 3.

For a few residues in the N-loop of apo-S100A1 (Ser 19, Glu 22, Lys 25), transverse relaxation rates R2 are markedly greater than for any other residues (Table SI2) that indicates additional motions in the micro- to millisecond time scale affecting the transverse relaxation rates. Such processes referred to as chemical exchange effects can result in the nonzero Rex term. A few residues in the N-loop of apoS100A1-βME, similarly as in the unmodified protein, display large values of Rex parameter (see Figure 7 and Table SI3).

Standard deviations (in parentheses) and number of amino acid residues used in calculation are also given. Only these residues, for which data are available for both proteins, have been taken into consideration. b Helix IV in apo-S100A1-βME extends up to Phe 89, but only residues corresponding to helix IV in apo-S100A1 have been used in calculations. c Data are not available for residues Lys 27 to Ser 29 of apo-S100A1.

For two residues, Gly 23 and Lys 25, Rex values at 9.4 T equal 3.6 ( 0.8 s-1 and 3.8 ( 0.4 s-1 in apo-S100A1-βME as compared with 1.1 ( 0.3 s-1 and 2.1 ( 0.5 s-1 in apoS100A1, respectively. Most probably, insignificant Rex value for Glu 22 in apo-S100A1-βME (0.4 ( 0.3 s-1) in contrast to apo-S100A1 (4.5 ( 0.6 s-1) is not due to the lack of an exchange process, but rather to a small chemical shift difference between exchanging conformations precluding manifestation of an exchange effect on the transverse relaxation rate. The Rex difference for Glu 22 correlates with the largest backbone amide 15N chemical shift difference of 7.2 ppm (122.5 and 115.3 ppm in apo-S100A1βME and apo-S100A1, respectively) which can be only compared with the 7.0 ppm difference at the modification site (Cys 85). Different chemical exchange rates and/or different populations of the exchanging states in both

Changes Induced in apo-S100A1 by Its Cys 85 Modification

Biochemistry, Vol. 47, No. 2, 2008 647 DISCUSSION

FIGURE 6: Order parameter values (S2) with corresponding error bars as the function of amino acid sequence for apo-S100A1-βME (upper part) and apo-S100A1 (central part). In the lower part the differences ∆S2 ) S2(apo-S100A1-βME) - S2(apo-S100A1) are shown. Black bars correspond to those residues that are in the R-helical conformation in both proteins. The data for the highly mobile terminal residues, Ser 2 and Trp 90 to Ser 93, are not shown since these residues cannot be regarded as a part of a rigid rotor.

FIGURE 7: Chemical exchange terms (Rex) at 9.4 T with corresponding error bars as a function of the amino acid sequence of apo-S100A1-βME (upper part) and apo-S100A1 (lower part). Black bars correspond to those residues that are in the R-helical conformation in both proteins. The data for the highly mobile terminal residues, Ser 2 and Trp 90 to Ser 93, are not shown since these residues cannot be regarded as a part of a rigid rotor.

proteins might be also responsible for the observed differences of Rex values.

From the data presented in this paper the following picture of conformational changes induced in S100A1 protein by β-Me modification of Cys 85 can be derived: (1) The primary effect of the modification is an elongation of helix IV. In apo-S100A1 (4) the last residue in the helical conformation is Ala 84, and the next, Cys 85, plays the role of the c-cap which NH group forms the hydrogen bond with the carbonyl of Leu 81 but which conformation (φ ) -136, ψ ) +30°) is drastically different from helical. In apoS100A1-βME five other residues at the C-terminus assume the helical conformation documented by NOE contacts and 3 J(HN,HR) scalar coupling constant measurements: Cys 85, Asn 86, Asn 87, Phe 88, and Phe 89. An analogous conformational transition at the C-terminus of helix IV of S100A1 protein arises from calcium coordination (5). What mechanism is responsible for helix IV elongation in apo-S100A1-βME? Most probably it is a hydrophobic interaction between the βME moiety and the aromatic rings of Phe 88 and 89. From NOE measurements we could not obtain direct evidence for these contacts because β-mercaptoethanol used by us was not labeled with 13C. Nevertheless, in the low-energy structures of S100A1-βME determined by us the distances between CR of βME (HO-CRH2-CβH2S-) from H2 and Hζ protons of the Phe 89 ring equal merely 4.3 and 4.8 Å, respectively. In our previous paper (6) we suggested that thionylation of the cysteine residue situated at the 12th position from the C calcium binding loop (Cys 85 in S100A1) regulates, maybe, biological functions not only of S100A1 but also of its many homologues, because cysteine is strictly preserved at this particular position in 7 out of 13 S100A1-like proteins isolated from various species. Remarkably, in all of them (but not in other S100 proteins) the 16th position (analogous to Phe 89 in S100A1) is occupied by phenylalanine2 and the 15th position either by phenylalanine or by tyrosine with only one exception of S100A11. It seems, therefore, that in all these proteins thionylation of the cysteine situated at the 12th position from the C calcium binding loop can induce the same effect of elongation of helix IV, as a consequence of hydrophobic interactions between a moiety of a small thiol molecule attached to the cysteine with aromatic side chains of residues situated at positions 15 and 16. (2) The C-terminal segment of helix IV created by the Cys 85 modification, formed by residues 85-89, attracts the N loop of the other subunit of the protein and drastically changes its spatial orientation (see Figure 5). Consequently, the intrinsic conformation of the loop is also changed in a dramatic way. This transition cannot be precisely described because neither in apo-S100A1 (5) nor in apo-S100A1-βME (this paper) has the exact conformation of the N-loop been determined. Nevertheless, the most important data characterizing conformational features of the N-loop in both proteins are available, Viz., φ and ψ angles of Leu 28 and Glu 32 and χ1 angle of Glu 32 have been evaluated (see Table 6) and the values of S2 and Rex parameters have been measured for many residues forming the N-loop (Figures 6 and 7 and Table SI3). 2 With the unique exception of chick S100A11 protein in which Leu was identified at this position although Phe is preserved in four other S100A11 proteins isolated from man, mouse, pig, and rabbit.

648 Biochemistry, Vol. 47, No. 2, 2008

Zhukov et al.

Table 6: Conformational Data Concerning the Structural Scaffold Formed by Leu 28 and Val 69 in S100A1 Protein and of Glu 32 and Glu 73 Residues Providing the Most Important Ligands of Calcium Ions Coordinated by Its N- and C-Loop, Respectively conformation

Leu 28 φ/ψ

Glu 32 φ/ψ/χ1

Val 69 φ/ψ

Glu 73 φ/ψ/χ1

apo-S100A1 (PDB code 1K2H) holo-S100A1 (PDB code 1ZFS) apo-S100A1-bMe (PDB code 2JPT)

+82°/-176° -82°/+167° -133°/+121°

-74°/-56°/-57° -62°/-50°/+48° -50°/-57°/+35°

-140°/+143° -122°/+151° -138°/+139°

-81°/-37°/+74° -60°/-40°/+62° -48°/-23°/+48°

Table 7: Angles between Axes of Helices I, II, and IVa in apo-S100A1b and apo-S100A1-βME (This Work)c

helices I f II I f IV II f IV I f I′ IV f IV′

∆ in apo-S100A1- apo-S100A1∆ in apo-S100A1 βME βME holo-S100A1 119 ( 3 105 ( 2 -46 ( 2 -167 ( 2 178 ( 2

113 ( 2 124 ( 2 -43 ( 3 -170 ( 3 143 ( 5

6(5 19 ( 4 within error within error -35 ( 6

12 ( 4 24 ( 4 17 ( 2 within error -24 ( 5

a Calculated using the Interhlx program (52). The sign convention as in ref (53). b The protein structure was taken from the Protein Data Bank (accession code no. 1K2H). c ∆: Respective changes of interhelical angles induced by the βME modification or calcium binding to S100A1 as reported in ref (5).

In apo-S100A1 the N-loop segment is quite rigid and by no means properly adjusted for coordination of the calcium ion. The conformations of Leu 28sthe central residue of the loopsand of Glu 32 providing the most important, bidentate ligand of Ca2+ ion are drastically different from those in holo-S100A1 (see Table 6). Upon βME modification of the protein the values of φ and ψ angles of both these residues and of χ1 angle of Glu 32 change dramatically and become much closer to those in holo-S100A1. It may be expected that such a conformational transition considerably increases the N-loop affinity of apo-S100A1 for calcium since the proper orientation of central residues of calcium binding loops and of the glutamate residues always present at their C-termini is believed to be of crucial importance for calcium coordination by all EF-hand proteins (50, 51). In S100A1-βME the precise adjustment of the N-loop conformation necessary for calcium binding is additionally facilitated by increased mobility of this protein segment not only in the picosecond time scale (compare the S2 parameter values given in Table 5) but also in the micro- to millisecond time scale as evidenced by the increased Rex values for Gly 23 and Lys 25 (see the section “Protein Backbone Dynamics”). (3) The enhancement of fast internal motions is observed also for a number of residues situated within the C-loop (Table 5) although its conformation does not change significantly upon βME modification (Table 6). (4) In the modified protein the orientation of helix IV is shifted in such a way that the helices forming the interface between the protein subunits (I, I′, IV, and IV′) assume exactly the same relative orientation as in holo-S100A1 protein (compare interhelical angles I/I′, IV/IV′, and I/IV listed in Table 7). CONCLUSIONS Thionylation of Cys 85 residue of apo-S100A1 by β-mercaptoethanol modifies, in many respects, the protein structure. Although it still remains “closed”, unable to bind target proteins, it becomes more similar to that of holo protein,

which explains why the affinity of S100A1 for calcium increases upon the thionylation. Nevertheless, on the basis of the results presented in this paper it is not possible to explain why the Ca2+-binding constants to thionylated S100A1 increase in such a drastic way (1 order of magnitude to its C-loops and as much as 4 orders of magnitude to N-loops) and why it changes dramatically the mechanism of metal binding from noncooperative to strongly cooperative one (6). This intriguing phenomenon can be elucidated by further systematic structural and dynamic studies of a series of properly chosen protein mutants. Such a method was previously applied, with great success, by Chazin and his co-workers to study the mechanism of calcium binding to calbindin D9k (54-56). ACKNOWLEDGMENT The access to a high field (17.6 T) NMR spectrometer obtained in frame of Large Scale Facility (LSF) in the Institute of Molecular Biotechnology (Jena, Germany) is acknowledged. SUPPORTING INFORMATION AVAILABLE Two tables containing backbone amide 15N relaxation rates and NOEs of apo-S100A1-βME and apo-S100A1 at two magnetic fields: 9.4 and 11.7 T. One table containing modelfree approach parameters for individual amino acid residues of both proteins. This material is available free of charge via the Internet at http://pubs.acs.org. REFERENCES 1. Donato, R. (2001) S100: A multigenic family of calciummodulated proteins of the EF-hand type with intracellular and extracellular functional roles, Int. J. Biochem. Cell Biol. 33, 637668. 2. Heizmann, C. W., Fritz, G., and Schafer, B. W. (2004) S100 proteins: Structure, functions and pathology, Front. Biosci. 7, d1356-d1368. 3. Moore, B. (1965) A soluble protein characteristic of the nervous system, Biochem. Biophys. Res. Commun. 19, 739-744. 4. Rustandi, R. R., Baldisseri, D. M., Inman, K. G., Nizner, P., Hamilton, S. M., Landar, A., Landar, A., Zimmer, D B., and Weber, D. J. (2002) Three-dimensional solution structure of the calcium-signaling protein apo-S100A1 as determined by NMR, Biochemistry 41, 788-796. 5. Wright, N. T., Varney, K. M., Ellis, K. C., Markowitz, J., Gitti, R. K., Zimmer, D. B., and Weber, D. J. (2005) The three-dimensional solution structure of Ca2+-bound S100A1 as determined by NMR spectroscopy, J. Mol. Biol. 353, 410426. 6. Goch, G., Vdovenko, S., Kozłowska, H., and Bierzyn´ski, A. (2005) Affinity of S100A1 protein for calcium increases dramatically upon glutathionylation, FEBS J. 272, 2557-2565. 7. Berridge, M. J., Lipp, P., and Bootman, M. D. (2000) The versality and universality of calcium signaling, Nat. ReV. Mol. Cell Biol. 1, 11-21. 8. Bolewska, K., Kozłowska, H., Goch, G., Mikołajek, B., and Bierzyn´ski, A. (1997) Molecular cloning and expression in Escherichia coli of a gene coding for bovine S100A1 protein

Changes Induced in apo-S100A1 by Its Cys 85 Modification and its Glu32 f Gln and Glu73 f Gln mutants, Acta Biochim. Polon. 44, 275-284. 9. Isobe, T, and Okuyama, T. (1981) The amino-acid sequence of the a subunit in bovine brain S-100a protein, Eur. J. Biochem. 116, 79-86. 10. Baldisseri, D. M., Rustandi, R. R., Zhang, Z., Tang, C., Bair, C. L., Landar, A., Zimmer, D. B., and Weber, D. J. (1999) 1H, 13C and 15N NMR sequence-specific resonance assignments for rat apo-S100A1(RR), J. Biomol. NMR 14, 91-92. 11. Marion, D., Ikura, M., Tschudin, R., and Bax, A. (1989) Rapid recording of 2D NMR spectra without phase cycling. Application to the study of hydrogen exchange in proteins, J. Magn. Reson. 85, 393-399. 12. Kay, L. E., Keifer, P., and Saarinen, T. (1992) Pure absorption gradient enhanced heteronuclear single quantum correlation spectroscopy with improved sensitivity, J. Am. Chem. Soc. 114, 10663-10665. 13. Wishart, D. S., Bigam, C. G., Yao, J., Abildgaad, F., Dyson, H. J., Oldfield, E., Markley, J. L., and Sykes, B. D. (1995) 1H, 13C, and 15N chemical shifts referencing in biomolecular NMR, J. Biomol. NMR 6, 135-140. 14. Delaglio, F., Grzesiek, S., Vuister, G. W., Zhu, G., Pfeifer, J., and Bax, A. (1995) A multidimensional spectral processing system based on UNIX pipes, J. Biomol. NMR 6, 277-293. 15. Bartels, C., Xia, T., Billeter, M., Gu¨ntert, P., and Wu¨thrich, K. (1995) The program XEASY for computer-supported NMR spectral analysis of biological macromolecules, J. Biomol. NMR 5, 1-10. 16. Goddard, T. D., and Kneller, D. G. SPARKY 3, University of California, San Francisco. 17. Wittekind, M., and Mueller, L. (1993) HNCACB, a high-sensitivity 3D NMR experiments to correlate amide-proton and nitrogen resonances with the alpha- and beta-carbon resonances in proteins, J. Magn. Reson. 101B, 201-205. 18. Grzesiek, S., and Bax, A. (1992) Correlating backbone amide and side chain resonances in larger proteins by multiple relayed triple resonance NMR, J. Am. Chem. Soc. 114, 62916293. 19. Muhandiram, D. R., and Kay, L. E. (1994) Gradient-enhanced triple-resonance three-dimensional NMR experiments with improved sensitivity, J. Magn. Reson. 103B, 203-216. 20. Bax, A., and Ikura, M. (1991) An efficient 3D NMR technique for correlating the proton and 15N backbone amide resonances with the R-carbon of the preceding residue in uniformly 15N/13C enriched proteins, J. Biomol. NMR 1, 99-104. 21. Ikura, M., Kay, L. E., and Bax, A. (1990) A novel approach for sequential assignments of 1H, 13C, and 15N spectra for proteins: heteronuclear triple-resonance three-dimensional NMR spectroscopy. Application to calmodulin, Biochemistry 29, 46594667. 22. Lo¨hr, F., and Ru¨terjans, H. (1995) A new triple-resonance experiment for the sequential assignment of backbone resonances in proteins, J. Biomol. NMR 6, 189-197. 23. Grzesiek, S., Anglister, J., and Bax, A. (1993) Correlation of backbone amide and aliphatic sidechain resonances in 13C/15Nenriched proteins by isotropic mixing of 13C magnetization, J. Magn. Reson. 101B, 114-119. 24. Farrow, N. A., Muhamdiram, R., Singer, A. U., Paskal, S. M., Kay, C. M., Gish, G., Shoelton, S. E., Pawson, T., Forman-Kay, J. D., and Kay, L. E. (1994) Backbone dynamics of a free and a phosphopeptide-complexed Src homology 2 domain studied by 15 N NMR relaxation, Biochemistry 33, 5984-6003. 25. Kay, L. E., Nicholson, L. K., Delaglio, F., Bax, A., and Torchia, D. A. (1992) Pulse schemes for removal of the effects of crosscorrelation between dipolar and chemical-shift anisotropy relaxation mechanisms on the measurement of heteronuclear T1 and T2 values in proteins, J. Magn. Reson. 97, 359-375. 26. Zhukov, I., and Ejchart, A. (1999) Factors improving the accuracy of determination of 15N relaxation parameters in proteins, Acta Biochim. Polon. 46, 665-671. 27. Shaka, A. J., Baker, P. B., and Freeman, R. (1985) Computeroptimized decoupling scheme for wideband applications at lowlevel operation, J. Magn. Reson. 64, 547-552. 28. Wu¨thrich, K. (1986) NMR of Proteins and Nucleic Acids, John Wiley, New York.

Biochemistry, Vol. 47, No. 2, 2008 649 29. Clore, G. M., Gronenborn, A. M., Nilges, M., and Ryan, C. A. (1987) Three-dimensional structure of potato carboxypeptidase inhibitor in solution. A study using nuclear magnetic resonance, distance geometry, and restrained molecular dynamics, Biochemistry 26, 8012-8013. 30. Nilges, M. (1993) A calculation strategy fir the structure determination of symmetric dimers by 1H NMR, Proteins: Struct., Funct., Genet. 17, 297-309. 31. Cornilescu, G., Delaglio, F., and Bax, A. (1999) Protein backbone angle restraints from searching a database for chemical shift and sequence homology, J. Biomol. NMR 13, 289-302. 32. Gu¨ntert, P., Braun, W., and Wu¨thrich, K. (1991) Efficient computation of three-dimensional structures in solution from nuclear magnetic resonance data using the program Diana, Caliba, Habas and Glomsa, J. Mol. Biol. 217, 517-530. 33. Gu¨ntert, P., Mumenthaler, C., and Wu¨thrich, K. (1997) Torsion angle dynamics from NMR structure calculation with a new program Dyana, J. Mol. Biol. 273, 283-298. 34. Bru¨nger, A. T. (1992) X-PLOR Version 3.1. A system for X-ray crystallography and NMR, Yale University, New Haven, CT. 35. Schwieters, C. D., Kuszewski, J. J., Tjandra, N., and Clore, G. M. (2003) The Xplor-NIH NMR molecular structure determination package, J. Magn. Reson. 160, 65-73. 36. Korzhnev, D. M., Billeter, M., Arseniev, A. S., and Orekhov, V. Y. (2001) NMR studies of brownian tumbling and internal motions in proteins, Prog. Nucl. Magn. Reson. Spectrosc. 38, 197-266. 37. Stone, M. J., Fairbrother, W. J., Palmer, A. G., III, Reizer, J., Saier, M. H., Jr., and Wright, P. E. (1992) Backbone dynamics of the Bacillus subtilis glucose permease IIA domain determined from 15N NMR relaxation measurements, Biochemistry 31, 43944406. 38. Lipari, G., and Szabo, A. (1982) Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules 1. Theory and range of validity. 2. Analysis of experimental results, J. Am. Chem. Soc. 104, 4546-4570. 39. Woessner, D. E. (1962) Nuclear spin relaxation in ellipsoid undergoing rotational Brownian motion, J. Chem. Phys. 37, 647654. 40. Baber, J. L., Szabo, A., and Tjandra, N. (2001) Analysis of slow interdomain motion of macromolecules using NMR relaxation data, J. Am. Chem. Soc. 123, 3953-3959. 41. Alexandrescu, A. T., and Shortle, D. (1994) Backbone dynamics of a highly disordered 131 residue fragment of staphylococcal nuclease, J. Mol. Biol. 242, 527-546. 42. Brutscher, B., Bru¨schweiler, R., and Ernst, R. R. (1997) Backbone dynamics and structural characterization of the partially folded A state of ubiquitin by 1H, 13C, and 15N nuclear magnetic resonance spectroscopy, Biochemistry 36, 13043-13053. 43. Laskowski, R., Rullmann, J., MacArthur, M., Kaptein, R., and Thornton, J. (1996) AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR, J. Biomol. NMR 8, 477-486. 44. Strynadka, N. C. J., and James, M. N. G. (1989) Crystal structures of the helix-loop-helix calcium-binding proteins, Annu. ReV. Biochem. 58, 951-998. 45. de Dios, A. C., and Oldfield, E. (1996) Chemical Shifts in Biochemical Systems, in Encyclopedia of NMR (Grant, D. M., and Harris, R. K., Eds.) Vol. 2, pp 1330-1336, John Wiley & Sons, Chichester. 46. Wishart, D. S., and Sykes, B. D. (1994) Chemical shifts as a tool for structure determination, Methods Enzymol. 239, 363-392. 47. Ayeda, A., Mulder, F. A. A., Yi, G.-S., Lu, Y., Kay, L. E., and Arrowsmith, C. H. (2001) Latent and active p53 are identical in conformation, Nat. Struct. Biol. 8, 756-760. 48. Mandel, A. M., Akke, M., and Palmer, A. G., III (1995) Backbone dynamics of Escherichia coli ribonuclease HI: correlation with structure and function in an active enzyme, J. Mol. Biol. 246, 144-163. 49. Inmam, K. G., Baldisseri, D. M., Miller, E, K., and Weber, D. J. (2001) Backbone dynamics of the calcium-signaling protein apoS100B as determined by 15N NMR relaxation, Biochemistry 40, 3439-3448. 50. Grabarek, Z. (2005) Structure of trapped intermediate of calmodulin: calcium regulation of EF-hand proteins from a new perspective. J. Mol. Biol. 346, 1351-1366. 51. Grabarek, Z. (2006) Structural basis for diversity of the EF-hand calcium-binding proteins. J. Mol. Biol. 359, 509-525. 52. Ikura, M. (1998) http://nmr.uhnres.utoronto.ca/ikura/interhlx/.

650 Biochemistry, Vol. 47, No. 2, 2008 53. Drohat, A. C., Amburgey, J. C., Abilgaard, F., Starich, M. R., Baldisseri, D., and Weber, D. J. (1996) Solution structure of rat apo-S100B(ββ) as determined by NMR spectroscopy, Biochemistry 35, 11577-11588. 54. Bunick, C. G., Nelson, M. R., Mangahas, S., Hunter, M. J., Sheehan, J. H., Mizoue, L. S., Bunick, G. J., and Chazin, W. J. (2004) Designing sequence to control protein function in an EFhand protein, J. Am. Chem. Soc. 126, 5990-5998.

Zhukov et al. 55. Nelson, M. R., Thulin, E., Fagan, P. A., Forsen, S., and Chazin, W. J. (2002) The EF-hand domain: A globally cooperative structural unit, Protein Sci. 11, 198-205. 56. Ma¨ler, L., Blankenship, J., Rance, M., and Chazin, W. J. (2000) Site-site communication in the EF-hand Ca2+-binding protein calbindin D9k, Nat. Struct. Biol. 7, 245-250. BI701762V