Search

Interfacial water molecules play an important role in many aspects of protein-DNA specificity and recognition. Yet they have been mostly neglected in the computational modeling of these complexes. We present here a solvated docking protocol that allows explicit inclusion of water molecules in the docking of protein-DNA complexes and demonstrate its feasibility on a benchmark of 30 high-resolution protein-DNA complexes containing crystallographically-determined water molecules at their interfaces. Our protocol is capable of reproducing the solvation pattern at the interface and recovers hydrogen-bonded water-mediated contacts in many of the benchmark cases. Solvated docking leads to an overall improvement in the quality of the generated protein-DNA models for cases with limited conformational change of the partners upon complex formation. The applicability of this approach is demonstrated on real cases by docking a representative set of 6 complexes using unbound protein coordinates, model-built DNA and knowledge-based restraints. As HADDOCK supports the inclusion of a variety of NMR restraints, solvated docking is also applicable for NMR-based structure calculations of protein-DNA complexes.

Inaccuracies in computational molecular modeling methods are often counterweighed by brute-force generation of a plethora of putative solutions. These are then typically sieved via structural clustering based on similarity measures such as the root mean square deviation (RMSD) of atomic positions. Albeit widely used, these measures suffer from several theoretical and technical limitations (e.g., choice of regions for fitting) that impair their application in multicomponent systems (N > 2), large-scale studies (e.g., interactomes), and other time-critical scenarios. We present here a simple similarity measure for structural clustering based on atomic contacts--the fraction of common contacts--and compare it with the most used similarity measure of the protein docking community--interface backbone RMSD. We show that this method produces very compact clusters in remarkably short time when applied to a collection of binary and multicomponent protein-protein and protein-DNA complexes. Furthermore, it allows easy clustering of similar conformations of multicomponent symmetrical assemblies in which chain permutations can occur. Simple contact-based metrics should be applicable to other structural biology clustering problems, in particular for time-critical or large-scale endeavors., (Copyright © 2012 Wiley Periodicals, Inc.)

The intrinsic flexibility of DNA and the difficulty of identifying its interaction surface have long been challenges that prevented the development of efficient protein-DNA docking methods. We have demonstrated the ability our flexible data-driven docking method HADDOCK to deal with these before, by using custom-built DNA structural models. Here we put our method to the test on a set of 47 complexes from the protein-DNA docking benchmark. We show that HADDOCK is able to predict many of the specific DNA conformational changes required to assemble the interface(s). Our DNA analysis and modelling procedure captures the bend and twist motions occurring upon complex formation and uses these to generate custom-built DNA structural models, more closely resembling the bound form, for use in a second docking round. We achieve throughout the benchmark an overall success rate of 94% of one-star solutions or higher (interface root mean square deviation ≤4 A and fraction of native contacts >10%) according to CAPRI criteria. Our improved protocol successfully predicts even the challenging protein-DNA complexes in the benchmark. Finally, our method is the first to readily dock multiple molecules (N > 2) simultaneously, pushing the limits of what is currently achievable in the field of protein-DNA docking.

There is a growing interest in structural studies of DNA by both experimental and computational approaches. Often, 3D-structural models of DNA are required, for instance, to serve as templates for homology modeling, as starting structures for macro-molecular docking or as scaffold for NMR structure calculations. The conformational adaptability of DNA when binding to a protein is often an important factor and at the same time a limitation in such studies. As a response to the demand for 3D-structural models reflecting the intrinsic plasticity of DNA we present the 3D-DART server (3DNA-Driven DNA Analysis and Rebuilding Tool). The server provides an easy interface to a powerful collection of tools for the generation of DNA-structural models in custom conformations. The computational engine beyond the server makes use of the 3DNA software suite together with a collection of home-written python scripts. The server is freely available at http://haddock.chem.uu.nl/dna without any login requirement.

We present a protein-DNA docking benchmark containing 47 unbound-unbound test cases of which 13 are classified as easy, 22 as intermediate and 12 as difficult cases. The latter shows considerable structural rearrangement upon complex formation. DNA-specific modifications such as flipped out bases and base modifications are included. The benchmark covers all major groups of DNA-binding proteins according to the classification of Luscombe et al., except for the zipper-type group. The variety in test cases make this non-redundant benchmark a useful tool for comparison and development of protein-DNA docking methods. The benchmark is freely available as download from the internet.

