Distortion of the three-dimensional structure of the vnd/NK-2 homeodomain bound to DNA induced by an embryonically lethal A35T point mutation

The three-dimensional solution structure obtained by NMR of the A35T mutant vnd/NK-2 homeodomain bound to the vnd/NK-2 consensus 16 bp DNA sequence was determined. This mutation to threonine from alanine in position 35 in helix II of the vnd/NK-2 homeodomain is associated with early embryonic lethality in Drosophila melanogaster. Although the unbound mutant protein is not structured, in the DNA-bound state it adopts the three-helix fold characteristic of all known homeodomains, but with alterations relative to the structure of the wild-type analogue. These structural modifications occur, and are accompanied by a 50-fold reduction in the DNA binding affinity, even though most of the protein-DNA interactions originally seen for the wild-type homeodomain are found likewise in the threonine analogue. Alterations include torsional angle changes in the loop between helix I and helix II, and in the turn between helix II and helix III, as well as in a distortion of the usual antiparallel orientation of helix I with respect to helix II. The alteration of the position of leucine 40 in the A35T mutant is proposed to explain the observed 1.27 ppm upfield shift of the corresponding amide proton resonance relative to the value observed for the wild-type analogue. A detailed comparison of the structures of the mutant A35T and wild-type vnd/NK-2 homeodomains bound to the cognate DNA is presented. The consequences of the structural alteration of the DNA-bound A35T mutant vnd/NK-2 protein may constitute the basis of the observed early embryonic lethality.     

The vnd/NK-2 homeodomain: thermodynamics of reversible unfolding and DNA binding for wild-type and with residue replacements H52R and H52R/T56W in helix III

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.     

The three-dimensional structure of the vnd/NK-2 homeodomain-DNA complex by NMR spectroscopy.

The three-dimensional solution structure obtained by NMR of the complex formed between the uniformly singly15N and doubly13C/15N-labeled vnd/NK-2 homeodomain and its consensus 16 base-pair DNA binding sequence was determined. This work was carried out using the accepted repertoire of experiments augmented with a novel implementation of the water flipback technique to enhance signals from exchangeable amide protons. The results using this new technique confirm the existence of hydrogen bonding between the invariant Asn51 and the second adenine of the DNA binding sequence, as seen in crystal structures of other homeodomain-DNA complexes, but never before detected by NMR. Hydrogen bonding by Arg5 and Lys3 in the minor groove of the DNA appears to be responsible for two unusually upfield-shifted ribose H1' resonances. The DNA duplex is nearly straight and its structure is primarily that of B -DNA. A detailed comparison is presented for all available homeodomain-DNA structures including the vnd/NK-2 DNA complex, which demonstrates that homology is maintained in the protein structure, whereas for the orientation of the homeodomain relative to DNA, small but significant variations are observed. Interactions are described involving certain residues in specific positions of the homeodomain, namely Leu7, Thr41, and Gln50 of vnd/NK-2, where single amino acid residue mutations lead to dramatic developmental alterations. The availability of our previously determined three- dimensional structure of the vnd/NK-2 homeodomain in the absence of DNA allows us to assess structural changes in the homeodomain induced by DNA binding.     

Quantitative measurement of water diffusion lifetimes at a protein/DNA interface by NMR

Hydration site lifetimes of slowly diffusing water molecules at the protein/DNA interface of the vnd/NK-2 homeodomain DNA complex were determined using novel three-dimensional NMR techniques. The lifetimes were calculated using the ratios of ROE and NOE cross-relaxation rates between the water and the protein backbone and side chain amides. This calculation of the lifetimes is based on a model of the spectral density function of the water-protein interaction consisting of three timescales of motion: fast vibrational/rotational motion, diffusion into/out of the hydration site, and overall macromolecular tumbling. The lifetimes measured ranged from approximately 400 ps to more than 5 ns, and nearly all the slowly diffusing water molecules detected lie at the protein/DNA interface. A quantitative analysis of relayed water cross-relaxation indicated that even at very short mixing times, 5 ms for ROESY and 12 ms for NOESY, relay of magnetization can make a small but detectable contribution to the measured rates. The temperature dependences of the NOE rates were measured to help discriminate direct dipolar cross-relaxation from chemical exchange. Comparison with several X-ray structures of homeodomain/DNA complexes reveals a strong correspondence between water molecules in conserved locations and the slowly diffusing water molecules detected by NMR. A homology model based on the X-ray structures was created to visualize the conserved water molecules detected at the vnd/NK-2 homeodomain DNA interface. Two chains of water molecules are seen at the right and left sides of the major groove, adjacent to the third helix of the homeodomain. Two water-mediated hydrogen bond bridges spanning the protein/DNA interface are present in the model, one between the backbone of Phe8 and a DNA phosphate, and one between the side chain of Asn51 and a DNA phosphate. The hydrogen bond bridge between Asn51 and the DNA might be especially important since the DNA contact made by the invariant Asn51 residue, seen in all known homeodomain/DNA structures, is critical for binding affinity and specificity.     

