Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

Clinical and Biochemical Function of Polymorphic NR0B1 GGAA-Microsatellites in Ewing Sarcoma: A Report from the Children's Oncology Group

Background

The genetics involved in Ewing sarcoma susceptibility and prognosis are poorly understood. EWS/FLI and related EWS/ETS chimeras upregulate numerous gene targets via promoter-based GGAA-microsatellite response elements. These microsatellites are highly polymorphic in humans, and preliminary evidence suggests EWS/FLI-mediated gene expression is highly dependent on the number of GGAA motifs within the microsatellite.

Objectives

Here we sought to examine the polymorphic spectrum of a GGAA-microsatellite within the NR0B1 promoter (a critical EWS/FLI target) in primary Ewing sarcoma tumors, and characterize how this polymorphism influences gene expression and clinical outcomes.

Results

A complex, bimodal pattern of EWS/FLI-mediated gene expression was observed across a wide range of GGAA motifs, with maximal expression observed in constructs containing 20–26 GGAA motifs. Relative to white European and African controls, the NR0B1 GGAA-microsatellite in tumor cells demonstrated a strong bias for haplotypes containing 21–25 GGAA motifs suggesting a relationship between microsatellite function and disease susceptibility. This selection bias was not a product of microsatellite instability in tumor samples, nor was there a correlation between NR0B1 GGAA-microsatellite polymorphisms and survival outcomes.

Conclusions

These data suggest that GGAA-microsatellite polymorphisms observed in human populations modulate EWS/FLI-mediated gene expression and may influence disease susceptibility in Ewing sarcoma.     

Structure specific DNA recognition by the SLX1–SLX4 endonuclease complex

The SLX1–SLX4 structure-specific endonuclease complex is involved in processing diverse DNA damage intermediates, including resolution of Holliday junctions, collapse of stalled replication forks and removal of DNA flaps. The nuclease subunit SLX1 is inactive on its own, but become activated upon binding to SLX4 via its conserved C-terminal domain (CCD). Yet, how the SLX1–SLX4 complex recognizes specific DNA structure and chooses cleavage sites remains unknown. Here we show, through a combination of structural, biochemical and computational analyses, that the SAP domain of SLX4 is critical for efficient and accurate processing of 5′-flap DNA. It binds the minor groove of DNA about one turn away from the flap junction, and the 5′-flap is implicated in binding the core domain of SLX1. This binding mode accounts for specific recognition of 5′-flap DNA and specification of cleavage site by the SLX1–SLX4 complex.     

The crystal structure of the ternary complex of T.thermophilus seryl-tRNA synthetase with tRNA(Ser) and a seryl-adenylate analogue reveals a conformational switch in the active site.

The low temperature crystal structure of the ternary complex of Thermus thermophilus seryl-tRNA synthetase with tRNA(Ser) (GGA) and a non-hydrolysable seryl-adenylate analogue has been refined at 2.7 angstrom resolution. The analogue is found in both active sites of the synthetase dimer but there is only one tRNA bound across the two subunits. The motif 2 loop of the active site into which the single tRNA enters interacts within the major groove of the acceptor stem. In particular, a novel ring-ring interaction between Phe262 on the extremity of this loop and the edges of bases U68 and C69 explains the conservation of pyrimidine bases at these positions in serine isoaccepting tRNAs. This active site takes on a significantly different ordered conformation from that observed in the other subunit, which lacks tRNA. Upon tRNA binding, a number of active site residues previously found interacting with the ATP or adenylate now switch to participate in tRNA recognition. These results shed further light on the structural dynamics of the overall aminoacylation reaction in class II synthetases by revealing a mechanism which may promote an ordered passage through the activation and transfer steps.     

DNA nicking by HinP1I endonuclease: bending, base flipping and minor groove expansion