Intrinsic flexibility of DNA has hampered the development of efficient protein-DNA docking methods. In this study we extend HADDOCK (High Ambiguity Driven DOCKing) [C. Dominguez, R. Boelens and A. M. J. J. Bonvin (2003) J. Am. Chem. Soc. 125, 1731-1737] to explicitly deal with DNA flexibility. HADDOCK uses non-structural experimental data to drive the docking during a rigid-body energy minimization, and semi-flexible and water refinement stages. The latter allow for flexibility of all DNA nucleotides and the residues of the protein at the predicted interface. We evaluated our approach on the monomeric repressor-DNA complexes formed by bacteriophage 434 Cro, the Escherichia coli Lac headpiece and bacteriophage P22 Arc. Starting from unbound proteins and canonical B-DNA we correctly predict the correct spatial disposition of the complexes and the specific conformation of the DNA in the published complexes. This information is subsequently used to generate a library of pre-bent and twisted DNA structures that served as input for a second docking round. The resulting top ranking solutions exhibit high similarity to the published complexes in terms of root mean square deviations, intermolecular contacts and DNA conformation. Our two-stage docking method is thus able to successfully predict protein-DNA complexes from unbound constituents using non-structural experimental data to drive the docking.

Recognition of the lac operator by the lac repressor involves specific interactions between residues in the repressor's recognition helix and bases in the DNA major groove. Tyr17 and Gln18, at positions 1 and 2 in the lac repressor recognition helix, can be exchanged for other amino acids to generate mutant repressors that display altered specificity. We have solved the solution structure of a protein-DNA complex of an altered-specificity mutant lac headpiece in which Tyr17 and Gln18 were exchanged for valine and alanine, respectively, as found in the recognition helix of the gal repressor. As previously described by Lehming et al. (EMBO J. 1987, 6, 3145-3153), this altered-specificity mutant of the lac repressor recognizes a variant lac operator that is similar to the gal operator Oe. The mutant lac headpiece showed the predicted specificity and is also able to mimic the gal repressor by recognizing and bending the natural gal operator Oe. These structural data show that, while most of the anchoring points that help the lac headpiece to assemble on the lac operator were preserved, a different network of protein-DNA interactions connecting Ala17 and Val18 to bases in the DNA major groove drives the specificity towards the altered operator.

The lac repressor-operator system is a model system for understanding protein-DNA interactions and allosteric mechanisms in gene regulation. Despite the wealth of biochemical data provided by extensive mutations of both repressor and operator, the specific recognition mechanism of the natural lac operators by lac repressor has remained elusive. Here we present the first high-resolution structure of a dimer of the DNA-binding domain of lac repressor bound to its natural operator 01. The global positioning of the dimer on the operator is dramatically asymmetric, which results in a different pattern of specific contacts between the two sites. Specific recognition is accomplished by a combination of elongation and twist by 48 degrees of the right lac subunit relative to the left one, significant rearrangement of many side chains as well as sequence-dependent deformability of the DNA. The set of recognition mechanisms involved in the lac repressor-operator system is unique among other protein-DNA complexes and presents a nice example of the adaptability that both proteins and DNA exhibit in the context of their mutual interaction.

The retinoid X receptor (RXR) is a prominent member of the nuclear receptor family of ligand-inducible transcription factors. Many proteins of this family exert their function as heterodimers with RXR as a common upstream partner. Studies of the DNA-binding domains of several nuclear receptors reveal differences in structure and dynamics, both between the different proteins and between the free- and DNA-bound receptor DBDs. We investigated the differences in dynamics between RXR free in solution and in complex with a 14 base-pair oligonucleotide, using (1)H and (15)N relaxation studies. Nano- to picosecond dynamics were probed on (15)N, employing Lipari-Szabo analysis with an axially symmetric tumbling model to estimate the exchange contributions to the transverse relaxation rates. Furthermore, milli- to microsecond dynamics were estimated qualitatively for (1)H and (15)N, using CPMG-HSQC and CPMG-T(2) measurements with differential pulse spacing. RXR shows hardly any nano- to picosecond time-scale internal motion. Upon DNA binding, the order parameters show a tiny increase. Dynamics in the milli- to microsecond time scale is more prevalent. It is localized in the first and second zinc fingers of the free RXR. Upon DNA-binding, exchange associated with specific/aspecific DNA-binding of RXR is observed throughout the sequence, whereas conformational flexibility of the D-box and the second zinc finger of RXR is greatly reduced. Since this DNA-binding induced folding transition occurs remote from the DNA in a region which is involved in protein-protein interactions, it may very well be related to the cooperativity of dimeric DNA binding.

The localisation and dynamics of sodium counterions around the DNA duplex d(AGCGTACTAGTACGCT)2 corresponding to the trp operator fragment used in the crystal structure of the half site complex has been studied by a 1.4 ns molecular dynamics simulation in explicit solvent. A continuous and well-defined counterion density is shown to be present around the minor groove, while density patches are found in the major groove in regions where DNA bending is observed. A residence time analysis reveals the dynamic nature of these distributions. The resulting picture agrees with previous theoretical and experimental studies of A-tract DNA sequences, and is consistent with the polyelectrolyte condensation model.

