Rotational diffusion anisotropy of proteins from simultaneous analysis of 15N and 13C alpha nuclear spin relaxation

Current methods of determining the rotational diffusion tensors of proteins in solution by NMR spectroscopy exclusively utilize relaxation rate constants for backbone amide 15N spins. However, the distributions of orientations of N-H bond vectors are not isotropic in many proteins, and correlations between bond vector orientations reduce the accuracy and precision of rotational diffusion tensors extracted from 15N spin relaxation data. The inclusion of both 13C alpha and 15N spin relaxation rate constants increases the robustness of the diffusion tensor analysis because the orientations of the C alpha-H alpha bond vectors differ from the orientations of the N-H bond vectors. Theoretical and experimental results for calbindin D9k, granulocyte colony stimulating factor, and ubiquitin, three proteins with different distributions of N-H and C alpha-H alpha bond vectors, are used to illustrate the advantages of the simultaneous utilization of 13C alpha and 15N relaxation data.

1H NMR assignments of apo calcyclin and comparative structural analysis with calbindin D9k and S100 beta

The homodimeric S100 protein calcyclin has been studied in the apo state by two-dimensional 1H NMR spectroscopy. Using a combination of scalar correlation and NOE experiments, sequence-specific 1H NMR assignments were obtained for all but one backbone and > 90% of the side-chain resonances. To our knowledge, the 2 x 90 residue (20 kDa) calcyclin dimer is the largest protein system for which such complete assignments have been made by purely homonuclear methods. Sequential and medium-range NOEs and slowly exchanging backbone amide protons identified directly the four helices and the short antiparallel beta-type interaction between the two binding loops that comprise each subunit of the dimer. Further analysis of NOEs enabled the unambiguous assignment of 556 intrasubunit distance constraints, 24 intrasubunit hydrogen bonding constraints, and 2 x 26 intersubunit distance constraints. The conformation of the monomer subunit was refined by distance geometry and restrained molecular dynamics calculations using the intrasubunit constraints only. Calculation of the dimer structure starting from this conformational ensemble has been reported elsewhere. The extent of structural homology among the apo calcyclin subunit, the monomer subunit of apo S100 beta, and monomeric apo calbindin D9k has been examined in detail by comparing 1H NMR chemical shifts and secondary structures. This analysis was extended to a comprehensive comparison of the three-dimensional structures of the calcyclin monomer subunit and calbindin D9k, which revealed greater similarity in the packing of their hydrophobic cores than was anticipated previously. Together, these results support the hypothesis that all members of the S100 family have similar core structures and similar modes of dimerization. Analysis of the amphiphilicity of Helix IV is used to explain why calbindin D9k is monomeric, but full-length S100 proteins form homodimers.

Sequence dependence and direct measurement of crossover isomer distribution in model Holliday junctions using NMR spectroscopy

A 32-base-pair model of the Holliday junction (HJ) intermediate in genetic recombination has been prepared and analyzed in-depth by 2D and 3D (1)H NMR spectroscopy. This HJ (J2P1) corresponds to a cyclic permutation of the base pairs at the junction relative to a previously studied HJ [J2; Chen, S.-M., & Chazin, W.J. (1994) Biochemistry 33, 11453-11459], designed to probe the effect of the sequence at the n - 1 position (where n is the residue directly at the branch point) on the stacking geometry. Observation of several interbase nuclear Overhauser effects (NOEs) clearly indicates a strong preference for the isomer opposite that observed for J2, confirming the dependence of stacking isomer preference on the sequence at the junction. As for other model HJs studied, a small equilibrium distribution of the alternate isomer could be identified. A sample of J2P1 was prepared with a single (15)N-labeled thymine residue at the branch point. 1D (15)n-filtered (1)H-detected experiments on this sample at low temperature give strong support for the co-existence of the two stacking isomers and provide a much more direct and accurate measure of the crossover isomer distribution. The comparative analysis of our immobile HJs and a model cruciform structure [Pikkemaat, J.A., van den Elst, H., van Boom, J.H., & Altona, C. (1994) Biochemistry 33, 14896-14907] sheds new light on the issue of the relevance of crossover isomer preference in vivo.

