Modelling Protein-DNA Interactions

Electrostatic hot spot on DNA-binding domains mediates phosphate desolvation and the pre-organization of specificity determinant side chains

A major obstacle towards elucidating the molecular basis of transcriptional regulation is the lack of a detailed understanding of the interplay between non-specific and specific protein–DNA interactions. Based on molecular dynamics simulations of C₂H₂ zinc fingers (ZFs) and engrailed homeodomain transcription factors (TFs), we show that each of the studied DNA-binding domains has a set of highly constrained side chains in preset configurations ready to form hydrogen bonds with the DNA backbone. Interestingly, those domains that bury their recognition helix into the major groove are found to have an electrostatic hot spot for Cl⁻ ions located on the same binding cavity as the most buried DNA phosphate. The spot is characterized by three protein hydrogen bond donors, often including two basic side chains. If bound, Cl⁻ ions, likely mimicking phosphates, steer side chains that end up forming specific contacts with bases into bound-like conformations. These findings are consistent with a multi-step DNA-binding mechanism in which a pre-organized set of TF side chains assist in the desolvation of phosphates into well defined sites, prompting the re-organization of specificity determining side chains into conformations suitable for the recognition of their cognate sequence.

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

Operator recognition by the phage 434 cI repressor: MD simulations of free and bound 50-bp DNA reveal important differences between the OR1 and OR2 sites

Using molecular dynamics simulations in explicit solvent, we investigated the behavior of a 50-bp DNA sequence containing the 434 bacteriophage operators OR1 and OR2 separated by an 8-bp spacer. Two simulations of 1 ns each were carried out, with DNA alone and with DNA complexed to dimers of the R1-69 DNA binding domain of the phage 434 cI repressor protein at the OR1 and OR2 sites. Strong correlations among average structural parameters are observed between our simulations and available experimental data for the bound OR1/OR2 subsites. In the free state, some differences appear between the three relevant fragments (OR1, the spacer, and OR2). Unbound OR1 exhibits a large, shallow major groove into which the base atoms protrude and is also bent toward the major groove. This structure is maintained because structural fluctuations are weak. Unbound OR2 resembles canonical B-DNA although the structural parameters show greater fluctuations, essentially due to a malleable step (the innermost CpA/TpG), absent in OR1. Complexation with the proteins slightly alters the base positions but strongly modifies the sugar and backbone motions. The most crucial repressor effects are changes in the flexibility of the OR1/OR2 sites. Structural fluctuations are enhanced for OR1, conferring a favorable energetic contribution to the OR1 binding, whereas they are reduced for OR2. Therefore, both structural and dynamic properties of DNA suggest OR1 is the most attractive site for the repressor, which may explain the different binding association constants observed for the OR1 and OR2 sites. Finally, we also investigated the impact of the protein on the DNA backbone dynamics and find that direct or indirect interactions facilitate the DNA structural variations required for achieving complementarity with the protein.