Stepwise induced fit in the pico- to nanosecond time scale governs the complexation of the even-skipped transcriptional repressor homeodomain to DNA

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Dynamic Local Polymorphisms in the Gbx1 Homeodomain Induced by DNA Binding

The Gastrulation Brain Homeobox 1 (Gbx1) gene encodes the Gbx1 homeodomain that targets TAATTA motifs in double-stranded DNA (dsDNA). Residues Glu17 and Arg52 in Gbx1 form a salt bridge, which is preserved in crystal structures and molecular dynamics simulations of homologous homeodomain-DNA complexes. In contrast, our nuclear magnetic resonance (NMR) studies show that DNA binding to Gbx1 induces dynamic local polymorphisms, which include breaking of the Glu17-Arg52 salt bridge. To study this interaction, we produced a variant with Glu17Arg and Arg52Glu mutations, which exhibited the same fold as the wild-type protein, but a 2-fold reduction in affinity for dsDNA. Analysis of the NMR structures of the Gbx1 homeodomain in the free form, the Gbx1[E17R,R52E] variant, and a Gbx1 homeodomain-DNA complex showed that stabilizing interactions of the Arg52 side chain with the DNA backbone are facilitated by transient breakage of the Glu17-Arg52 salt bridge in the DNA-bound Gbx1.

Binding polymorphism in the DNA bound state of the Pdx1 homeodomain

The subtle effects of DNA-protein recognition are illustrated in the homeodomain fold. This is one of several small DNA binding motifs that, in spite of limited DNA binding specificity, adopts crucial, specific roles when incorporated in a transcription factor. The homeodomain is composed of a 3-helix domain and a mobile N-terminal arm. Helix 3 (the recognition helix) interacts with the DNA bases through the major groove, while the N-terminal arm becomes ordered upon binding a specific sequence through the minor groove. Although many structural studies have characterized the DNA binding properties of homeodomains, the factors behind the binding specificity are still difficult to elucidate. A crystal structure of the Pdx1 homeodomain bound to DNA (PDB 2H1K) obtained previously in our lab shows two complexes with differences in the conformation of the N-terminal arm, major groove contacts, and backbone contacts, raising new questions about the DNA recognition process by homeodomains. Here, we carry out fully atomistic Molecular Dynamics simulations both in crystal and aqueous environments in order to elucidate the nature of the difference in binding contacts. The crystal simulations reproduce the X-ray experimental structures well. In the absence of crystal packing constraints, the differences between the two complexes increase during the solution simulations. Thus, the conformational differences are not an artifact of crystal packing. In solution, the homeodomain with a disordered N-terminal arm repositions to a partially specific orientation. Both the crystal and aqueous simulations support the existence of different stable binding conformers identified in the original crystallographic data with different degrees of specificity. We propose that protein-protein and protein-DNA interactions favor a subset of the possible conformations. This flexibility in DNA binding may facilitate multiple functions for the same transcription factor.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

M.TaqI facilitates the base flipping via an unusual DNA backbone conformation

MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.

The N6-Methyl Group of Adenine Further Increases the BI Stability of DNA Compared to C5-Methyl Groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.

Water-mediated contacts in thetrp-repressor operator complex recognition process

Water-mediated contacts are known as an important recognition tool in trp-repressor operator systems. One of these contacts involves two conserved base pairs (G(6).C(-6) and A(5). T(-5)) and three amino acids (Lys 72, Ile 79, and Ala 80). To investigate the nature of these contacts, we analyzed the X-ray structure (PDB code: 1TRO) of the trp-repressor operator complex by means of molecular dynamics simulations. This X-ray structure contains two dimers that exhibit structural differences. From these two different starting structures, two 10 ns molecular dynamics simulations have been performed. Both of our simulations show an increase of water molecules in the major groove at one side of the dimer, while the other side remains unchanged compared to the X-ray structure. Though the maximum residence time of the concerned water molecules decreases with an increase of solvent at the interface, these water molecules continue to play an important role in mediating DNA-protein contacts. This is shown by new stable amino acids-DNA distances and a long water residence time compared to free DNA simulation. To maintain stability of the new contacts, the preferential water binding site on O6(G6) is extended. This extension agrees with mutation experiment data on A5 and G6, which shows different relative affinity due to mutation on these bases [A. Joachimiak, T. E. Haran, P. B. Sigler, EMBO Journal 1994, Vol. 13, No. (2) pp. 367-372]. Due to the rearrangements in the system, the phosphate of the base G6 is able to interconvert to the B(II) substate, which is not observed on the other half side of the complex. The decrease of the number of hydrogen bonds between protein and DNA backbone could be the initial step of the dissociation process of the complex, or in other words an intermediate complex conformation of the association process. Thus, we surmise that these features show the importance of water-mediated contacts in the trp-repressor operator recognition process.