The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Protein--DNA contacts in the structure of a homeodomain--DNA complex determined by nuclear magnetic resonance spectroscopy in solution.

The 1:1 complex of the mutant Antp(C39----S) homeodomain with a 14 bp DNA fragment corresponding to the BS2 binding site was studied by nuclear magnetic resonance (NMR) spectroscopy in aqueous solution. The complex has a molecular weight of 17,800 and its lifetime is long compared with the NMR chemical shift time scale. Investigations of the three-dimensional structure were based on the use of the fully 15N-labelled protein, two-dimensional homonuclear proton NOESY with 15N(omega 2) half-filter, and heteronuclear three-dimensional NMR experiments. Based on nearly complete sequence-specific resonance assignments, both the protein and the DNA were found to have similar conformations in the free form and in the complex. A sufficient number of intermolecular 1H-1H Overhauser effects (NOE) could be identified to enable a unique docking of the protein on the DNA, which was achieved with the use of an ellipsoid algorithm. In the complex there are intermolecular NOEs between the elongated second helix in the helix-turn-helix motif of the homeodomain and the major groove of the DNA. Additional NOE contacts with the DNA involve the polypeptide loop immediately preceding the helix-turn-helix segment, and Arg5. This latter contact is of special interest, both because Arg5 reaches into the minor groove and because in the free Antp(C39----S) homeodomain no defined spatial structure could be found for the apparently flexible N-terminal segment comprising residues 0-6.     

Flexibility and plasticity of human centrin 2 binding to the xeroderma pigmentosum group C protein (XPC) from nuclear excision repair

Human centrin 2 is a component of the nucleotide excision repair system, as a subunit of the heterotrimer including xeroderma pigmentosum group C protein (XPC) and hHR23B. The C-terminal domain of centrin (C-HsCen2) binds strongly a peptide from the XPC protein (P1-XPC: N(847)-R(863)). Here, we characterize the solution Ca(2+)-dependent structural and molecular features of the C-HsCen2 in complex with P1-XPC, mainly using NMR spectroscopy and molecular modeling. The N-terminal half of the peptide, organized as an alpha helix is anchored into a deep hydrophobic cavity of the protein, because of three bulky hydrophobic residues in position 1-4-8 and electrostatic contacts with the centrin helix E. Investigation of the whole centrin interactions shows that the N-terminal domain of the protein is not involved in the complex formation and is structurally independent from the peptide-bound C-terminal domain. The complex may exist in three different binding conformations corresponding to zero, one, and two Ca(2+)-bound states, which may exchange with various rates and have distinct structural stability. The various features of the intermolecular interaction presented here constitute a centrin-specific mode for the target binding.     

Intermolecular and intramolecular readout mechanisms in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.     

Geometric analysis and comparison of protein-DNA interfaces: why is there no simple code for recognition?

Structural studies of protein-DNA complexes have shown that there are many distinct families of DNA-binding proteins, and have shown that there is no simple "code" describing side-chain/base interactions. However, systematic analysis and comparison of protein-DNA complexes has been complicated by the diversity of observed contacts, the sheer number of complexes currently available and the absence of any consistent method of comparison that retains detailed structural information about the protein-DNA interface. To address these problems, we have developed geometric methods for characterizing the local structural environment in which particular side-chain/base interactions are observed. In particular, we develop methods for analyzing and comparing spatial relationships at the protein-DNA interface. Our method involves attaching local coordinate systems to the DNA bases and to the C(alpha) atoms of the peptide backbone (these are relatively rigid structural units). We use these tools to consider how the position and orientation of the polypeptide backbone (with respect to the DNA) helps to determine what contacts are possible at any given position in a protein-DNA complex. Here, we focus on base contacts that are made in the major groove, and we use spatial relationships in analyzing: (i) the observed patterns of side-chain/base interactions; (ii) observed helix docking orientations; (iii) family/subfamily relationships among DNA-binding proteins; and (iv) broader questions about evolution, altered specificity mutants and the limits for the design of new DNA-binding proteins. Our analysis, which highlights differences in spatial relationships in different complexes and at different positions in a complex, helps explain why there is no simple, general code for protein-DNA recognition.     

Dynamics and binding mode of Hoechst 33258 to d(GTGGAATTCCAC)2 in the 1:1 solution complex as determined by two-dimensional 1H NMR.

We have investigated the interaction of the bisbenzimidazole derivative Hoechst 33258 with the self-complementary dodecadeoxynucleotide duplex d(GTGGAATTCCAC)2 using one-dimensional (1D) and two-dimensional (2D) proton nuclear magnetic resonance (1H NMR) spectroscopy. To monitor the extent of complex formation, we used the imino proton region of the 1D 1H NMR spectra acquired in H2O solution. These spectra show that the DNA duplex loses its inherent C2v symmetry upon addition of the drug, indicating that the two molecules form a kinetically stable complex on the NMR time scale (the lifetime of the complex has been measured to be around 450 ms). We obtained sequence-specific assignments for all protons of the ligand and most protons of each separate strand of the oligonucleotide duplex using a variety of homonuclear 2D 1H NMR experiments. The aromatic protons of the DNA strands, which are symmetrically related in the free duplex, exhibit exchange cross peaks in the complex. This indicates that the drug binds in two equivalent sites on the 12-mer, with an exchange rate constant of 2.2 +/- 0.2 s-1. Twenty-five intermolecular NOEs were identified, all involving adenine 2 and sugar 1' protons of the DNA and protons in all four residues of the ligand, indicating that Hoechst 33258 is located in the minor groove at the AATT site. Only protons along the same edge of the two benzimidazole moieties of the drug show NOEs to DNA protons at the bottom of the minor groove. Using molecular mechanics, we have generated a unique model of the complex using distance constraints derived from the intermolecular NOEs. We present, however, evidence that the piperazine group may adopt at least two locally different conformations when the drug is bound to this dodecanucleotide.     