Solution structure of (Cd2+)1-calbindin D9k reveals details of the stepwise structural changes along the Apo-->(Ca2+)II1-->(Ca2+)I,II2 binding pathway

The three-dimensional solution structure of (Cd2+)1-calbindin D9k has been determined by distance geometry, restrained molecular dynamics and relaxation matrix calculations using experimental constraints obtained from two-dimensional 1H and 15N-1H NMR spectroscopy. The final input data consisted of 1055 NOE distance constraints and 71 dihedral angle constraints, corresponding to 15 constraints per residue on average. The resulting ensemble of 24 structures has no distance or dihedral angle constraints consistently violated by more than 0.07 A and 1.8 degrees, respectively. The structure is characteristic of an EF-hand protein, with two helix-loop-helix calcium binding motifs joined by a flexible linker, and a short anti-parallel beta-type interaction between the two ion-binding sites. The four helices are well defined with a root mean square deviation from the mean coordinates of 0.35 A for the backbone atoms. The structure of the half-saturated cadmium state was compared with the previously determined solution structures of the apo and fully calcium saturated calbindin D9k. The comparisons were aided by introducing the ensemble averaged distance difference matrix as a tool for analyzing differences between two ensembles of structures. Detailed analyses of differences between the three states in backbone and side-chain dihedral angles, hydrogen bonds, interatomic distances, and packing of the hydrophobic core reveal the reorganization of the protein that occurs upon ion binding. Overall, it was found that (Cd2+)1-calbindin D9k, representing the half-saturated calcium state with an ion in site II, is structurally more similar to the fully calcium-saturated state than the apo state. Thus, for the binding sequence apo-->(Ca2+)II1-->(Ca2+)I,II2, the structural changes occurring upon ion binding are most pronounced for the first binding step, an observation that bears significantly on the molecular basis for cooperative calcium binding in calbindin D9k.

Releasing the calcium trigger

NMR structures of calmodulin, troponin C and related proteins are providing the atomic details of the conformational changes that transduce Ca2+ signals into mechanical or metabolic responses.

The structure of calcyclin reveals a novel homodimeric fold for S100 Ca(2+)-binding proteins

The S100 calcium-binding proteins are implicated as effectors in calcium-mediated signal transduction pathways. The three-dimensional structure of the S100 protein calcyclin has been determined in solution in the apo state by NMR spectroscopy and a computational strategy that incorporates a systematic docking protocol. This structure reveals a symmetric homodimeric fold that is unique among calcium-binding proteins. Dimerization is mediated by hydrophobic contacts from several highly conserved residues, which suggests that the dimer fold identified for calcyclin will serve as a structural paradigm for the S100 subfamily of calcium-binding proteins.

Determination of the solution structure of Apo calbindin D9k by NMR spectroscopy

The three-dimensional structure of apo calbindin D9k has been determined using constraints generated from nuclear magnetic resonance spectroscopy. The family of solution structures was calculated using a combination of distance geometry, restrained molecular dynamics, and hybrid relaxation matrix analysis of the nuclear Overhauser effect (NOE) cross-peak intensities. Errors and inconsistencies in the input constraints were identified using complete relaxation matrix analyses based on the results of preliminary structure calculations. The final input data consisted of 994 NOE distance constraints and 122 dihedral constraints, aided by the stereospecific assignment of the resonances from 21 beta-methylene groups and seven isopropyl groups of leucine and valine residues. The resulting family of 33 structures contain no violation of the distance constraints greater than 0.17 A or of the dihedral angle constraints greater than 10 degrees. The structures consist of a well-defined, antiparallel four-helix bundle, with a short anti-parallel beta-interaction between the two unoccupied calcium-binding loops. The root-mean-square deviation from the mean structure of the backbone heavy-atoms for the well-defined helical residues is 0.55 A. The remainder of the ion-binding loops, the linker loop connecting the two sub-domains of the protein, and the N and C termini exhibit considerable disorder between different structures in the ensemble. A comparison with the structure of the (Ca2+)2 state indicates that the largest changes associated with ion-binding occur in the middle of helix IV and in the packing of helix III onto the remainder of the protein. The change in conformation of these helices is associated with a subtle reorganization of many residues in the hydrophobic core, including some side-chains that are up to 15 A from the ion-binding site.