HinP1I recognizes and cleaves the palindromic tetranucleotide sequence G↓CGC in DNA. We report three structures of HinP1I–DNA complexes: in the presence of Ca²⁺ (pre-reactive complex), in the absence of metal ion (binary complex) and in the presence of Mg²⁺ (post-reactive complex). HinP1I forms a back-to-back dimer with two active sites and two DNA duplexes bound on the outer surfaces of the dimer facing away from each other. The 10 bp DNA duplexes undergo protein-induced distortions exhibiting features of A-, B- and Z-conformations: bending on one side (by intercalation of a phenylalanine side chain into the major groove), base flipping on the other side of the recognition site (by expanding the step rise distance of the local base pair to Z-form) and a local A-form conformation between the two central C:G base pairs of the recognition site (by binding of the N-terminal helix in the minor groove). In the pre- and post-reactive complexes, two metals (Ca²⁺ or Mg²⁺) are found in the active site. The enzyme appears to cleave DNA sequentially, hydrolyzing first one DNA strand, as seen in the post-reactive complex in the crystalline state, and then the other, as supported by the observation that, in solution, a nicked DNA intermediate accumulates before linearization.     

Incoming nucleotide binds to Klenow ternary complex leading to stable physical sequestration of preceding dNTP on DNA

Klenow–DNA complex is known to undergo a rate-limiting, protein conformational transition from an ‘open’ to ‘closed’ state, upon binding of the ‘correct’ dNTP at the active site. In the ‘closed’ state, Mg²⁺ mediates a rapid chemical step involving nucleophilic displacement of pyrophosphate by the 3′ hydroxyl of the primer terminus. The enzyme returns to the ‘open’ state upon the release of PPi and translocation permits the next round of reaction. To determine whether Klenow can translocate to the next site on the addition of the next dNTP, without the preceding chemical step, we studied the ternary complex (Klenow–DNA–dNTP) in the absence of Mg²⁺. While the ternary complex is proficient in chemical addition of dNTPs in Mg²⁺, as revealed by primer extensions, the same in Mg²⁺-deficient conditions lead to non-covalent (physical) sequestration of first two ‘correct’ dNTPs in the ternary complex. Moreover, the second dNTP traps the first one in the DNA-helix of the ternary complex. Such a dNTP–DNA complex is found to be stable even after the dissociation of Klenow. This reveals the novel state of the dNTP–DNA complex where the complementary base is stacked in a DNA-helix non-covalently, without the phosphodiester linkage. Further, shuttling of the DNA between the polymerase and the exonuclease site mediates the release of such a DNA complex. Interestingly, Klenow in such a Mg²⁺-deficient ternary complex exhibits a ‘closed’ conformation.     

Elucidation of the mode of interaction in the UP1-telomerase RNA-telomeric DNA ternary complex which serves to recruit telomerase to telomeric DNA and to enhance the telomerase activity

We found that UP1, a proteolytic product of heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1), both enhances and represses the telomerase activity. The formation of the UP1–telomerase RNA–telomeric DNA ternary complex was revealed by a gel retardation experiment. The interactions in the ternary and binary complexes were elucidated by NMR. UP1 has two nucleic acid-binding domains, BD1 and BD2. In the UP1–telomerase RNA binary complex, both BD1 and BD2 interact with telomerase RNA. Interestingly, when telomeric DNA was added to the binary complex, telomeric DNA bound to BD1 in place of telomerase RNA. Thus, BD1 basically binds to telomeric DNA, while BD2 mainly binds to telomerase RNA, which resulted in the formation of the ternary complex. Here, UP1 bridges telomerase and telomeric DNA. It is supposed that UP1/hnRNP A1 serves to recruit telomerase to telomeric DNA through the formation of the ternary complex. A model has been proposed for how hnRNP A1/UP1 contributes to enhancement of the telomerase activity through recruitment and unfolding of the quadruplex of telomeric DNA.     

Caught in the act: visualization of an intermediate in the DNA base-flipping pathway induced by HhaI methyltransferase

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket (‘base flipping’) is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI–DNA–AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.     

Investigation of a conserved stacking interaction in target site recognition by the U1A protein

