The HMG domain of lymphoid enhancer factor 1 bends DNA and facilitates assembly of functional nucleoprotein structures.

The high mobility group (HMG) domain is a DNA-binding motif that is associated with several eukaryotic regulatory proteins, including the lymphoid enhancer-binding factor LEF-1 and the testis-determining factor SRY. Here, we provide evidence that DNA binding by the HMG domain of LEF-1 involves primarily minor groove contacts and induces a bend of approximately 130 degrees in the DNA helix. Bending was also found to accompany sequence-specific DNA binding by the SRY-HMG domain. Examining possible regulatory roles of HMG domain-induced DNA bends, we found that LEF-1 can function in a manner similar to bacterial integration host factor and facilitate communication between widely separated protein-binding sites in a recombination assay. Together with the previous observation that LEF-1 by itself is unable to augment basal promoter activity, these data suggest that HMG domain proteins can serve as "architectural" elements in the assembly of higher-order nucleoprotein structures.     

Effect of G40R mutation on the binding of human SRY protein to DNA: a molecular dynamics view

Molecular dynamics simulation was conducted to investigate the reason why the mutant G40R of hSRY protein has a low affinity for DNA. Compared with the previous dynamics results of the wild-type hSRY-HMG-DNA complex, the results of molecular dynamics simulation on the mutant G40R hSRY-HMG-DNA system demonstrated that the whole structure of DNA (especially the second strand) had a major deviation away from the short arm of the HMG box. Consequently, the DNA and the mutant protein could not specifically recognize each other, that is, very different, and low-occupancy, direct, and water-mediated hydrogen bonds were detected at the protein-DNA interface, no conformational changes occurred at the loop region around Met9 during the simulation, and residue IIe13 did not intercalate between the bases of A5 and A6. These results indicated that the mutant G40R did not form a specific complex with the DNA target, hence led to complete gonadal dysgenesis. From the simulation, we realized that the residue Gly40 played a critical structural role in the hSRY-DNA recognition. It might be a structural supporting point of DNA binding because of the absence of a side chain. The reason for the difficulty of the mutant G40R to form a complex with DNA might be that the long and positively charged side chain of Arg40 by its bulk and positive charge hindered the DNA's access to the active sites of the protein.     

T antigen origin-binding domain of simian virus 40: determinants of specific DNA binding

To better understand origin recognition and initiation of DNA replication, we have examined by NMR complexes formed between the origin-binding domain of SV40 T antigen (T-ag-obd), the initiator protein of the SV40 virus, and cognate and noncognate DNA oligomers. The results reveal two structural effects associated with "origin-specific" binding that are absent in nonspecific DNA binding. The first is the formation of a hydrogen bond (H-bond) involving His 203, a residue that genetic studies have previously identified as crucial to both specific and nonspecific DNA binding in full-length T antigen. In free T-ag-obd, the side chain of His 203 has a pK(a) value of approximately 5, titrating to the N(epsilon)(1)H tautomer at neutral pH (Sudmeier, J. L., et al. (1996) J. Magn. Reson., Ser. B 113, 236-247). In complexes with origin DNA, His 203 N(delta)(1) becomes protonated and remains nontitrating as the imidazolium cation at all pH values from 4 to 8. The H-bonded N(delta1)H resonates at 15.9 ppm, an unusually large N-H proton chemical shift, of a magnitude previously observed only in the catalytic triad of serine proteases at low pH. The formation of this H-bond requires the middle G/C base pair of the recognition pentanucleotide, GAGGC. The second structural effect is a selective distortion of the A/T base pair characterized by a large (0.6 ppm) upfield chemical-shift change of its Watson-Crick proton, while nearby H-bonded protons remain relatively unaffected. The results indicate that T antigen, like many other DNA-binding proteins, may employ "catalytic" or "transition-state-like" interactions in binding its cognate DNA (Jen-Jacobson, L. (1997) Biopolymers 44, 153-180), which may be the solution to the well-known paradox between the relatively modest DNA-binding specificity exhibited by initiator proteins and the high specificity of initiation.     