Quantitative measurements of the cooperativity in an EF-hand protein with sequential calcium binding

Positive cooperativity, defined as an enhancement of the ligand affinity at one site as a consequence of binding the same type of ligand at another site, is a free energy coupling between binding sites. It can be present both in systems with sites having identical ligand affinities and in systems where the binding sites have different affinities. When the sites have widely different affinities such that they are filled with ligand in a sequential manner, it is often difficult to quantify or even detect the positive cooperativity, if it occurs. This study presents verification and quantitative measurements of the free energy coupling between the two calcium binding sites in a mutant form of calbindin D9k. Wild-type calbindin D9k binds two calcium ions with similar affinities and positive cooperativity--the free energy coupling, delta delta G, is around -8 kJ.mol-1 (Linse S, et al., 1991, Biochemistry 30: 154-162). The mutant, with the substitution Asn 56-->Ala, binds calcium in a sequential manner. In the present work we have taken advantage of the variations among different metal ions in terms of their preferences for the two binding sites in calbindin D9k. Combined studies of the binding of Ca2+, Cd2+, and La3+ have allowed us to conclude that in this mutant delta delta G < -6.4 kJ.mol-1, and that Cd2+ and La3+ also bind to this protein with positive cooperativity. The results justify the use of the (Ca2+)1 state of the Asn 56-->Ala mutant, as well as the (Cd2+)1 state of the wild type, as models for the half-saturated states along the two pathways of cooperative Ca2+ binding in calbindin D9k.

Characterization of the N-terminal half-saturated state of calbindin D9k: NMR studies of the N56A mutant

Calbindin D9k is a small EF-hand protein that binds two calcium ions with positive cooperativity. The molecular basis of cooperativity for the binding pathway where the first ion binds in the N-terminal site (1) is investigated by NMR experiments on the half-saturated state of the N56A mutant, which exhibits sequential yet cooperative binding (Linse S, Chazin WJ, 1995, Protein Sci 4:1038-1044). Analysis of calcium-induced changes in chemical shifts, amide proton exchange rates, and NOEs indicates that ion binding to the N-terminal binding loop causes significant changes in conformation and/or dynamics throughout the protein. In particular, all three parameters indicate that the hydrophobic core undergoes a change in packing to a conformation very similar to the calcium-loaded state. These results are similar to those observed for the (Cd2+)1 state of the wild-type protein, a model for the complementary half-saturated state with an ion bound in the C-terminal site (II). Thus, with respect to cooperativity in either of the binding pathways, binding of the first ion drives the conformation and dynamics of the protein far toward the (Ca2+)2 state, thereby facilitating binding of the second ion. Comparison with the half-saturated state of the analogous E65Q mutant confirms that mutation of this critical bidentate calcium ligand at position 12 of the consensus EF-hand binding loop causes very significant structural perturbations. This result has important implications regarding numerous studies that have utilized mutation of this critical residue for site deactivation.

The effect of protein concentration on ion binding

