Molecular dynamics simulation reveals sequence-intrinsic and protein-induced geometrical features of the OL1 DNA operator

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

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.

Molecular dynamics simulations of papilloma virus E2 DNA sequences: dynamical models for oligonucleotide structures in solution

The specificity of papilloma virus E2 protein-DNA binding depends critically upon the sequence of a region of the DNA not in direct contact with the protein, and represents one of the simplest known examples of indirect readout. A detailed characterization of this system in solution is important to the further investigation hypothesis of a structural code for DNA recognition by regulatory proteins. In the crystalline state, the E2 DNA oligonucleotide sequence, d(ACCGAATTCGGT), exhibits three different structural forms. We report herein studies of the structure of E2 DNA in solution based on a series of molecular dynamics (MD) simulations including counterions and water, utilizing both the canonical and various crystallographic structures as initial points of departure. All MDs converged on a single dynamical structure of d(ACCGAATTCGGT) in solution. The predicted structure is in close accord with two of the three crystal structures, and indicates that a significant kink in the double helix at the central ApT step in the other crystal molecule may be a packing effect. The dynamical fine structure was analyzed on the basis of helicoidal parameters. The calculated curvature in the sequence was found to originate primarily from YPR steps in the regions flanking the central AATT tract. In order to study the role of structural adaptation of the DNA in the binding process, a subsequent simulation on the 16-mer cognate sequence d(CAACCGAATTCGGTTG) was initiated from the crystallographic coordinates of the bound DNA in the crystal structure of the protein DNA complex. MD simulations starting with the protein-bound form relaxed rapidly back to the dynamical structure predicted from the previous simulations on the uncomplexed DNA. The MD results show that the bound form E2 DNA is a dynamically unstable structure in the absence of protein, and arises as a consequence of both structural changes intrinsic to the sequence and induced by the interaction with protein.

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.