Cardiac-specific Nkx2.5 homeodomain: conformational stability and specific DNA binding of Nkx2.5(C56S)

The cardiac-specific Nkx2.5 homeodomain has been expressed as a 79-residue protein with the oxidizable Cys(56) replaced with Ser. The Nkx2.5 or Nkx2.5(C56S) homeodomain is 73% identical in sequence to and has the same NMR structure as the vnd (ventral nervous system defective)/NK-2 homeodomain of Drosophila when bound to the same specific DNA. The thermal unfolding of Nkx2.5(C56S) at pH 6.0 or 7.4 is a reversible, two-state process with unit cooperativity, as measured by differential scanning calorimetry (DSC) and far-UV circular dichroism. Adding 100 mM NaCl to Nkx2.5(C56S) at pH 7.4 increases T(m) from 44 to 54 +/- 0.2 degrees C and DeltaH from 34 to 45 +/- 2 kcal/mol (giving a DeltaC(p) of approximately 1.2 kcal K(-)(1) mol(-)(1) for homeodomain unfolding). DSC profiles of Nkx2.5 indicate fluctuating nativelike structures at <37 degrees C. Titrations of specific 18 bp DNA with Nkx2.5(C56S) in buffer at pH 7.4 with 100 mM NaCl yield binding constants of 2-6 x 10(8) M(-)(1) from 10 to 37 degrees C and a stoichiometry of 1:1 for homeodomain binding DNA, using isothermal titration calorimetry. The DNA binding reaction of Nkx2.5 is enthalpically controlled, and the temperature dependence of DeltaH gives a DeltaC(p) of -0.18 +/- 0.01 kcal K(-)(1) mol(-)(1). This corresponds to 648 +/- 36 A(2) of buried apolar surface upon Nkx2.5(C56S) binding duplex B-DNA. Thermodynamic parameters differ for Nkx2.5 and vnd/NK-2 homeodomains binding specific DNA. Unbound NK-2 is more flexible than Nkx2.5.     

Smoluchowski dynamics of the vnd/NK-2 homeodomain from Drosophila melanogaster: second-order maximum correlation approximation

The mode coupling diffusion theory is applied to the derivation of local dynamics in proteins in solution. The rotational dynamics of the bonds along the protein sequence are calculated and compared to the experimentally measured nmr (15)N spin-lattice relaxation time T(1), at 36.5, 60.8, and 81.1 MHz of the vnd/NK-2 homeodomain from Drosophila melanogaster. The starting point for the calculations is the experimental three-dimensional solution structure of the homeodomain determined by multidimensional nmr spectroscopy. The higher order mode-coupling computations are compared also with the recently published first-order approximation calculations. The more accurate calculations improve substantially the first-order ORZLD calculations and show that the role of the strength of the hydrodynamic interactions becomes crucial to fix the order of magnitude of the rotational dynanics for these very compact molecules characterized by partial screening of the internal atoms to water. However, the relative mobility of the bonds along the sequence and the differential fluctuations depend only weakly on the hydrodynamic strength but strongly on the geometry of the three-dimensional structure and on the statistics incorporated into the theory. Both rigid and fluctuating dynamic models are examined, with fluctuations evaluated using molecular dynamics simulations. The comparison with nmr data shows that mode coupling diffusion accounts for the T(1) relaxation pattern at low frequency where the rotational tumbling dominates. An important contribution of internal motions in the nanosecond time scale is seen at high frequencies and is discussed in terms of diffusive concepts.     

Functional consequences of preorganized helical structure in the intrinsically disordered cell-cycle inhibitor p27(Kip1).

p27(Kip1) contributes to cell-cycle regulation by inhibiting cyclin-dependent kinase (Cdk) activity. The p27 Cdk-inhibition domain has an ordered conformation comprising an alpha-helix, a 3(10) helix, and beta-structure when bound to cyclin A-Cdk2. In contrast, the unbound p27 Cdk-inhibition domain is intrinsically disordered (natively unfolded) as shown by circular dichroism spectroscopy, lack of chemical-shift dispersion, and negative heteronuclear nuclear Overhauser effects. The intrinsic disorder is not due to the excision of the Cdk-inhibition domain from p27, since circular dichroism spectra of the full-length protein are also indicative of a largely unfolded protein. Both the inhibition domain and full-length p27 are active as cyclin A-Cdk2 inhibitors. Using circular dichroism and proline mutagenesis, we demonstrate that the unbound p27 Cdk-inhibition domain is not completely unfolded. The domain contains marginally stable helical structure that presages the alpha-helix, but not the 3(10) helix, adopted upon binding cyclin A-Cdk2. Increasing or reducing the stability of the partially preformed alpha-helix in the isolated p27 domain with alanine or proline substitutions did not affect formation of the p27-inhibited cyclin A-Cdk2 complex in energetic terms. However, stabilization of the helix with alanine hindered kinetically the formation of the inhibited complex, suggesting that p27 derives a kinetic advantage from intrinsic structural disorder.     

Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease

This paper describes the use of novel two-dimensional nuclear magnetic resonance (NMR) pulse sequences to provide insight into protein dynamics. The sequences developed permit the measurement of the relaxation properties of individual nuclei in macromolecules, thereby providing a powerful experimental approach to the study of local protein mobility. For isotopically labeled macromolecules, the sequences enable measurements of heteronuclear nuclear Overhauser effects (NOE) and spin-lattice (T1) and spin-spin (T2) 15N or 13C relaxation times with a sensitivity similar to those of many homonuclear 1H experiments. Because T1 values and heteronuclear NOEs are sensitive to high-frequency motions (10(8)-10(12) s-1) while T2 values are also a function of much slower processes, it is possible to explore dynamic events occurring over a large time scale. We have applied these techniques to investigate the backbone dynamics of the protein staphylococcal nuclease (S. Nase) complexed with thymidine 3',5'-bisphosphate (pdTp) and Ca2+ and labeled uniformly with 15N. T1, T2, and NOE values were obtained for over 100 assigned backbone amide nitrogens in the protein. Values of the order parameter (S), characterizing the extent of rapid 1H-15N bond motions, have been determined. These results suggest that there is no correlation between these rapid small amplitude motions and secondary structure for S. Nase. In contrast, 15N line widths suggest a possible correlation between secondary structure and motions on the millisecond time scale. In particular, the loop region between residues 42 and 56 appears to be considerably more flexible on this slow time scale than the rest of the protein.     

Backbone dynamics of the glucocorticoid receptor DNA-binding domain.

The extent of rapid (picosecond) backbone motions within the glucocorticoid receptor DNA-binding domain (GR DBD) has been investigated using proton-detected heteronuclear NMR spectroscopy on uniformly 15N-labeled protein fragments containing the GR DBD. Sequence-specific 15N resonance assignments, based on two- and three-dimensional heteronuclear NMR spectra, are reported for 65 of 69 backbone amides within the segment C440-A509 of the rat GR in a protein fragment containing a total of 82 residues (MW = 9200). Individual backbone 15N spin-lattice relaxation times (T1), rotating-frame spin-lattice relaxation times (T1 rho), and steady-state (1H)-15N nuclear Overhauser effects (NOEs) have been measured at 11.74 T for a majority of the backbone amide nitrogens within the segment C440-N506. T1 relaxation times and NOEs are interpreted in terms of a generalized order parameter (S2) and an effective correlation time (tau e) characterizing internal motions in each backbone amide using an optimized value for the correlation time for isotropic rotational motions of the protein (tau R = 6.3 ns). Average S2 order parameters are found to be similar (approximately 0.86 +/- 0.07) for various functional domains of the DBD. Qualitative inspection as well as quantitative analysis of the relaxation and NOE data suggests that the picosecond flexibility of the DBD backbone is limited and uniform over the entire protein, with the possible exception of residues S448-H451 of the first zinc domain and a few residues for which relaxation and NOE parameters were not obtained. in particular, we find no evidence for extensive rapid backbone motions within the second zinc domain. Our results therefore suggest that the second zinc domain is not disordered in the uncomplexed state of DBD, although the possibility of slowly exchanging (ordered) conformational states cannot be excluded in the present analysis.     