The concentration of protein in a solution has been found to have a significant effect on ion binding affinity. It is well known that an increase in ionic strength of the solvent medium by addition of salt modulates the ion-binding affinity of a charged protein due to electrostatic screening. In recent Monte Carlo simulations, a similar screening has been detected to arise from an increase in the concentration of the protein itself. Experimental results are presented here that verify the theoretical predictions; high concentrations of the negatively charged proteins calbindin D9k and calmodulin are found to reduce their affinity for divalent cations. The Ca(2+)-binding constant of the C-terminal site in the Asn-56 --> Ala mutant of calbindin D9k has been measured at seven different protein concentrations ranging from 27 microM to 7.35 mM by using 1H NMR. A 94% reduction in affinity is observed when going from the lowest to the highest protein concentration. For calmodulin, we have measured the average Mg(2+)-binding constant of sites I and II at 0.325, 1.08, and 3.25 mM protein and find a 13-fold difference between the two extremes. Monte Carlo calculations have been performed for the two cases described above to provide a direct comparison of the experimental and simulated effects of protein concentration on metal ion affinities. The overall agreement between theory and experiment is good. The results have important implications for all biological systems involving interactions between charged species.

Two-dimensional 1H NMR studies of immobile Holliday junctions: nonlabile proton assignments and identification of crossover isomers

The nonlabile protons of two 32-base-pair models of the Holliday junction intermediate in genetic recombination have been studied by two-dimensional 1H nuclear magnetic resonance (NMR) spectroscopy. The sequence of these models is designed to fully inhibit branch migration of the junction and to probe the possible sequence dependence of these four-arm DNA structures. Overlap of resonances in homonuclear two-dimensional nuclear Overhauser enhancement (NOE) spectra necessitates the use of a multipathway approach for obtaining sequence-specific assignments, wherein all possible NOE connectivities are analyzed in parallel. Using this strategy, 1H resonance assignments were obtained for virtually all nonlabile base protons and C1', C2', and C3' sugar protons. Several unambiguous cross-arm NOE connectivities were identified, directly establishing the stacking arrangements of each contiguous (two-arm) helical domain. The distribution of the two possible stacking isomers is distinctly different for the two junctions studied, thereby indicating that the relative stability of the isomers is dependent on the sequence at the junction.

High level expression in Escherichia coli and characterization of the EF-hand calcium-binding protein caltractin

Caltractin is a member of the calmodulin superfamily of Ca(2+)-binding proteins that was originally cloned at the DNA level from the unicellular green alga Chlamydomonas reinhardtii. Human and mouse homologs to algal caltractin have been recently characterized. In the studies reported here, recombinant Chlamydomonas caltractin was expressed at high levels in Escherichia coli and purified to homogeneity. The use of the ompT-host BL21 proved critical for obtaining high yields of homogeneous full-length protein. Growth and purification protocols were optimized to allow reproducible and efficient production of tens of milligrams of pure protein from 1-liter cultures. Caltractin has a distinct UV spectrum which is largely dominated by the fine structure due to the 9 Phe residues. Unlike other members of the same protein family, the UV and the CD spectra do not change upon addition of Ca2+ to the apoprotein. However, the 1H NMR spectrum shows distinct changes upon Ca2+ binding, which are indicative of structural and/or dynamic changes largely reminiscent of other members of the calmodulin superfamily. Ca2+ binding measurements demonstrated the binding of four Ca2+ ions to caltractin with two higher affinity (Kd = 1.2 x 10(-6) M) and two lower affinity (Kd = 1.6 x 10(-4) M) sites. Caltractin is highly stable in both the apo- and the Ca(2+)-loaded states. The unusual stability of apocaltractin makes this protein highly suited for structural studies by multidimensional NMR aimed at understanding the structural and dynamic consequences of Ca2+ binding, and the molecular basis of Ca2+ signal transduction.

Signal transduction versus buffering activity in Ca(2+)-binding proteins

The three-dimensional structure of calbindin D9k in the absence of Ca2+ has been determined using NMR spectroscopy in solution, allowing the first direct analysis of the consequences of Ca2+ binding for a member of the calmodulin superfamily of proteins. The overall response in calbindin D9k is much attenuated relative to the current model for calmodulin and troponin C. These results demonstrate a novel mechanism for modulating the conformational response to Ca(2+)-binding in calmodulin superfamily proteins and provide insights into how their Ca(2+)-binding domains can be fine-tuned to remain essentially intact or respond strongly to ion binding, in relation to their functional requirements.

