Conformational changes of the BS2 operator DNA upon complex formation with the Antennapedia homeodomain studied by NMR with 13C/15N-labeled DNA

Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

Slow Rotational Dynamics of Cytosine NH2 Groups in Double-Stranded DNA

Aromatic NH₂ groups are essential as hydrogen-bond donors in secondary structures of DNA and RNA. Although rapid rotations of NH₂ groups of adenine and guanine bases were previously characterized, there has been a lack of quantitative information about slow rotations of cytosine NH₂ groups in Watson-Crick base pairs. In this study, using an NMR method we had recently developed, we determined the kinetic rate constants and energy barriers for cytosine NH₂ rotations in a 15-base-pair DNA duplex. Our data show that the rotational dynamics of cytosine NH₂ groups depend on local environments. Qualitative correlation between the ranges of ¹⁵N chemical shifts and rotational time scales for various NH₂ groups of nucleic acids and proteins illuminates a relationship between the partial double-bond character of the C-N bond and the time scale for NH₂ rotations.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

Internal Motions of Basic Side Chains of the Antennapedia Homeodomain in the Free and DNA-Bound States

Basic side chains play crucial roles in protein-DNA interactions. In this study, using NMR spectroscopy, we investigated the dynamics of Arg and Lys side chains of the fruit fly Antennapedia homeodomain in the free state and in the complex with target DNA. We measured 15N relaxation for Arg and Lys side chains at two magnetic fields, from which generalized order parameters for the cationic groups were determined. Mobility of the R5 side chain, which makes hydrogen bonds with a thymine base in the DNA minor groove, was greatly dampened. Several Lys and Arg side chains that form intermolecular ion pairs with DNA phosphates were found to retain high mobility with the order parameter being <0.6 in the DNA-bound state. Interestingly, some of the interfacial cationic groups in the complex were more mobile than in the free protein. The retained or enhanced mobility of the Arg and Lys side chains in the complex should mitigate the overall loss of conformational entropy in the protein-DNA association and allow dynamic molecular recognition.

Entropic Enhancement of Protein-DNA Affinity by Oxygen-to-Sulfur Substitution in DNA Phosphate

Dithioation of DNA phosphate is known to enhance binding affinities, at least for some proteins. We mechanistically characterized this phenomenon for the Antennapedia homeodomain-DNA complex by integrated use of fluorescence, isothermal titration calorimetry, NMR spectroscopy, and x-ray crystallography. By fluorescence and isothermal titration calorimetry, we found that this affinity enhancement is entropy driven. By NMR, we investigated the ionic hydrogen bonds and internal motions of lysine side-chain NH3(+) groups involved in ion pairs with DNA. By x-ray crystallography, we compared the structures of the complexes with and without dithioation of the phosphate. Our NMR and x-ray data show that the lysine side chain in contact with the DNA phosphate becomes more dynamic upon dithioation. Our thermodynamic, structural, and dynamic investigations collectively suggest that the affinity enhancement by the oxygen-to-sulfur substitution in DNA phosphate is largely due to an entropic gain arising from mobilization of the intermolecular ion pair at the protein-DNA interface.

GFT projection NMR for efficient (1)H/ (13)C sugar spin system identification in nucleic acids

A newly implemented G-matrix Fourier transform (GFT) (4,3)D HC(C)CH experiment is presented in conjunction with (4,3)D HCCH to efficiently identify (1)H/(13)C sugar spin systems in (13)C labeled nucleic acids. This experiment enables rapid collection of highly resolved relay 4D HC(C)CH spectral information, that is, shift correlations of (13)C-(1)H groups separated by two carbon bonds. For RNA, (4,3)D HC(C)CH takes advantage of the comparably favorable 1'- and 3'-CH signal dispersion for complete spin system identification including 5'-CH. The (4,3)D HC(C)CH/HCCH based strategy is exemplified for the 30-nucleotide 3'-untranslated region of the pre-mRNA of human U1A protein.

Solution structure of the K50 class homeodomain PITX2 bound to DNA and implications for mutations that cause Rieger syndrome

We have determined the solution structure of a complex containing the K50 class homeodomain Pituitary homeobox protein 2 (PITX2) bound to its consensus DNA site (TAATCC). Previous studies have suggested that residue 50 is an important determinant of differential DNA-binding specificity among homeodomains. Although structures of several homeodomain-DNA complexes have been determined, this is the first structure of a native K50 class homeodomain. The only K50 homeodomain structure determined previously is an X-ray crystal structure of an altered specificity mutant, Engrailed Q50K (EnQ50K). Analysis of the NMR structure of the PITX2 homeodomain indicates that the lysine at position 50 makes contacts with two guanines on the antisense strand of the DNA, adjacent to the TAAT core DNA sequence, consistent with the structure of EnQ50K. Our evidence suggests that this side chain may make fluctuating interactions with the DNA, which is complementary to the crystal data for EnQ50K. There are differences in the tertiary structure between the native K50 structure and that of EnQ50K, which may explain differences in affinity and specificity between these proteins. Mutations in the human PITX2 gene are responsible for Rieger syndrome, an autosomal dominant disorder. Analysis of the residues mutated in Rieger syndrome indicates that many of these residues are involved in DNA binding, while others are involved in formation of the hydrophobic core of the protein. Overall, the role of K50 in homeodomain recognition is further clarified, and the results indicate that native K50 homeodomains may exhibit differences from altered specificity mutants.

Specific DNA recognition by theAntp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.

Solid-phase synthesis of selectively labeled DNA: applications for multidimensional nuclear magnetic resonance spectroscopy

The solid-phase chemical synthesis method has a strong advantage over the enzymatic method for preparing selectively labeled DNA oligomers. Atom-specific and fully labeled 2'-deoxynucleosides are economically prepared with routinely available isotope precursors using this synthetic route. Special DNA oligomers prepared by advanced labeling techniques are needed for advanced NMR applications, and chemical synthesis is the method of choice to respond to such demands. As a summary of this chapter, two tables are given. Table I lists the labeled nucleosides reported to be available by chemical syntheses. Table II lists the NMR studies using labeled DNA oligomers that were prepared by chemical syntheses.

An efficient NMR experiment for analyzing sugar-puckering in unlabeled DNA: application to the 26-kDa dead ringer-DNA complex

We present a new NMR experiment for estimating the type and degree of sugar-puckering in high-molecular-weight unlabeled DNA molecules. The experiment consists of a NOESY sequence preceded by a constant-time scalar coupling period. Two subexperiments are compared, each differing in the amount of time the (3)J(H3'H2') and (3)J(H3'H2") couplings are active on the H3' magnetization. The resultant data are easy to analyze, since a comparison of the signal intensities of any resolved NOE cross peak originating from H3' atoms of the duplex can be used to estimate the sum of the (3)J(H3'H2') and (3)J(H3'H2") couplings and thus the puckering type of the deoxyribose ring. Isotope filters to eliminate signals of the (13)C-labeled component in the F1-dimension are implemented, facilitating analyses of high-molecular-weight protein-DNA complexes containing (13)C-labeled protein and unlabeled DNA. The utility of the experiment is demonstrated on the 26-kDa Dead Ringer protein-DNA complex and reveals that the DNA uniformly adopts the S-type configuration when bound to protein.