Experimental constraints on quaternary structure in Alzheimer's beta-amyloid fibrils

We describe solid-state nuclear magnetic resonance (NMR) measurements on fibrils formed by the 40-residue beta-amyloid peptide associated with Alzheimer's disease (Abeta(1-40)) that place constraints on the identity and symmetry of contacts between in-register, parallel beta-sheets in the fibrils. We refer to these contacts as internal and external quaternary contacts, depending on whether they are within a single molecular layer or between molecular layers. The data include (1) two-dimensional 13C-13C NMR spectra that indicate internal quaternary contacts between side chains of L17 and F19 and side chains of I32, L34, and V36, as well as external quaternary contacts between side chains of I31 and G37; (2) two-dimensional 15N-13C NMR spectra that indicate external quaternary contacts between the side chain of M35 and the peptide backbone at G33; (3) measurements of magnetic dipole-dipole couplings between the side chain carboxylate group of D23 and the side chain amine group of K28 that indicate salt bridge interactions. Isotopic dilution experiments allow us to make distinctions between intramolecular and intermolecular contacts. On the basis of these data and previously determined structural constraints from solid-state NMR and electron microscopy, we construct full molecular models using restrained molecular dynamics simulations and restrained energy minimization. These models apply to Abeta(1-40) fibrils grown with gentle agitation. We also present evidence for different internal quaternary contacts in Abeta(1-40) fibrils grown without agitation, which are morphologically distinct.     

The flexible loop L1 of the H3K4 demethylase JARID1B ARID domain has a crucial role in DNA-binding activity

JARID1B, a member of the JmjC demethylase family, has a crucial role in H3K4me3 demethylation. The ARID domain is a potential DNA-binding domain of JARID1B. Previous studies indicate that a GC-rich DNA motif is the specific target of the ARID domain. However, the details of the interaction between the ARID domain and duplex DNA require further study. Here, we utilized NMR spectroscopy to assign the backbone amino acids and mapped the DNA-binding sites of the human JARID1B ARID domain. Perturbations to ¹H-¹⁵N correlation spectra revealed that the flexible loop L1 of ARID was the main DNA-binding interface. EMSA and intrinsic fluorescence experiments demonstrated that mutations on loop L1 strongly reduced the DNA-binding activity of JARID1B ARID. Furthermore, transfection of mutant forms resulted in a distinct loss of intrinsic H3K4 demethylase activity, implying that the flexible loop L1 made a major contribution to sustaining the DNA-binding ability of JARID1B ARID domain.     

Structural basis for the binding affinity of a homologous series of synthetic phenoxazone drugs with DNA: NMR and molecular mechanics analysis

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5'-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear (1)H-(1)H 2D-NOESY and heteronuclear (1)H-(31)P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH(2) groups in the dimethylaminoalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH(3))(2) groups of the side chains.     

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 NH₃⁺ 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.     

An effective approach for generating a three-Cys<Subscript>2</Subscript>His<Subscript>2</Subscript> zinc-finger-DNA complex model by docking

Background  

Determination of protein-DNA complex structures with both NMR and X-ray crystallography remains challenging in many cases. High Ambiguity-Driven DOCKing (HADDOCK) is an information-driven docking program that has been used to successfully model many protein-DNA complexes. However, a protein-DNA complex model whereby the protein wraps around DNA has not been reported. Defining the ambiguous interaction restraints for the classical three-Cys₂His₂ zinc-finger proteins that wrap around DNA is critical because of the complicated binding geometry. In this study, we generated a Zif268-DNA complex model using three different sets of ambiguous interaction restraints (AIRs) to study the effect of the geometric distribution on the docking and used this approach to generate a newly reported Sp1-DNA complex model.     

Specificity and Affinity of Lac Repressor for the Auxiliary Operators O2 and O3 Are Explained by the Structures of Their Protein-DNA Complexes

The structures of a dimeric mutant of the Lac repressor DNA-binding domain complexed with the auxiliary operators O2 and O3 have been determined using NMR spectroscopy and compared to the structures of the previously determined Lac-O1 and Lac-nonoperator complexes. Structural analysis of the Lac-O1 and Lac-O2 complexes shows highly similar structures with very similar numbers of specific and nonspecific contacts, in agreement with similar affinities for these two operators. The left monomer of the Lac repressor in the Lac-O3 complex retains most of these specific contacts. However, in the right half-site of the O3 operator, there is a significant loss of protein-DNA contacts, explaining the low affinity of the Lac repressor for the O3 operator. The binding mode in the right half-site resembles that of the nonspecific complex. In contrast to the Lac-nonoperator DNA complex where no hinge helices are formed, the stability of the hinge helices in the weak Lac-O3 complex is the same as in the Lac-O1 and Lac-O2 complexes, as judged from the results of hydrogen/deuterium experiments.