The structure and hydration of the DNA duplex d-(AGCGTACTAGTACGCT)2 corresponding to the trp operator fragment used in the crystal structure of the half site complex (PDB entry 1TRR) was studied by a 1.4 ns molecular dynamics simulation in water. The simulation, starting from a B-DNA conformation, used a non-bonded cutoff of 1.4 nm with a reaction field correction and resulted in a stable trajectory. The average DNA conformation obtained was closer to the ones found in the crystal structures of the complexes (PDB entries 1TRO and 1TRR) than to the crystal structure of unbound trp operator (Nucleic Acid Database entry BDJ061). The DNA hydration was characterized in terms of hydrogen bond percentages and corresponding residence times. The residence times of water molecules within 0.35 nm of the DNA non-exchangeable protons were calculated for comparison with NMR measurements of intermolecular water-DNA NOEs and nuclear magnetic relaxation dispersion measurements. No significant difference was found between major and minor groove hydration. The DNA donors and acceptors were hydrogen bonded to water molecules for 77(+/-19)% of the time on average. The average residence time of the hydrogen bonded water molecules was 11(+/-11) ps with a maximum of 223 ps. When all water molecules within NOE distance (0.35 nm) of non-exchangeable protons were considered, the average residence times increased to an average of 100(+/-4) ps and a maximum of 608 ps. These results agree with the experimental NMR results of Sunnerhagen et al. which did not show any evidence for water molecules bound with more than 1 ns residence time on the DNA surface. The exchange of hydration water from the DNA occurred in the major groove primarily through direct exchange with the bulk solvent, while access to and from the minor groove frequently proceeded via pathways involving ribose O3' and O4' and phosphate O2P oxygen atoms. The most common water diffusion pathways in the minor groove were perpendicular to the groove direction. In general, water molecules visited only a limited number of sites in the DNA grooves before exiting. The hydrogen bonding sites, where hydrogen bonds could be formed with donor and acceptor groups of the DNA, were filled with water molecules with an average B-factor value of 0.58 mn2. No special values were observed at any of the sites, where water molecules were observed both in the trp repressor/operator co-crystals and in the crystal structure of unbound DNA., (Copyright 1998 Academic Press.)

The present NMR study investigates the residence times of the hydration water molecules associated with uncomplexed trp operator DNA in solution by measuring intermolecular nuclear Overhauser effects (NOE) between water and DNA protons, and the nuclear magnetic relaxation dispersion (NMRD) of the water 2H and 17O resonances. Both methods indicate that the hydration water molecules exchange with bulk water on the sub-nanosecond time scale at 4 degreesC. No evidence was obtained for water molecules bound with longer residence times. In particular, the water molecules at the sites of interfacial hydration in the trp repressor/operator complex do not seem kinetically stabilized in the uncomplexed DNA. Analysis of the crystal structures of two different trp repressor/operator complexes shows very similar structural environments for the water molecules mediating specific contacts between the protein and the DNA, whereas much larger variations are observed for the location of corresponding water molecules detected in the crystal structure of an uncomplexed trp operator DNA duplex. Therefore, it appears unlikely that the hydration characteristics of the uncomplexed DNA target would be a major determinant of trp repressor/operator recognition., (Copyright 1998 Academic Press.)

The three-dimensional structure of the DNA-binding domain of the human retinoic acid receptor-beta (hRAR-beta) has been determined by nuclear magnetic resonance spectroscopy in conjunction with distance geometry, restrained molecular dynamics and iterative relaxation matrix calculations. A total of 1244 distance restraints were obtained from NOE intensities, of which 448 were intra-residue and 796 inter-residue restraints. In addition 23 chi and 30 phi dihedral angle restraints were obtained from J-coupling data. The two 'zinc-finger' regions of the 80-amino acid residue protein are followed by two alpha-helices that cross each other perpendicularly. There is a short stretch of b-sheet near the N-terminus. The alpha-helical core of the protein is well determined with a backbone root-mean-square deviation (r.m.s.d.) with respect to the average of 0.18 A and 0.37 A when the side chains of residues 31, 32, 36, 61, 62, 65 and 69 are included. The r.m.s.d. for the backbone of residues 5-80 is 0.76 A. For the first finger (residues 8-28), the r.m.s.d. of the backbone is 0.79 A. For the second finger (residues 44-62) the r.m.s.d. is 0.64 A. The overall structure is similar to that of the corresponding domain of the glucocorticoid receptor, although the C-terminal part of the protein is different. The second alpha-helix is two residues shorter and is followed by a well-defined region of extended backbone structure.