Three highly conserved aromatic residues in RNA recognition motifs (RRM) participate in stacking interactions with RNA bases upon binding RNA. We have investigated the contribution of one of these aromatic residues, Phe56, to the complex formed between the N-terminal RRM of the spliceosomal protein U1A and stem–loop 2 of U1 snRNA. Previous work showed that the aromatic group is important for high affinity binding. Here we probe how mutation of Phe56 affects the kinetics of complex dissociation, the strength of the hydrogen bonds formed between U1A and the base that stacks with Phe56 (A6) and specific target site recognition. Substitution of Phe56 with Trp or Tyr increased the rate of dissociation of the complex, consistent with previously reported results. However, substitution of Phe56 with His decreased the rate of complex association, implying a change in the initial formation of the complex. Simultaneous modification of residue 56 and A6 revealed energetic coupling between the aromatic group and the functional groups of A6 that hydrogen bond to U1A. Finally, mutation of Phe56 to Leu reduced the ability of U1A to recognize stem–loop 2 correctly. Taken together, these experiments suggest that Phe56 contributes to binding affinity by stacking with A6 and participating in networks of energetically coupled interactions that enable this conserved aromatic amino acid to play a complex role in target site recognition.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Simulations of DNA topoisomerase 1B bound to supercoiled DNA reveal changes in the flexibility pattern of the enzyme and a secondary protein–DNA binding site

Human topoisomerase 1B has been simulated covalently bound to a negatively supercoiled DNA minicircle, and its behavior compared to the enzyme bound to a simple linear DNA duplex. The presence of the more realistic supercoiled substrate facilitates the formation of larger number of protein–DNA interactions when compared to a simple linear duplex fragment. The number of protein–DNA hydrogen bonds doubles in proximity to the active site, affecting all of the residues in the catalytic pentad. The clamp over the DNA, characterized by the salt bridge between Lys369 and Glu497, undergoes reduced fluctuations when bound to the supercoiled minicircle. The linker domain of the enzyme, which is implicated in the controlled relaxation of superhelical stress, also displays an increased number of contacts with the minicircle compared to linear DNA. Finally, the more complex topology of the supercoiled DNA minicircle gives rise to a secondary DNA binding site involving four residues located on subdomain III. The simulation trajectories reveal significant changes in the interactions between the enzyme and the DNA for the more complex DNA topology, which are consistent with the experimental observation that the protein has a preference for binding to supercoiled DNA.     

Sequence specific DNA binding of Ets-1 transcription factor: molecular dynamics study on the Ets domain--DNA complexes