(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles

We report the backbone dynamics of monomeric phospholamban in dodecylphosphocholine micelles using (1)H/(15)N heteronuclear NMR spectroscopy. Phospholamban is a 52-amino acid membrane protein that regulates Ca-ATPase in cardiac muscle. Phospholamban comprises three structural domains: a transmembrane domain from residues 22 to 52, a connecting loop from 17 to 21, and a cytoplasmic domain from 1 to 16 that is organized in an "L"-shaped structure where the transmembrane and the cytoplasmic domain form an angle of approximately 80 degrees (Zamoon et al., 2003; Mascioni et al., 2002). T(1), T(2), and (1)H/(15)N nuclear Overhauser effect values measured for the amide backbone resonances were interpreted using the model-free approach of Lipari and Szabo. The results point to the existence of four dynamic domains, revealing the overall plasticity of the cytoplasmic helix, the flexible loop, and part of the transmembrane domain (residues 22-30). In addition, using Carr-Purcell-Meiboom-Gill-based experiments, we have characterized phospholamban dynamics in the micros-ms timescale. We found that the majority of the residues in the cytoplasmic domain, the flexible loop, and the first ten residues of the transmembrane domain undergo dynamics in the micros-ms range, whereas minimal dynamics were detected for the transmembrane domain. Hydrogen/deuterium exchange factors measured at different temperatures support the existence of slow motion in both the loop and the cytoplasmic helix. We propose that these dynamic properties are critical factors in the biomolecular recognition of phospholamban by Ca-ATPase and other interacting proteins such as protein kinase A and protein phosphatase 1.     

Serine 16 phosphorylation induces an order-to-disorder transition in monomeric phospholamban

Phospholamban (PLB) is a 52 amino acid membrane-endogenous regulator of the sarco(endo)plasmic calcium adenosinetriphosphatase (SERCA) in cardiac muscle. PLB's phosphorylation and dephosphorylation at S16 modulate its regulatory effect on SERCA by an undetermined mechanism. In this paper, we use multidimensional (1)H/(15)N solution NMR methods to establish the structural and dynamics basis for PLB's control of SERCA upon S16 phosphorylation. For our studies, we use a monomeric, fully active mutant of PLB, where C36, C41, and C46 have been mutated to A36, F41, and A46, respectively. Our data show that phosphorylation disrupts the "L-shaped" structure of monomeric PLB, causing significant unwinding of both the cytoplasmic helix (domain Ia) and the short loop (residues 17-21) connecting this domain to the transmembrane helix (domains Ib and II). Concomitant with this conformational transition, we also find pronounced changes in both the pico- to nanosecond and the micro- to millisecond time scale dynamics. The (1)H/(15)N heteronuclear NOE values for residues 1-25 are significantly lower than those of unphosphorylated PLB, with slightly lower NOE values in the transmembrane domain, reflecting less restricted motion throughout the whole protein. These data are supported by the faster spin-lattice relaxation rates (R(1)) present in both the cytoplasmic and loop regions and by the enhanced spin-spin transverse relaxation rates (R(2)) observed in the transmembrane domain. These results demonstrate that while S16 phosphorylation induces a localized structural transition, changes in PLB's backbone dynamics are propagated throughout the protein backbone. We propose that the regulatory mechanism of PLB phosphorylation involves an order-to-disorder transition, resulting in a decrease in the PLB inhibition of SERCA.     

The stability of the cytochrome c scaffold as revealed by NMR spectroscopy

NMR spectroscopy was used to study the effect of guanidinium chloride on the unfolding of horse heart and yeast iso-1 cytochrome c under mild alkaline conditions. The structural changes on the horse heart protein were detected through NOESY (Nuclear Overhauser Effect SpectroscopY) experiments whereas (15)N-(1)H heteronuclear NMR was used to monitor the behavior of the yeast protein. The latter represents the first characterization through (15)N-(1)H heteronuclear NMR spectroscopy of the guanidinium chloride induced unfolding of mitochondrial cytochrome c. The presence of denaturants decreases the temperature at which the native Met80 axial ligand is displaced from the iron center under the present mild alkaline conditions. The process can be described in terms of protein fragments behaving as unfolding units of different stability. The comparison between the two proteins indicates that the loop+helix connecting the proximal and distal sites, as well as the long Met80-containing loop immediately after a short helix, are structural characteristics of mitochondrial cytochrome c that appear to be responsible for the Met80-iron(III) bond fragility.     

Effects of varying the local propensity to form secondary structure on the stability and folding kinetics of a rapid folding mixed alpha/beta protein: characterization of a truncation mutant of the N-terminal domain of the ribosomal protein L9.

The N-terminal domain of the ribosomal protein L9 forms a split betaalphabeta structure with a long C-terminal helix. The folding transitions of a 56 residue version of this protein have previously been characterized, here we report the results of a study of a truncation mutant corresponding to residues 1-51. The 51 residue protein adopts the same fold as the 56 residue protein as judged by CD and two-dimensional NMR, but it is less stable as judged by chemical and thermal denaturation experiments. Studies with synthetic peptides demonstrate that the C-terminal helix of the 51 residue version has very little propensity to fold in isolation in contrast to the C-terminal helix of the 56 residue variant. The folding rates of the two proteins, as measured by stopped-flow fluorescence, are essentially identical, indicating that formation of local structure in the C-terminal helix is not involved in the rate-limiting step of folding.