Effects of ion binding on the backbone dynamics of calbindin D9k determined by 15N NMR relaxation

The backbone dynamics of apo- and (Cd2+)1-calbindin D9k have been characterized by 15N nuclear magnetic resonance spectroscopy. Spin-lattice and spin-spin relaxation rate constants and steady-state [1H]-15N nuclear Overhauser effects were measured at a magnetic field strength of 11.74 T by two-dimensional, proton-detected heteronuclear NMR experiments using 15N-enriched samples. The relaxation parameters were analyzed using a model-free formalism that characterizes the dynamics of the N-H bond vectors in terms of generalized order parameters and effective correlation times. The data for the apo and (Cd2+)1 states were compared to those for the (Ca2+)2 state [Kördel, J., Skelton, N. J., Akke, M., Palmer, A. G., & Chazin, W. J. (1992) Biochemistry 31, 4856-4866] to ascertain the effects on ion ligation on the backbone dynamics of calbindin D9k. The two binding loops respond differently to ligation by metal ions: high-frequency (10(9)-10(12) s-1) fluctuations of the N-terminal ion-binding loop are not affected by ion binding, whereas residues G57, D58, G59, and E60 in the C-terminal ion-binding loop have significantly lower order parameters in the apo state than in the metal-bound states. The dynamical responses of the four helices to binding of ions are much smaller than that for the C-terminal binding loop, with the strongest effect on helix III, which is located between the linker loop and binding site II. Significant fluctuations on slower time scales also were detected in the unoccupied N-terminal ion-binding loop of the apo and (Cd2+)1 states; the apparent rates were greater for the (Cd2+)1 state. These results on the dynamical response to ion binding in calbindin D9k provide insights into the molecular details of the binding process and qualitative evidence for entropic contributions to the cooperative phenomenon of calcium binding for the pathway in which the ion binds first in the C-terminal site.

High-resolution structure of calcium-loaded calbindin D9k

The three-dimensional solution structure of calcium-loaded calbindin D9k has been determined using experimental constraints obtained from nuclear magnetic resonance spectroscopy. A total of 1176 constraints (16 per residue overall, 32 per residue for the core residues) was used for the final refinement, including 1002 distance and 174 dihedral angle constraints. In addition, 23 hydrogen bond constraints were used for the generation of initial structures. Stereospecific assignments were made for 37 of 61 (61%) prochiral methylene protons and the methyl groups of all three valine residues and five out of 12 leucine residues. These constraints were used as input for a series of calculations of three-dimensional structures using a combination of distance geometry and restrained molecular dynamics. The 33 best structures selected for further analysis have no distance constraint violations greater than 0.3 A and good local geometries as reflected by low total energies (< or = -1014 kcal/mol in the AMBER 4.0 force field). The core of the protein consists of four well-defined helices with root-mean-square deviations from the average of 0.45 A for the N, C alpha and C' backbone atoms. These helices are packed in an antiparallel fashion to form two helix-loop-helix calcium-binding motifs, termed EF-hands. The two EF-hands are joined at one end by a ten-residue linker segment, and at the other by a short beta-type interaction between the two calcium-binding loops. Overall, the average solution structure of calbindin D9k is very similar to the crystal structure, with a pairwise root-mean-square deviation of 0.85 A for the N, C alpha and C' backbone atoms of the four helices. The differences that are observed between the solution and the crystal structures are attributed to specific crystal contacts, increased side-chain flexibility in solution, or artifacts arising from molecular dynamics refinement of the solution structures in vacuo.

Two-dimensional 1H nuclear magnetic resonance studies of the half-saturated (Ca2+)1 state of calbindin D9k. Further implications for the molecular basis of cooperative Ca2+ binding