Molecular dynamics (MD) simulations for Ets-1 ETS domain-DNA complexes were performed to investigate the mechanism of sequence-specific recognition of the GGAA DNA core by the ETS domain. Employing the crystal structure of the Ets-1 ETS domain-DNA complex as a starting structure we carried out MD simulations of: (i). the complex between Ets-1 ETS domain and a 14 base-pair DNA containing GGAA core sequence (ETS-GGAA); (ii). the complex between the ETS domain and a DNA having single base-pair mutation, GGAG sequence (ETS-GGAG); and (iii). the 14 base-pair DNA alone (GGAA). Comparative analyses of the MD structures of ETS-GGAA and ETS-GGAG reveal that the DNA bending angles and the ETS domain-DNA phosphate interactions are similar in these complexes. These results support that the GGAA core sequence is distinguished from the mutated GGAG sequence by a direct readout mechanism in the Ets-1 ETS domain-DNA complex. Further analyses of the direct contacts in the interface between the helix-3 region of Ets-1 and the major groove of the core DNA sequence clearly show that the highly conserved arginine residues, Arg391 and Arg394, play a critical role in binding to the GGAA core sequence. These arginine residues make bidentate contacts with the nucleobases of GG dinucleotides in GGAA core sequence. In ETS-GGAA, the hydroxyl group of Tyr395 is hydrogen bonded to N7 nitrogen of A(3) (the third adenosine in the GGAA core), while the hydroxyl group makes a contact with N4 nitrogen of C(4') (the complementary nucleotide of the fourth guanosine G(4) in the GGAG sequence) in the ETS-GGAG complex. We have found that this difference in behavior of Tyr395 results in the relatively large motion of helix-3 in the ETS-GGAG complex, causing the collapse of bidentate contacts between Arg391/Arg394 and the GG dinucleotides in the GGAG sequence.     

Role of the protein in the DNA sequence specificity of the cleavage site stabilized by the camptothecin topoisomerase IB inhibitor: a metadynamics study

Camptothecin (CPT) is a topoisomerase IB (TopIB) selective inhibitor whose derivatives are currently used in cancer therapy. TopIB cleaves DNA at any sequence, but in the presence of CPT the only stabilized protein–DNA covalent complex is the one having a thymine in position −1 with respect to the cleavage site. A metadynamics simulation of two TopIB–DNA–CPT ternary complexes differing for the presence of a thymine or a cytosine in position −1 indicates the occurrence of two different drug’s unbinding pathways. The free-energy difference between the bound state and the transition state is large when a thymine is present in position −1 and is strongly reduced in presence of a cytosine, in line with the different drug stabilization properties of the two systems. Such a difference is strictly related to the changes in the hydrogen bond network between the protein, the DNA and the drug in the two systems, indicating a direct role of the protein in determining the specificity of the cleavage site sequence stabilized by the CPT. Calculations carried out in presence of one compound of the indenoisoquinoline family (NSC314622) indicate a comparable energy difference between the bound and the transition state independently of the presence of a thymine or a cytosine in position −1, in line with the experimental results.     

Replicative helicases can translocate through abasic site-induced covalent topoisomerase IV-DNA complexes

Some topoisomerase inhibitors trap covalent topoisomerase–DNA complexes as topoisomerase–drug–DNA ternary complexes. Ternary complex formation results in inhibition of DNA replication and generation of permanent double-strand breaks. Recent demonstrations of the stimulation of covalent topoisomerase–DNA complex formation by DNA lesions suggest that DNA damage may act as an endogenous topoisomerase poison. We have investigated the effects of abasic (AP) sites on topoisomerase IV (Topo IV). AP sites can stimulate the formation of covalent Topo IV–DNA complexes when they are located either within the 4 base overhang generated by DNA scission or immediately 5′ to the point of scission (the –1 position). Thus, the AP site acts as a position-specific, endogenous topoisomerase poison. Both EDTA and salt can reverse covalent Topo IV–DNA complexes induced by AP sites located within the 4 base overhang. Interestingly, an AP site at the –1 position inhibits EDTA-mediated reversal of formation of the covalent Topo IV–DNA complex. Furthermore, we find that, unlike quinolone-induced covalent Topo IV–DNA complexes, AP site-induced covalent Topo IV–DNA complexes do not inhibit the helicase activities of the DnaB and T7 Gene 4 proteins. These results suggest that the AP site-induced poisoning of Topo IV does not arrest replication fork progression.     

Crystal structures of a ddATP-, ddTTP-, ddCTP, and ddGTP- trapped ternary complex of Klentaq1: insights into nucleotide incorporation and selectivity

The mechanism by which DNA polymerase I enzymes function has been the subject of extensive biochemical and structural studies. We previously determined the structure of a ternary complex of the large fragment of DNA polymerase I from Thermus aquaticus (Klentaq1) bound to a primer/template DNA and a dideoxycytidine 5'-triphosphate (ddCTP). In this report, we present the details of the 2.3-A resolution crystal structures of three additional ternary complexes of Klentaq1 bound to a primer/template DNA and a dideoxyguanosine 5'-triphosphate (ddGTP), a dideoxythymidine 5'-triphosphate (ddTTP), or a dideoxyadenosine 5'-triphosphate (ddATP). Comparison of the active site of the four ternary complexes reveals that the protein residues around the nascent base pair (that formed between the incoming dideoxynucleoside triphosphate [ddNTP] and the template base) form a snug binding pocket into which only a correct Watson-Crick base pair can fit. Except in the ternary complex bound to dideoxyguanosine 5'-triphosphate, there are no sequence specific contacts between the protein side chains and the nascent base pair, suggesting that steric constraints imposed by the protein onto the nascent base pair is the major contributor to nucleotide selectivity at the polymerase active site. The protein around the polymerase active site also shows plasticity, which may be responsible for the substrate diversity of the enzyme. Two conserved side chains, Q754 and R573, form hydrogen bonds with the N3 atom in the purine base and O2 atom in the pyrimidine base at the minor groove side of the base pair formed by the incorporated ddNMP and the corresponding template base in all the four ternary complexes. These hydrogen-bonding interactions may provide a means of detecting misincorporation at this position.     

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.