DNA strand arrangement within the SfiI-DNA complex: atomic force microscopy analysis

The SfiI restriction enzyme binds to DNA as a tetramer holding two usually distant DNA recognition sites together before cleavage of the four DNA strands. To elucidate structural properties of the SfiI-DNA complex, atomic force microscopy (AFM) imaging of the complexes under noncleaving conditions (Ca2+ instead of Mg2+ in the reaction buffer) was performed. Intramolecular complexes formed by protein interaction between two binding sites in one DNA molecule (cis interaction) as well as complexes formed by the interaction of two sites in different molecules (trans interaction) were analyzed. Complexes were identified unambiguously by the presence of a tall spherical blob at the DNA intersections. To characterize the path of DNA within the complex, the angles between the DNA helices in the proximity of the complex were systematically analyzed. All the data show clear-cut bimodal distributions centered around peak values corresponding to 60 degrees and 120 degrees. To unambiguously distinguish between the crossed and bent models for the DNA orientation within the complex, DNA molecules with different arm lengths flanking the SfiI binding site were designed. The analysis of the AFM images for complexes of this type led to the conclusion that the DNA recognition sites within the complex are crossed. The angles of 60 degrees or 120 degrees between the DNA helices correspond to a complex in which one of the helices is flipped with respect to the orientation of the other. Complexes formed by five different recognition sequences (5'-GGCCNNNNNGGCC-3'), with different central base pairs, were also analyzed. Our results showed that complexes containing the two possible orientations of the helices were formed almost equally. This suggests no preferential orientation of the DNA cognate site within the complex, suggesting that the central part of the DNA binding site does not form strong sequence specific contacts with the protein.     

An asymmetric structure of the Bacillus subtilis replication terminator protein in complex with DNA

In Bacillus subtilis, the termination of DNA replication via polar fork arrest is effected by a specific protein:DNA complex formed between the replication terminator protein (RTP) and DNA terminator sites. We report the crystal structure of a replication terminator protein homologue (RTP.C110S) of B. subtilis in complex with the high affinity component of one of its cognate DNA termination sites, known as the TerI B-site, refined at 2.5 A resolution. The 21 bp RTP:DNA complex displays marked structural asymmetry in both the homodimeric protein and the DNA. This is in contrast to the previously reported complex formed with a symmetrical TerI B-site homologue. The induced asymmetry is consistent with the complex's solution properties as determined using NMR spectroscopy. Concomitant with this asymmetry is variation in the protein:DNA binding pattern for each of the subunits of the RTP homodimer. It is proposed that the asymmetric "wing" positions, as well as other asymmetrical features of the RTP:DNA complex, are critical for the cooperative binding that underlies the mechanism of polar fork arrest at the complete terminator site.     

Insights into the mechanisms of the selectivity filter of Escherichia coli aquaporin Z

Aquaporin Z (AQPZ) is a tetrameric protein that forms water channels in the cell membrane of Escherichia coli. The histidine residue (residue 174) in the selectivity filter (SF) region plays an important role in the transport of water across the membrane. In this work, we perform equilibrium molecular dynamics (MD) simulations to illustrate the gating mechanism of the SF and the influences of residue 174 in two different protonation states: Hsd174 with the proton at Nδ, and Hse174 with the proton at Nε. We calculate the pore radii in the SF region versus the simulation time. We perform steered MD to compute the free-energy profile, i.e., the potential of mean force (PMF) of a water molecule through the SF region. We conduct a quantum mechanics calculation of the binding energy of one water molecule with the residues in the SF region. The hydrogen bonds formed between the side chain of Hsd174 and the side chain of residue 189 (Arg189) play important roles in the selectivity mechanism of AQPZ. The radii of the pores, the hydrogen-bond analysis, and the free energies show that it is easier for water molecules to permeate through the SF region of AQPZ with residue 174 in the Hse state than in the Hsd state.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

The DNA bend angle and binding affinity of an HMG box increased by the presence of short terminal arms.