Calbindin D9k exhibits cooperative binding of two calcium ions, hence study of the half-saturated states of the protein is critical to understanding the binding process. However, the half-saturated states are not significantly populated under equilibrium conditions. To circumvent this problem, an absolutely conserved glutamic acid residue in the C-terminal binding site (site II) has been mutated to glutamine (E65Q), causing a substantial reduction in calcium affinity and permitting detailed two-dimensional 1H NMR analysis of calbindin D9k with a calcium ion bound only in the N-terminal EF-hand. Complete 1H resonance assignments have been obtained for (Ca2+)1 E65Q, as well as near complete assignments for the apo and (Ca2+)2 states. A value of 1.1(+/- 0.2) x 10(3) M-1 has been determined for the calcium binding constant in site II, from an analysis of the chemical shift changes in response to titration with calcium. The elements of secondary structure and global folding patterns were identified from nuclear Overhauser effects, backbone spin-spin coupling constants and the exchange rates of backbone amide protons. Although the mutation has only very small effects on the secondary structure and global fold of the protein, it so drastically lowers affinity for Ca2+ in the C-terminal site that (Ca2+)2 E65Q does not correspond to a standard (Ca2+)2 state. From the analysis of the half-saturated state, it is apparent that some reorganization of the structure and changes in the internal dynamics of calbindin D9k does occur for each step of the apo-->(Ca2+)1(I)-->(Ca2+)2 binding pathway. When the first ion is bound to the N-terminal EF-hand, that half of the molecule adopts a conformation and dynamic state similar to the fully calcium-loaded protein state, whereas only minor changes occur in the C-terminal EF-hand. It is only upon binding of the second calcium ion that the C-terminal EF-hand switches over to the fully calcium-loaded state. Together with the results from our earlier study of the apo-->(Ca2+)1(II)-->(Ca2+)2 binding pathway, these findings indicate that changes in protein conformation and dynamics associated with Ca2+ binding contribute to the observed positive cooperativity, and that the molecular details of the cooperative binding events are different for the two binding pathways.

Two-dimensional 1H NMR studies of 32-base-pair synthetic immobile Holliday junctions: complete assignments of the labile protons and identification of the base-pairing scheme

The labile protons of two 32-base-pair, four-arm models of immobile Holliday junctions have been studied by two-dimensional 1H nuclear magnetic resonance (NMR) spectroscopy. Overlap of resonances in the imino-imino region of two-dimensional nuclear Overhauser enhancement (NOE) spectra necessitates the use of a multi-pathway approach for obtaining sequence-specific assignments wherein all possible NOE connectivities to the labile protons are utilized, including those from the 2H of adenine, 5CH3 of thymine, and 5H of cytosine. Resonance assignments are obtained for all slowly exchanging imino and cytosine amino protons. Base-pairing up to and including the junction point is found in all four arms of both Holliday junctions. Several cross-arm NOE connectivities are identified and can be used to infer the geometry of the helical stacking domains. The two Holliday junctions studied, which differ only by the exchange of two base pairs at the branch point, appear to have opposite arm stacking geometries. These assignments form an important part of the critical background for detailed NMR analysis of Holliday junction structure and dynamics.

Nuclear magnetic resonance studies of the internal dynamics in Apo, (Cd2+)1 and (Ca2+)2 calbindin D9k. The rates of amide proton exchange with solvent

The backbone dynamics of the EF-hand Ca(2+)-binding protein, calbindin D9k, has been investigated in the apo, (Cd2+)1 and (Ca2+)2 states by measuring the rate constants for amide proton exchange with solvent. 15N-1H correlation spectroscopy was utilized to follow direct 1H-->2H exchange of the slowly exchanging amide protons and to follow indirect proton exchange via saturation transfer from water to the rapidly exchanging amide protons. Plots of experimental rate constants versus intrinsic rate constants have been analyzed to give qualitative insight into the opening modes of the protein that lead to exchange. These results have been interpreted within the context of a progressive unfolding model, wherein hydrophobic interactions and metal chelation serve to anchor portions of the protein, thereby damping fluctuations and retarding amide proton exchange. The addition of Ca2+ or Cd2+ was found to retard the exchange of many amide protons observed to be in hydrogen-bonding environments in the crystal structure of the (Ca2+)2 state, but not of those amide protons that were not involved in hydrogen bonds. The largest changes in rate constant occur for residues in the ion-binding loops, with substantial effects also found for the adjacent residues in helices I, II and III, but not helix IV. The results are consistent with a reorganization of the hydrogen-bonding networks in the metal ion-binding loops, accompanied by a change in the conformation of helix IV, as metal ions are chelated. Further analysis of the results obtained for the three states of metal occupancy provides insight into the nature of the changes in conformational fluctuations induced by ion binding.

15N NMR assignments and chemical shift analysis of uniformly labeled 15N calbindin D9k in the apo, (Cd2+)1 and (Ca2+)2 states

15N has been uniformly incorporated into the EF-hand Ca(2+)-binding protein calbindin D9k so that heteronuclear experiments can be used to further characterize the structure and dynamics of the apo, (Cd2+)1 and (Ca2+)2 states of the protein. The 15N NMR resonances were assigned by 2D 15N-resolved 1H experiments, which also allowed the identification of a number of sequential and medium-range 1H-1H contacts that are obscured by chemical shift degeneracy in homonuclear experiments. The 15N chemical shifts are analyzed with respect to correlations with protein secondary structure. In addition, the changes in 15N chemical shift found for the apo----(Cd2+)1----(Ca2+)2 binding sequence confirm that the effects on the protein are mainly associated with chelation of the first ion.

Backbone dynamics of calcium-loaded calbindin D9k studied by two-dimensional proton-detected 15N NMR spectroscopy

Backbone dynamics of calcium-loaded calbindin D9k have been investigated by two-dimensional proton-detected heteronuclear nuclear magnetic resonance spectroscopy, using a uniformly 15N enriched protein sample. Spin-lattice relaxation rate constants, spin-spin relaxation rate constants, and steady-state [1H]-15N nuclear Overhauser effects were determined for 71 of the 72 backbone amide 15N nuclei. The relaxation parameters were analyzed using a model-free formalism that incorporates the overall rotational correlation time of the molecule, and a generalized order parameter (S2) and an effective internal correlation time for each amide group. Calbindin D9k contains two helix-loop-helix motifs joined by a linker loop at one end of the protein and a beta-type interaction between the two calcium-binding loops at the other end. The amplitude of motions for the calcium-binding loops and the helices are similar, as judged from the average S2 values of 0.83 +/- 0.05 and 0.85 +/- 0.04, respectively. The linker region joining the two calcium-binding subdomains of the molecule has a significantly higher flexibility, as indicated by a substantially lower average S2 value of 0.59 +/- 0.23. For residues in the linker loop and at the C-terminus, the order parameter is further decomposed into separate order parameters for motional processes on two distinct time scales. The effective correlation times are significantly longer for helices I and IV than for helices II and III or for the calcium-binding loops. Residue by residue comparisons reveal correlations of the order parameters with both the crystallographic B-factors and amide proton exchange rates, despite vast differences in the time scales to which these properties are sensitive. The order parameters are also utilized to distinguish regions of the NMR-derived three-dimensional structure of calbindin D9k that are poorly defined due to inherently high flexibility, from poorly defined regions with average flexibility but a low density of structural constraints.

Two-dimensional NMR studies of d(GGTTAATGCGGT).d(ACCGCATTAACC) complexed with the minor groove binding drug SN-6999