The HMG box of human LEF-1 (hLEF-1, formerly TCF1alpha) has been expressed in four forms: a parent box of 81 amino acids and constructs having either a 10 amino acid C-terminal extension, a 9 amino acid N-terminal extension, or both. These four species have been compared for DNA binding and bending ability using a 28 bp recognition sequence from the TCR alpha-chain enhancer. In the bending assay, whereas the parent box and that with the N-terminal extension bent the DNA by 57/58 degrees, the box extended at the C-terminus bent the DNA by 77/78 degrees, irrespective of the presence or absence of the N-terminal extension. A 6- fold increase in DNA affinity also resulted from addition of both terminal extensions. These observations redefine the functional boundaries of the HMG box. The structure of a mouse LEF-1/DNA complex recently published [Love et al. (1995) Nature 376, 791-795] implies that the higher DNA affinity and in particular the increased bend angle observed are consequences, at least in part, of the C-terminal extension spanning the major groove on the inside of the DNA bend.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

The intercalation of imidazoacridinones into DNA induces conformational changes in their side chain

Imidazoacridinones (IAs) are a new group of highly active antitumor compounds. The intercalation of the IA molecule into DNA is the preliminary step in the mode of action of these compounds. There are no experimental data about the structure of an intercalation complex formed by imidazoacridinones. Therefore the design of new potentially better compounds of this group should employ the molecular modelling techniques. The results of molecular dynamics simulations performed for four IA analogues are presented. Each of the compounds was studied in two systems: i) in water, and ii) in the intercalation complex with dodecamer duplex d(GCGCGCGCGCGC)2. Significant differences in the conformation of the side chain in the two environments were observed for all studied IAs. These changes were induced by electrostatic as well as van der Waals interactions between the intercalator and DNA. Moreover, the results showed that the geometry of the intercalation complex depends on: i) the chemical constitution of the side chain, and ii) the substituent in position 8 of the ring system.     

Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain in complex with DNA and free in solution.

Molecular dynamics simulations have been performed on the glucocorticoid receptor DNA binding domain (GR DBD) in aqueous solution as a dimer in complex with DNA and as a free monomer. In the simulated complex, we find a slightly increased bending of the DNA helix axis compared with the crystal structure in the spacer region of DNA between the two half-sites that are recognized by GR DBD. The bend is mainly caused by an increased number of interactions between DNA and the N-terminal extended region of the sequence specifically bound monomer. The recognition helices of GR DBD are pulled further into the DNA major groove leading to a weakening of the intrahelical hydrogen bonds in the middle of the helices. Many ordered water molecules with long residence times are found at the intermolecular interfaces of the complex. The hydrogen-bonding networks (including water bridges) on either side of the DNA major groove involve residues that are highly conserved within the family of nuclear receptors. Very similar hydrogen-bonding networks are found in the estrogen receptor (ER) DBD in complex with DNA, which suggests that this is a common feature for proper positioning of the recognition helix in ER DBD and GR DBD.     

Solution structure and base perturbation studies reveal a novel mode of alkylated base recognition by 3-methyladenine DNA glycosylase I

The specific recognition mechanisms of DNA repair glycosylases that remove cationic alkylpurine bases in DNA are not well understood partly due to the absence of structures of these enzymes with their cognate bases. Here we report the solution structure of 3-methyladenine DNA glycosylase I (TAG) in complex with its 3-methyladenine (3-MeA) cognate base, and we have used chemical perturbation of the base in combination with mutagenesis of the enzyme to evaluate the role of hydrogen bonding and pi-cation interactions in alkylated base recognition by this DNA repair enzyme. We find that TAG uses hydrogen bonding with heteroatoms on the base, van der Waals interactions with the 3-Me group, and conventional pi-pi stacking with a conserved Trp side chain to selectively bind neutral 3-MeA over the cationic form of the base. Discrimination against binding of the normal base adenine is derived from direct sensing of the 3-methyl group, leading to an induced-fit conformational change that engulfs the base in a box defined by five aromatic side chains. These findings indicate that base specific recognition by TAG does not involve strong pi-cation interactions, and suggest a novel mechanism for alkylated base recognition and removal.     

Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.