The dodecadeoxynucleotide duplex d(GGTTAATGCGGT).d(ACCGCATTAACC) and its 1:1 complex with the minor groove binding drug SN-6999 have been prepared and studied by two-dimensional 1H nuclear magnetic resonance spectroscopy. Complete sequence-specific assignments have been obtained for the free duplex by standard methods. The line widths of the resonances in the complex are greater than those observed for the free duplex, which complicates the assignment process. Extensive use of two-quantum spectroscopy was required to determine the scalar correlations for identifying all of the base proton and most of the 1'H-2'H-2''H spin subsystems for the complex. This permitted unambiguous sequence-specific resonance assignments for the complex, which provides the necessary background for a detailed comparison of the structure of the duplex, with and without bound drug. A series of intermolecular NOEs between drug and DNA were identified, providing sufficient structural constraints to position the drug in the minor groove of the duplex. However, the combination of NOEs observed can only be rationalized by a model wherein the drug binds in the minor groove of the DNA in both orientations relative to the long helix axis and exchanges rapidly between the two orientations. The drug binds primarily in the segment of five consecutive dA-dT base pairs d(T3T4A5A6T7).d(A18T19T20A21A22), but surprisingly strong interactions are found to extend one residue in the 3' direction along each strand to G8 and C23. The observation of intermolecular contacts to residues neighboring the AT-rich region demonstrates that the stabilization of the bis(quaternary ammonium) heterocycle family of AT-specific, minor groove binding drugs is not based exclusively on interactions with dA-dT base pairs.

Three-dimensional solution structure of Ca(2+)-loaded porcine calbindin D9k determined by nuclear magnetic resonance spectroscopy

The three-dimensional solution structure of native, intact porcine calbindin D9k has been determined by distance geometry and restrained molecular dynamics calculations using distance and dihedral angle constraints obtained from 1H NMR spectroscopy. The protein has a well-defined global fold consisting of four helices oriented in a pairwise antiparallel manner such that two pairs of helix-loop-helix motifs (EF-hands) are joined by a linker segment. The two EF-hands are further coupled through a short beta-type interaction between the two Ca(2+)-binding loops. Overall, the structure is very similar to that of the highly homologous native, minor A form of bovine calbindin D9k determined by X-ray crystallography [Szebenyi, D. M. E., & Moffat, K. (1986) J. Biol. Chem. 261, 8761-8776]. A model structure built from the bovine calbindin D9k crystal structure shows several deviations larger than 2 A from the experimental distance constraints for the porcine protein. These structural differences are efficiently removed by subjecting the model structure to the experimental distance and dihedral angle constraints in a restrained molecular dynamics protocol, thereby generating a model that is very similar to the refined distance geometry derived structures. The N-terminal residues of the intact protein that are absent in the minor A form appear to be highly flexible and do not influence the structure of other regions of the protein. This result is important because it validates the conclusions drawn from the wide range of studies that have been carried out on minor A forms rather than the intact calbindin D9k.

Conformational studies of the duplex d-(CCAAAAATTTCC).d-(GGAAATTTTTGG) containing a (dA)5 tract using two-dimensional 1H-n.m.r. spectroscopy

A variety of two-dimensional nuclear magnetic resonance (n.m.r.) techniques were utilized for structural studies of the dodecadeoxy-nucleotide duplex d-(CCAAAAATTTCC).d-(GGAAATTTTTGG). Overall, the duplex was found to adopt a conformation of the B-DNA type. Direct measurement of 1H-1H coupling constants of the deoxyribose rings and comparison of nuclear Overhauser effects (NOEs) involving characteristic intra- and inter-nucleotide proton-proton connectivities indicate that the sugar rings are all in predominantly C2'-endo/S-type conformation. The appearance of NOEs between adenine 2H and sugar l'H resonances indicates a narrowing of the minor groove at the 3' end of the (dA)n tracts. Several lines of evidence indicate that a conformational change in the DNA duplex is introduced by the (dA)5 tract, in support of the junction model for curved DNA. The results are compared with previous structural studies of the evidence for curvature in duplex DNA containing (dA)n tracts.