Structures of repressors bound to direct and inverted DNA repeats exp
http://www.100md.com
《核酸研究医学期刊》
Institut für Chemie und Biochemie/Kristallographie, Freie Universit?t Berlin Takustr. 6, 14195 Berlin, Germany 1Departamento de Biotecnología Microbiana, Centro Nacional de Biotecnología CSIC, 28049 Madrid, Spain
*To whom correspondence should be addressed. Tel: +49 30 838 53412; Fax: +49 30 838 56702; Email: saenger@chemie.fu-berlin.de
ABSTRACT
Repressor regulates transcription of genes required for copy number control, accurate segregation and stable maintenance of inc18 plasmids hosted by Gram-positive bacteria. belongs to homodimeric ribbon-helix-helix (RHH2) repressors typified by a central, antiparallel ?-sheet for DNA major groove binding. Homodimeric 2 binds cooperatively to promotors with 7 to 10 consecutive non-palindromic DNA heptad repeats (5'-A/TATCACA/T-3', symbolized by ) in palindromic inverted, converging () or diverging () orientation and also, unique to 2 and contrasting other RHH2 repressors, to non-palindromic direct () repeats. Here we investigate with crystal structures how 2 binds specifically to heptads in minimal operators with () and () repeats. Since the pseudo-2-fold axis relating the monomers in 2 passes the central C–G base pair of each heptad with 0.3 ? downstream offset, the separation between the pseudo-2-fold axes is exactly 7 bp in (), 0.6 ? shorter in () but would be 0.6 ? longer in (). These variations grade interactions between adjacent 2 and explain modulations in cooperative binding affinity of 2 to operators with different heptad orientations.
INTRODUCTION
Gene expression in prokaryotes is primarily regulated by helix–turn–helix proteins that bind specifically to palindromic operators whereas recognition of arrays of direct or inverted repeats by transcriptional, homodimeric ribbon-helix-helix (RHH2) repressors like protein is less frequent (1). Structures are known for RHH2 repressors Arc (2), CopG (3) and MetJ (4) bound to their cognate operators that are bent by 50° to 60°. MetJ2 binds symmetrically to two to five consecutive 8 bp long palindromic repeats. By contrast, CopG2 and Arc2 bind asymmetrically to half sites of palindromic operators that are spaced by 10 and 11 bp, respectively. When bound to these operators, interactions between adjacent RHH2 contribute to high affinity and cooperative association.
Repressor is a global regulator of and encoded by broad-host-range and low-copy number plasmids belonging to the inc18 family that are stably maintained in Gram-positive bacteria (5–7). was originally isolated from Streptococcus pyogenes plasmid pSM19035 where 2 controls promoter regions located upstream of genes involved in plasmid copy number control (PcopS), plasmid partitioning (P) and post-segregational killing (P) if the plasmid is lost. These promoters comprise arrays of ten, nine or seven consecutive 7 bp repeats (heptads, symbolized by ), organized as: PcopS, (); P, () and P, () (1), see Supplementary Figure 9.
Binding of 2 to a single heptad or to heptads separated by one or more additional base pair is poor (kD >500 nM), but tight if operators include at least two consecutive heptads and tightens further with increasing number of heptads. In addition, the affinity depends on heptad arrangement as shown by 6-fold reduced affinity of 2 for diverging repeats () (kD 120 nM) compared to heptads in direct () or converging () arrangement (kD 20 nM) (8). Multiple repeat binding sites are also found for eukaryotic operators that interact cooperatively with monomeric and therefore asymmetric transcription factors (1). However, these repeats show different base pair spacings and protein–protein interactions in direct and inverted orientation (9).
In wild type (wt) 2, the N-terminal 20 residues of the 71 residues long monomers are probably unstructured as suggested by secondary structure prediction (10) and were cleaved during crystallization (202). The structure features a typical RHH-fold comprising a 2-fold symmetrical ?-sheet with antiparallel pairing of residues 28–32 of each monomer followed by -helices 1 (34–46) and 2 (51–66) (11).
We describe here the crystal structures of an N-terminal deletion mutant (see Results) with 19 residues removed, hereafter 19, in complex with two minimal operators comprising two heptads in () and () arrangement. It was of interest to elucidate the structural determinants for high specificity, affinity and cooperative binding of repressor to minimal binding sites and to extrapolate these to natural operators with different heptad arrangements.
Unintentionally, both complexes cocrystallized with free operator DNA. Free ()-DNA allowed us to compare structural changes in DNA induced by repressor binding, whereas free ()-DNA was ill-defined in the electron density and could not be fully modeled. The complex between 192 and ()-DNA dissociated during gel filtration and could not be crystallized.
MATERIALS AND METHODS
Plasmid construction
For expression of 19 in Escherichia coli, wt gene missing the first 19 codons was cloned into NcoI–BamHI-cleaved pET28a (Novagen) to render pET28a-19. The described (1) pHP14-borne gene (pHP14) was modified to pHP1419 containing promoter P, the ribosomal binding site and the Met start codon fused to codon 20 of gene. pHP14 mutants pHP14Thr29Ala and pHP14His38Val were generated by site-directed mutagenesis. The plasmids were transferred into Bacillus subtilis strain BG511 (P:lacZ, recA4) as described (1).
?-galactosidase assay
?-galactosidase assays (Table 1) were performed as described (1) except that the centrifuged B.subtilis cells were resuspended and lysed by the addition of 0.1% SDS (final concentration 0.0025%) and chloroform.
Preparation of protein–DNA complexes
192 was expressed in E.coli according to (1), the cell paste was resuspended in buffer A and lysed (French Press). The crude extract was processed (1), except that after the phosphocellulose step the fractions were pooled, diluted 5-fold with buffer A and loaded on a POROS 20 HE column (Applied Biosystems). 192 was eluted with 50–1000 mM NaCl gradient in buffer A. Concentrated fractions were gel filtrated on Superdex75 (GE Healthcare) run with buffer B .
Complementary oligonucleotides were purified by high-performance liquid chromatography (HPLC), mixed at 1:1 molar ratio, hybridized and purified using a MonoQ column (GE Healthcare). Eluted DNA was dialyzed against 20 mM Tris–HCl (pH 7.5), 100 mM KCl and 30 mM NaCl. 192 was added at 2.1:1 molar ratio and purified by gelfiltration (SuperdexS75). Fractions of 192/DNA were concentrated to 10.5 mg/ml for crystallization.
Crystallization
Using hanging-drop vapor diffusion, crystals with space group C2 grew from drops made of 2 μl of 2-() solution or 2 μl of 2-() (bp A9–T6' exchanged by bp G9–C6') solution and 2–3 μl precipitant solution .
Crystals of 2-() grew in space group P21 under similar conditions when precipitant solution contained 150 mM Na/KPO4 (pH 7.0), 2.4 M Na2-malonate and 2 to 3% 2-methyl-2,4-pentanediol. In all cases crystal quality was improved by micro-seeding.
Data collection, structure determination and refinement
X-ray data were collected at 100 K at the Protein Structure Factory beamline BL14.1 of Free University Berlin at BESSY and processed with HKL2K (12); Table 2.
The structure of 2-() was determined by molecular replacement in PHASER (13) with 202 (PDB code 1IRQ) modeled to 8 bp idealized B-DNA. After manual building of 2-() and restrained refinement in REFMAC5 (14), Fo-Fc maps showed additional electron density for another DNA molecule but not for additional 192. This second DNA molecule was build manually starting from ideal B-form DNA. In the final model the asymmetric unit consists of one 2-() and one free ()-DNA. 19 molecules A' and B could be modeled with all residues 19–71, but A and B' only with residues 23–71 and 25–71, respectively, see Figures 1 and 2 for assignment of A, A', B, B'.
Figure 1 DNA used for cocrystallization and arrangement of (A) 2-() and (B) 2-() and free DNAs in the crystal asymmetric units (large boxes). Subunits sharing the same protein–DNA interactions 19A/B shown in blue and 19A'/B' in white ovals. Solid bars indicate subunits interacting through helices 1, 1'. The pseudo-2-fold axis relating 192 and bound heptads in 2-() indicated by black ellipses in (B). Heptads with sequence 5'-AATCACA/T-3' outlined by arrows, and nucleotides are numbered 0 to 17 and –1' to 16', respectively. Top strands of free DNA in green. Grey nucleotides C17 in both structures and 4 to 15 bp of free DNA in 2-() could not be modeled due to poor electron density.
Figure 2 (A) Structure of 2-(). DNA backbone trace in light grey for top strand and dark grey for bottom strand, see Figure 1. 19 monomers A/A' and B/B' in light and dark green and blue, respectively, helices 1 and 2 labeled with white letters. Helices 1' of A' and 1 of B related by a pseudo-2-fold axis perpendicular to the paper plane (red ellipse) form the 192···192 interface. (B) Superimposition of 19A/A' of 2-() (green) and 19A/A' of 2-() (grey) to show slight positional differences of 19 B/B' associated with palindromic symmetry in (). Pseudo-2-fold axes relating monomers in 192 indicated by dashed lines colored green for 2-() and red for 2-(). Helices 1, 1' involved in 192···192 interactions superimpose well allowing cooperative binding in both complexes. DNA in 2-() is locally kinked by 12° at the centre (bp G11–C4') of heptad A8–T14, and 192 are 0.6 ? closer (vertical distance between pseudo-2-fold axes) than in 2-(), see text. (C) Phosphate backbone of the 14 bp operator regions of operator DNA free () in green and 192-bound () in blue superimposed on ideal B-DNA (grey).
The structure determination of 2-() mutant with bp A9–T6' replaced by bp G9–C6' used difference Fourier technique applied to the isomorphous crystal structure of 2-(), see Table 2.
The structure of 2-() was determined in MOLREP (15) using 2-() as search model. One 2-() was found, and after restrained refinement sparse electron density indicated only four additional bases for free ()-DNA that could not be modeled completely, Figure 1B. Molecule B' could be modeled with all residues, B with residues 25–71 and A, A' with residues 22–71.
For refinement of all three structures in Refmac5 TLS groups were assigned and refined for each polypeptide chain and oligonucleotide, see Table 2 for statistics. No non-crystallographic symmetry was used during refinements. Model quality was examined by Whatcheck and Procheck (16) showing that , torsion angles of most amino acids in all three structures are within the most favored, some are in additionally allowed and none are in forbidden areas of the Ramachandran plot. Figures were generated with MOLSCRIPT (17) and Raster3D (18). Analysis of DNA parameters used program Curves (19).
RESULTS
192 protein
Cocrystallization of wt 2 with operator DNA yielded only crystals with poor X-ray diffraction, but was successful with 19. In vitro, 192 binds specifically to promoter PcopS with 2-fold lower affinity (kD 12 nM) (see Supplementary Figure 10) compared to wt 2 (kD 6 nM) (8), and likewise plasmid-borne 19 gene product represses P utilization in vivo 2-fold weaker compared to wt gene (Table 1). This suggests that even without the N-terminal 19 residues, 2 still binds strongly (only 2-fold weaker) and specifically to DNA heptads in both, in vitro and in vivo gene regulation.
Table 1 Utilization of P:lacZ in the presence of variants in B.subtilis cells
Crystal unit cells contain free and 192 bound operator DNA
The crystal structures of 192 bound to two minimal operators formed by 17 bp DNAs with C, G overhangs and comprising direct () and inverted () heptads (Figures 1 and 2) were determined by molecular replacement at 2.45 and 2.6 ? resolution, respectively (Table 2). The asymmetric units of both complexes contain two 192 bound to operator DNAs (2-() and 2-()) which in turn interact with the ends to free operator DNAs () and (), respectively, to form pseudo-continuous DNA (Figures 1 and 3).
Table 2 Crystallographic data and refinement statistics
Figure 3 2F0-Fc electron density distribution contoured at 1.2 level; view perpendicular to crystallographic b, c planes with the c-axes vertically oriented. (A) Parallel free DNAs (yellow) in 2-()() are at 5 ? shortest distance to each other and stack at both ends with 192-bound DNAs (magenta). Free DNA is not involved in any further crystal contacts. (B) Discontinous electron density for free DNA in 2-()(). Modeled complexes of 2-() and modeled bases of free DNA (see Figure 1B) in yellow. The central parts of free DNA marked by dotted lines (4 to 15 bp, Figure 1B) could not be modeled due to patchy electron density. The shortest distance between free DNAs is 8 ? and they are at an angle of 40°.
On the ‘left’ sides of free DNA (Figure 1), nucleotides C17 of free and 192-bound DNA are not in helical arrangement and not seen in the electron density as they are disordered but were confirmed by MALDI-TOF-spectrometry of dissolved crystals (data not shown). Both, 2-bound and free DNAs, stack with bp G16–C-1' that are related by pseudo-2-fold symmetry.
On the ‘right’ sides of free DNA (Figure 1), the two 3'-G16' overhangs lie in the minor groove of the adjacent duplex and interact with both 5'-G0 to form two consecutive G16'*(G0–C15') base-triplets (Supplementary Figure 11) with similar geometry as reported (20–22).
In the crystals of 192 bound to () in space group C2 and bound to () in space group P21 (Table 2), both 2–DNA complexes interact by protein–protein contacts to form layers parallel to the crystallographic a, b planes. The crystallographic a, b-axes in space group P21 correspond to b, a-axes in C2 thus reflecting space group and lattice packing similarities. The pseudo-continuous DNA helices (Figures 1 and 3) are oriented in c-direction. The doubled c-axis in C2 relative to that in P21 is due to the C-centering, and the DNA helices are parallel to the c-axis at shortest inter-helix distance of 5 ? in 2-()() (base-triplets indicated by ), possibly stabilized by bridging water. By contrast, in 2-()() the pseudo-continuous DNA helices are parallel to the crystallographic a, c plane but inclined at an angle of 40° towards the c-axis and at least 8 ? apart, and poor electron density (Figure 3B) indicates that they are partially disordered as shown by B factors >100 ?2. Consequently, bp 4 to 15 of free ()-DNA in 2-()() could not be modeled (Figures 1B and 3B).
In both crystal unit cells (space groups C2 and P21, Table 2), the DNA-bound 192 show minor structural changes compared to the X-ray structure of free 202 (11). This concerns the loop connecting -helices 1 and 2 and to a minor extent the ?-sheet (residues 27–32, 0.8 ? r.m.s. deviation).
Since gel permeation chromatography of 2-() at 20 mM Tris–HCl (pH 7.5), 150 mM NaCl indicated an apparent molecular mass of 40 kDa (calculated Mw = 33.4 kDa), we assume that under crystallization conditions with 2.4 M Na2-malonate, the complexes associate pairwise through base-triplets G16'*(G0–C15') to form the complex 2-()()-2. If one 2-() of this complex forms layers mediated by protein–protein contacts as in the present structures, it appears that the other 2-() has to release 2 for packing reasons, thereby giving rise to the crystallized 2-()(). The other crystal contact with only DNA–DNA stacking interactions is more flexible than the base-triplets and permits formation of a regular crystal lattice. The same crystallization scenario applies for 2-().
Protein–DNA interactions
Formation of G16'*(G0–C15') base-triplets induces distortions at the 5' ends of heptads A1–A7 as indicated by different, partly water-mediated interactions to 192 compared to heptads A8–A14 in () and A8–T14 in (), Figures 4 and 5. For this reason, we focus here on protein–DNA interactions for the less distorted dimer 19B/B' bound to heptad A8–A14 of 2-() (Figures 2A and 4(left)). In both structures, for each 192-bound heptad the direct (not water-mediated) protein–DNA contacts are comparable.
Figure 4 Schematic representation of interactions between 192 and DNA in 2-() (left) and 2-() (right). The orientation of heptads is indicated by blue arrows and palindromic symmetry by the black ellipse. Yellow bases specifically interact with 192. Residues labeled green for 19 A/B and black for 19 A'/B'. Hydrogen bonds to phosphate oxygens (red) in thin black lines, specific hydrogen bonds to bases in the major groove in black arrows, blue lines indicate water (blue dots) mediated hydrogen bonds. K28 (K28') and Arg31 (Arg31') interact differently in subunits A/B and A'/B'.
Figure 5 Stereo view of specific interactions between 19 B/B' and bound heptad A8–A14 in 2-(). Water molecules; red spheres, hydrogen bonds; dashed lines, phosphorus atoms; magenta, directly contacted bp C11–G4' and base G2' in yellow.
In the major grooves, base pair specific interactions are formed with Thr29 and Arg31 located on the ?-sheet. Thr29O and Thr29'O of 19B and B' bind specifically to the central bp G4'–C11, and Arg31N, hydrogen bond with base G2' (Figure 4(left)). In contrast, the corresponding Arg31' of 19B' hydrogen bonds with N' to Thr29O of 19B and with N1',2' through three water molecules to bases G4', A5' and A9, Figures 4 and 5.
To see whether 192 would bind symmetrically to a palindromic heptad featuring two G to provide both, Arg31 and Arg31', with potential binding partners, we determined the 2.9 ? resolution crystal structure of 2 in complex with a mutated () where bp A9–T6' was replaced by G9–C6' (heptad sequence: 8AGTCACA14, Table 2). However, the same interaction pattern was found as in the original heptad. This agrees with similar (2-fold weaker) binding of wt 2 to mutated () with the same replacement in the first heptad (1AGTCACA7) compared to the original operator (8). A referee suggested to test the binding of 192 to heptads with pseudo-palindromic symmetry, 5'-TGTCACA-3'. In view of the binding geometry of 192 to 8AGTCACA14, we question whether this would provide novel knowledge. This is because the 5'-AAT- or 5'-AGT-ends of the heptads contact 192 in all cases exclusively through unspecific interactions with phosphate groups or are mediated by water molecules (Figure 4A and B). Hence, base pair exchanges in this part of the heptads should not significantly affect binding of 192.
Backbone phosphates of all four heptads contact helices 1 and 2 of 19 subunits with pseudo-2-fold symmetry (Figure 4). In 19B, the 5'-phosphate of A9 caps the N-terminus of 2 by hydrogen bonding to peptide amides of V51 and K52 in a pattern known for RHH2 proteins and other repressors (4,23), and His37N, K41N of 1 bind to the 5'-phosphate of T10. Corresponding residues of 19B' interact with 5'-phosphates of G2' and T3'. However, K28N located on the ?-sheet of 19B forms a salt-bridge with the 5'-phosphate of G4' whereas K28'N on 19B' is >5 ? away from the corresponding 5'-phosphate of C11.
The asymmetry in the binding of each 192 to its particular heptad is reflected by superimposition of C atoms of monomers (A on A' and B on B', Figure 2A) in each 192 showing 0.6 ? r.m.s. deviation partly associated with structural differences in the loops connecting helices 1 and 2 (residues 46–48) and in the ?-strands. When 192 dimers are superimposed on each other in the same orientation (A on B and A' on B'), i.e. the bound heptads point in the same direction, r.m.s. deviation of only 0.3 ? confirms that all 192-heptad interactions are similar.
Thr29 is essential for specific operator binding
To test the importance of Thr29 for heptad sequence recognition in vivo studies were conducted showing that 2Thr29Ala failed to completely repress promoter P utilization (Table 1). Binding of 2Thr29Ala to PcopS operator DNA embedded in 300 bp DNA was tested by electrophoretic mobility shift assays (EMSA) (Supplementary Figure 10). Protein–DNA complexes formed with 192 but Thr29Ala required 100-fold higher concentration (compared to 192) that yielded prominent but unspecific binding to PcopS as confirmed by DNase I footprinting (data not shown).
192···192 interactions
Protein–DNA interactions in both complexes bury 1610 ?2 of solvent accessible surface area. Another 550 ?2 are buried by interaction of pseudo-2-fold axis related 1 helices of adjacent monomers A' and B in 2-() and A' and B' in 2-(), respectively (Figure 1 and Figure 2A and B). In detail, bifurcated hydrogen bonds between His38N and Ala45'O/Lys46'O are augmented by hydrophobic contacts between Ile42, Ala45 of both subunits (Figure 6). Positions of His38 and His38' are identical in both complexes whereas the hydrophobic side-chain of Ile42 adopts different rotamers without affecting the size of the buried interface. These interactions ensure cooperative binding when several 2 associate with multiple heptad repeats as found in natural operators.
Figure 6 Interactions at the 192···192 interface (A) between helices 1' and 1 of subunits A' and B of 2-(), (see Figure 2A). The 192 are related by a pseudo-2-fold axis perpendicular to the paper plane (red ellipse). Residues with contacts <4 ? to neighbor helix are depicted with side-chains, colors as in Figure 2. Bifurcated hydrogen bonds between His38 (His38') and Ala45O/Lys46O to Ala45'O/Lys46'O of the other subunit in dashed lines. (B) Superimposition as in Figure 2B with view on the dimer–dimer interface to show that dimer–dimer interactions are comparable in both structures.
To test the importance of His38, plasmid-borne His38Val was constructed to remove the bifurcated hydrogen bonds but to maintain the hydrophobic character of the interface. 2His38Val repressed P utilization in vivo with 7-fold lower efficiency than wt 2 (Table 1), indicating the important role of His38 for cooperative binding between 2 and multiple consecutive heptads.
Conformation of free and 192-bound DNA
The DNAs in 2-() and 2-() are nearly straight B-form with average helical twist of 36° (range 23° to 43°; Supplementary Figure 12), but show distinct features. The major groove width in 2-() shows strong modulation depending on nucleotide sequence, being 13 to 14 ? except for the 7AAAT10 tract (11 to 12 ?) compared to 11.7 ? for ideal B-DNA (24). In contrast, the major groove width in 2-() is more continuous (13 to 14 ?) due to the palindromic symmetry (Figure 7A). In the 7AAAT10 tract of 2-() and in the 7AAGT10 tract of 2-() the minor groove is narrowed to 2.7 and 4 ?, respectively, compared to 5.7 ? in ideal B-DNA (24) due to an average negative base pair inclination of –8° (Figure 7C), as frequently found for A-tracts (25). This is associated in 7AAAT10 of 2-() with binding of a single spine of four water molecules, one in each base step as reported in (26). Strong opposite buckles of the C–G bp in the 4CAC6 (–15°, 15°) and 11CAC13 (–10°, 15°) segments of 2-() (Figure 7B) widen the major groove to 13 ? to accommodate the ?-sheet, and adjacent (up- and downstream) minor grooves are narrowed (Figure 7A). Similarly, in 2-() the C–G bp in segment 4CAC6 show comparable buckles of –10°, 19° but in 9GTG11 the buckles are reduced to –8°, 5°. Of the four heptads, three feature two subsequent CA steps that are known for their ability to bend B-DNA through positive roll (20,27,28). However, the CA steps do not exhibit unusual structure except for a moderate slide movement of 1 ? at the central C–G bp where 192 binds to both bases (Supplementary Figure 12).
Figure 7 Comparison of operator DNA free and in complex with 192. (A) Minor (m) and major (M) groove widths along the central 11 bp of each operator, as measured by shortest P-P distances (less two phosphate group radii, 5.8 ?) indicated by horizontal dotted black lines for ideal B-DNA (24), 5.7 ? and 11.7 ?, respectively, Note minor groove narrowing in the 7AAAT10 tract of in 192-bound and free () DNA and reduced narrowing in 7AAGT10 tract of (192)2-(). (B) Base pair buckle (°) and (C) inclination angles (°) of each operator.
In the crystal lattice, free ()-DNA contacts 192-bound ()-DNA with both ends to form pseudo-continuous helices (Figure 1) and is not distorted by any further crystal contacts (Figure 3A). Interestingly, free ()-DNA shows comparable structural features as found for 192-bound ()-DNA (Figure 2C). In the 7AAAT10 tract of free ()-DNA, base pair inclination angles are negative, and the minor groove is narrowed to 3.5 ? whereas the major groove is overall widened to an average of 12.5 ?. Additionally, strong opposite buckles are observed for C–G base pair in both CAC segments that are comparable to 192-bound ()-DNA (Figure 7B). These observations suggest that the described significant deviations from ideal B-DNA are a consequence of nucleotide sequence and not induced by binding to 192 (Figures 2C and 7).
However, conformational changes upon repressor binding concern base pair parameters helical twist, roll and slide (Supplementary Figure 12). In free ()-DNA these parameters cluster around values assigned to ideal B-DNA, but show a negative roll of –6° and positive slide of 1 ? for the central G11'–C4 and G4'–C11 bp of both heptads, and helical twist angles for bp A5–T10' and A12–T3' decrease to 25° in 192-bound ()-DNA. Since these deviations are similar for both heptads of 2-() they can be attributed to 192 binding. Such comparison is not possible with free ()-DNA as it could not be modeled completely.
Comparison of 2-() and 2-()
The pseudo-2-fold axis relating the monomers in 192 passes with 0.3 ? downstream (5'3') offset through the central G–C bp of each heptad. Consequently, these symmetry axes are separated by 7 bp in 2-() but they are 0.6 ? closer in 2-(), see Figure 2B. Despite the different dimer–dimer separations helices 1 and 1' forming the 192...192 interfaces superimpose well in both complexes, consistent with similar dissociation constants (kD 20 nM). In contrast, the separation between the two 192 will be 0.6 ? wider in diverging heptads 2-(). Assuming that the interaction pattern between 192 and heptads with () orientation is similar as with () and () heptad orientations, the expected 0.6 ? longer 1···1 contacts are probably less favorable and diminish cooperativity. This agrees with the 6-fold weaker affinity of wt 2 to heptads in () arrangement, the finding that 2-() dissociated during gel filtration, and the drastically reduced binding affinity to heptad repeats spaced by one or more additional base pair (8).
Whereas DNA in 2-() is nearly straight, the superimposition in Figure 2B reveals a kink of 12° at the G11–C4' bp at the centre of the heptad bound to dimer B/B' of 2-(). This kink is associated with 12° rotation of dimer B/B' of 2-(), the rotational pivot point being located in the 192···192 interface. By virtue of these two motions, the position and orientation of helix 1' of subunit B' remain almost as in 2-() and similar 192···192 interactions explain comparable dissociation constants of both complexes.
Structural model of 2-bound to natural operators
Extrapolation of the structures of 2-() and 2-() allowed modeling of 192 in complex with natural P promoter (Figure 8). The model implies that wt 2 binds as left-handed matrix to right-handed, straight B-type operator DNA, each 192 being displaced relative to its neighbor by 7 bp and rotated by 252°. Figure 8B shows that the negatively charged sugar-phosphate backbone of DNA faces positively charged surface of 192.
Figure 8 Model of seven 192-bound heptads () of the natural promoter P. (A) DNA in space filling and 192 as orange/red ribbons, (B) DNA in stick and 192 in surface representation colored according to electrostatic potential (negative and positive charges red and blue, respectively). The model is based on the structures of 2-() and 2-(). Repressors form a left-handed protein-matrix winding around the nearly straight operator.
DISCUSSION
Implications of 192-DNA structures for regulation of transcription
The inc18 family plasmids harbor genes to control their copy number, accurate segregation and stable maintenance during cell division. Since expression of these genes is regulated by the common 2 repressor, a unique mechanism has evolved to fine-tune repressor affinity for the different operators. How this is achieved is shown by the present study. It clearly indicates that the pseudo-symmetric 2 binds with 0.3 ? downstream (5'3') offset relative to the center G–C base pair of the cognate heptad. Since the operators are nearly straight B-DNA, different heptad numbers and orientations lead to different distances between 1 helices of adjacent 2, thereby modulating cooperative interactions between 2 and different operators. The ability to bind to palindromic as well as to non-palindromic operators is a unique feature of 2 and is not shared by other member of the RHH2 family. We associate this with the interactions between 192 that are related by a pseudo-2-fold rotation axis (Figure 2A and Figure 6A and B) so that they interact comparably, no matter what the orientations of adjacent heptads are, (), () or (). It is unlikely that the deleted N-termini would contribute to 192 ···192 interactions as the N-termini of 192 point away from the 192 ···192 interface (Figures 2A and 6A).
When RHH2 bind to DNA, we have to consider two different binding characteristics, the ‘variable’ binding of flexible side-chains of ?-sheets and 1 helices to major groove and phosphate groups, respectively, and ‘stiff’ capping of N-termini of 2 helices by rigid main-chain NH hydrogen bonding to phosphates. Although Arc2, CopG2 and MetJ2 bend minimal cognate operators with two repressor binding sites by 50 to 60°, their cores do not rearrange significantly compared to free repressors and are structurally similar . In contrast, Arc2, CopG2 and MetJ2 superimpose on DNA-bound 192 with higher r.m.s. deviations of 2.3 to 2.8 ? because the ?-sheet of 192 protrudes 3 ? less from the repressor surface and is closer to the N-termini of helices 2 (Supplementary Figure 13). Due to this geometry of 192, the N-termini of helices 2 are in the correct position to clamp the phosphate backbones of straight operator DNA when the ?-sheet is inserted into the major groove (Figures 2 and 4). In contrast, the other three repressors have to bend DNA around the more protruding ?-sheet to place the phosphate backbone in hydrogen bonding distance to N-H groups of the N-termini of their helices 2. It is notable that the distances between phosphates bound by 2 helices is 5 bp in 2 (Figure 4) but 6 bp in the other three RHH2 since DNA has to follow a longer path when bending around the protruding ?-sheet.
The N-termini of RHH2 proteins are associated with different functions
The N-terminus in CopG2 has no obvious function, in MetJ2 it is involved in binding the corepressor S-adenosylmethionine whereas in free Arc2 it is disordered but forms a 310 helix upon and contributes to DNA-binding. The N-terminal residues in 2 do not contribute to DNA-binding in vitro and in vivo as shown in the present work. However, we have recently shown that protein , a homolog to ParA proteins and involved in active plasmid partitioning during cell division, is activated by wt 2 but not by 192 (A. Cicek, F. Pratto, W. Weihofen, J.C. Alonso and W. Saenger, unpublished data). This suggests that 2 is the yet missing ParB protein of inc18 family plasmids in the known ParA/ParB family of plasmid partitioning systems (29).
Direct and indirect readout of DNA sequence
Cocrystallized 192-bound and free ()-DNA show similar phosphate backbone conformation (Figure 2C) with significant deviations from ideal B-DNA (Figure 7), indicating that conformation is predominantly dependent on DNA sequence and not induced by 192 binding. This provides a good example for the ‘indirect’ readout of local DNA conformation by 192 that depends on the particular nucleotide sequence (30). The combination with ‘direct’ readout resulting from interactions of 2 repressor amino acids with heptad bases increases protein–DNA-binding specificity and affinity. This view is consistent with studies in which the heptad nucleotide sequence was mutated base pair by base pair (8). Mutations at heptad positions 2, 3 and 5 (5'-AATCACA-3') resulted in at least 4-fold weaker binding to wt 2, although these base pair are not directly contacted by repressor side-chains contrasting direct read out positions 4 and 6 that were much more sensitive against mutations.
Protein Data Bank accession codes
Atomic coordinates and structure factors have been deposited with the Protein Data Bank under accession codes 2bnw, 2bnz and 2cax.
SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.
ACKNOWLEDGEMENTS
The authors acknowledge beamtime and technical support at the Protein Structure Factory beamlines of Free University Berlin at BESSY/Berlin. This work was partially supported by Fonds der Chemischen Industrie to W.S., BMC2003-00150 from Dirección General de Investigación, Ministerio de Educación y Ciencia to J.C.A., and QLK3-CT-2001-00277 from the European Community to J.C.A. and W.S. Funding to pay the Open Access publication charges for this article was provided by Free University Berlin.
REFERENCES
de la Hoz, A.B., Ayora, S., Sitkiewicz, I., Fernandez, S., Pankiewicz, R., Alonso, J.C., Ceglowski, P. (2000) Plasmid copy-number control and better-than-random segregation genes of pSM19035 share a common regulator Proc. Natl Acad. Sci. USA, 97, 728–733 .
Raumann, B.E., Rould, M.A., Pabo, C.O., Sauer, R.T. (1994) DNA recognition by ?-sheets in the Arc repressor-operator crystal structure Nature, 367, 754–757 .
Gomis-Ruth, F.X., Sola, M., Acebo, P., Parraga, A., Guasch, A., Eritja, R., Gonzalez, A., Espinosa, M., del Solar, G., Coll, M. (1998) The structure of plasmid-encoded transcriptional repressor CopG unliganded and bound to its operator EMBO J, . 17, 7404–7415 .
Somers, W.S. and Phillips, S.E. (1992) Crystal structure of the met repressor-operator complex at 2.8 A resolution reveals DNA recognition by ?-strands Nature, 359, 387–393 .
Brantl, S. (1994) The copR gene product of plasmid pIP501 acts as a transcriptional repressor at the essential repR promoter Mol. Microbiol, . 14, 473–483 .
Pujol, C., Chedin, F., Ehrlich, S.D., Janniere, L. (1998) Inhibition of a naturally occurring rolling-circle replicon in derivatives of the -replicating plasmid pIP501 Mol. Microbiol, . 29, 709–718 .
Ceglowski, P., Boitsov, A., Chai, S., Alonso, J.C. (1993) Analysis of the stabilization system of pSM19035-derived plasmid pBT233 in Bacillus subtilis Gene, 136, 1–12 .
de la Hoz, A.B., Pratto, F., Misselwitz, R., Speck, C., Weihofen, W., Welfle, K., Saenger, W., Welfle, H., Alonso, J.C. (2004) Recognition of DNA by protein from the broad-host range Streptococcus pyogenes plasmid pSM19035: analysis of binding to operator DNA with one to four heptad repeats Nucleic Acids Res, . 32, 3136–3147 .
Glass, C.K. (1994) Differential recognition of target genes by nuclear receptor monomers, dimers, and heterodimers, Endocr. Rev, . 15, 391–407 .
Misselwitz, R., de la Hoz, A.B., Ayora, S., Welfle, K., Behlke, J., Murayama, K., Saenger, W., Alonso, J.C., Welfle, H. (2001) Stability and DNA-binding properties of the regulator protein from the broad-host range Streptococcus pyogenes plasmid pSM19035 FEBS Lett, . 505, 436–440 .
Murayama, K., Orth, P., de la Hoz, A.B., Alonso, J.C., Saenger, W. (2001) Crystal structure of transcriptional repressor encoded by Streptococcus pyogenes plasmid pSM19035 at 1.5 A resolution J. Mol. Biol, . 314, 789–796 .
Otwinowski, Z. and Minor, W. (1997) Processing of X-ray diffraction data collected in oscillation mode In Carter, C.W., Jr and Sweet, R.M. (Eds.). Macromolecular Crystallography Part A, Academic Press pp. 307–326 .
Storoni, L.C., McCoy, A.J., Read, R.J. (2004) Likelihood-enhanced fast rotation functions Acta Crystallogr. D Biol. Crystallogr, . 60, 432–438 .
Murshudov, G.N., Vagin, A.A., Dodson, E.J. (1997) Refinement of macromolecular structures by the maximum- likelihood method Acta. Crystallogr. D Biol. Crystallogr, . 53, 240–255 .
Vagin, A. and Teplyakov, A. (1997) MOLREP: an automated program for molecular replacement J. Appl. Cryst, . 30, 1022–1025 .
Laskowski, R.A., MacArthur, M.W., Moss, D.S., Thornton, J.M. (1993) PROCHECK: a program to check the stereochemical quality of protein structures J. Appl. Cryst, . 26, 283–291 .
Kraulis, P.J. (1991) MOLSCRIPT: a program to produce both detailed and schematic plots of protein structures J. Appl. Cryst, . 24, 946–950 .
Merritt, E.A. and Bacon, D.J. (1997) Raster3D: photorealistic molecular graphics Meth. Enzymol, . 277, 505–524 .
Lavery, R. and Sklenar, H. (1988) The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids J. Biomol. Struct. Dyn, . 6, 63–91 .
Schultz, S.C., Shields, G.C., Steitz, T.A. (1991) Crystal structure of a CAP–DNA complex: the DNA is bent by 90 degrees Science, 253, 1001–1007 .
Spink, N., Nunn, C.M., Vojtechovsky, J., Berman, H.M., Neidle, S. (1995) Crystal structure of a DNA decamer showing a novel pseudo four-way helix–helix junction Proc. Natl Acad. Sci. USA, 92, 10767–10771 .
Nunn, C.M., Garman, E., Neidle, S. (1997) Crystal structure of the DNA decamer d(CGCAATTGCG) complexed with the minor groove binding drug netropsin Biochemistry, 36, 4792–4799 .
Aggarwal, A.K., Rodgers, D.W., Drottar, M., Ptashne, M., Harrison, S.C. (1988) Recognition of a DNA operator by the repressor of phage 434: a view at high resolution Science, 242, 899–907 .
Saenger, W. Principles of Nucleic Acid Structure, (1984) NY Springer-Verlag .
Shatzky-Schwartz, M., Arbuckle, N.D., Eisenstein, M., Rabinovich, D., Bareket-Samish, A., Haran, T.E., Luisi, B.F., Shakked, Z. (1997) X-ray and solution studies of DNA oligomers and implications for the structural basis of A-tract-dependent curvature J Mol. Biol, . 267, 595–623 .
Han, G.W., Kopka, M.L., Cascio, D., Grzeskowiak, K., Dickerson, R.E. (1997) Structure of a DNA analog of the primer for HIV-1 RT second strand synthesis J. Mol. Biol, . 269, 811–826 .
Chen, S., Gunasekera, A., Zhang, X., Kunkel, T.A., Ebright, R.H., Berman, H.M. (2001) Indirect readout of DNA sequence at the primary-kink site in the CAP–DNA complex: alteration of DNA binding specificity through alteration of DNA kinking J. Mol. Biol, . 314, 75–82 .
El Hassan, M.A. and Calladine, C.R. (1998) Two distinct modes of protein-induced bending in DNA J.Mol. Biol, . 282, 331–343 .
Pogliano, J. (2002) Dynamic cellular location of bacterial plasmids Curr. Opin. Microbiol, . 5, 586–590 .
Otwinowski, Z., Schevitz, R.W., Zhang, R.G., Lawson, C.L., Joachimiak, A., Marmorstein, R.Q., Luisi, B.F., Sigler, P.B. (1988) Crystal structure of trp repressor/operator complex at atomic resolution Nature, 335, 321–329 .(Wilhelm Andreas Weihofen, Aslan Cicek, F)
*To whom correspondence should be addressed. Tel: +49 30 838 53412; Fax: +49 30 838 56702; Email: saenger@chemie.fu-berlin.de
ABSTRACT
Repressor regulates transcription of genes required for copy number control, accurate segregation and stable maintenance of inc18 plasmids hosted by Gram-positive bacteria. belongs to homodimeric ribbon-helix-helix (RHH2) repressors typified by a central, antiparallel ?-sheet for DNA major groove binding. Homodimeric 2 binds cooperatively to promotors with 7 to 10 consecutive non-palindromic DNA heptad repeats (5'-A/TATCACA/T-3', symbolized by ) in palindromic inverted, converging () or diverging () orientation and also, unique to 2 and contrasting other RHH2 repressors, to non-palindromic direct () repeats. Here we investigate with crystal structures how 2 binds specifically to heptads in minimal operators with () and () repeats. Since the pseudo-2-fold axis relating the monomers in 2 passes the central C–G base pair of each heptad with 0.3 ? downstream offset, the separation between the pseudo-2-fold axes is exactly 7 bp in (), 0.6 ? shorter in () but would be 0.6 ? longer in (). These variations grade interactions between adjacent 2 and explain modulations in cooperative binding affinity of 2 to operators with different heptad orientations.
INTRODUCTION
Gene expression in prokaryotes is primarily regulated by helix–turn–helix proteins that bind specifically to palindromic operators whereas recognition of arrays of direct or inverted repeats by transcriptional, homodimeric ribbon-helix-helix (RHH2) repressors like protein is less frequent (1). Structures are known for RHH2 repressors Arc (2), CopG (3) and MetJ (4) bound to their cognate operators that are bent by 50° to 60°. MetJ2 binds symmetrically to two to five consecutive 8 bp long palindromic repeats. By contrast, CopG2 and Arc2 bind asymmetrically to half sites of palindromic operators that are spaced by 10 and 11 bp, respectively. When bound to these operators, interactions between adjacent RHH2 contribute to high affinity and cooperative association.
Repressor is a global regulator of and encoded by broad-host-range and low-copy number plasmids belonging to the inc18 family that are stably maintained in Gram-positive bacteria (5–7). was originally isolated from Streptococcus pyogenes plasmid pSM19035 where 2 controls promoter regions located upstream of genes involved in plasmid copy number control (PcopS), plasmid partitioning (P) and post-segregational killing (P) if the plasmid is lost. These promoters comprise arrays of ten, nine or seven consecutive 7 bp repeats (heptads, symbolized by ), organized as: PcopS, (); P, () and P, () (1), see Supplementary Figure 9.
Binding of 2 to a single heptad or to heptads separated by one or more additional base pair is poor (kD >500 nM), but tight if operators include at least two consecutive heptads and tightens further with increasing number of heptads. In addition, the affinity depends on heptad arrangement as shown by 6-fold reduced affinity of 2 for diverging repeats () (kD 120 nM) compared to heptads in direct () or converging () arrangement (kD 20 nM) (8). Multiple repeat binding sites are also found for eukaryotic operators that interact cooperatively with monomeric and therefore asymmetric transcription factors (1). However, these repeats show different base pair spacings and protein–protein interactions in direct and inverted orientation (9).
In wild type (wt) 2, the N-terminal 20 residues of the 71 residues long monomers are probably unstructured as suggested by secondary structure prediction (10) and were cleaved during crystallization (202). The structure features a typical RHH-fold comprising a 2-fold symmetrical ?-sheet with antiparallel pairing of residues 28–32 of each monomer followed by -helices 1 (34–46) and 2 (51–66) (11).
We describe here the crystal structures of an N-terminal deletion mutant (see Results) with 19 residues removed, hereafter 19, in complex with two minimal operators comprising two heptads in () and () arrangement. It was of interest to elucidate the structural determinants for high specificity, affinity and cooperative binding of repressor to minimal binding sites and to extrapolate these to natural operators with different heptad arrangements.
Unintentionally, both complexes cocrystallized with free operator DNA. Free ()-DNA allowed us to compare structural changes in DNA induced by repressor binding, whereas free ()-DNA was ill-defined in the electron density and could not be fully modeled. The complex between 192 and ()-DNA dissociated during gel filtration and could not be crystallized.
MATERIALS AND METHODS
Plasmid construction
For expression of 19 in Escherichia coli, wt gene missing the first 19 codons was cloned into NcoI–BamHI-cleaved pET28a (Novagen) to render pET28a-19. The described (1) pHP14-borne gene (pHP14) was modified to pHP1419 containing promoter P, the ribosomal binding site and the Met start codon fused to codon 20 of gene. pHP14 mutants pHP14Thr29Ala and pHP14His38Val were generated by site-directed mutagenesis. The plasmids were transferred into Bacillus subtilis strain BG511 (P:lacZ, recA4) as described (1).
?-galactosidase assay
?-galactosidase assays (Table 1) were performed as described (1) except that the centrifuged B.subtilis cells were resuspended and lysed by the addition of 0.1% SDS (final concentration 0.0025%) and chloroform.
Preparation of protein–DNA complexes
192 was expressed in E.coli according to (1), the cell paste was resuspended in buffer A and lysed (French Press). The crude extract was processed (1), except that after the phosphocellulose step the fractions were pooled, diluted 5-fold with buffer A and loaded on a POROS 20 HE column (Applied Biosystems). 192 was eluted with 50–1000 mM NaCl gradient in buffer A. Concentrated fractions were gel filtrated on Superdex75 (GE Healthcare) run with buffer B .
Complementary oligonucleotides were purified by high-performance liquid chromatography (HPLC), mixed at 1:1 molar ratio, hybridized and purified using a MonoQ column (GE Healthcare). Eluted DNA was dialyzed against 20 mM Tris–HCl (pH 7.5), 100 mM KCl and 30 mM NaCl. 192 was added at 2.1:1 molar ratio and purified by gelfiltration (SuperdexS75). Fractions of 192/DNA were concentrated to 10.5 mg/ml for crystallization.
Crystallization
Using hanging-drop vapor diffusion, crystals with space group C2 grew from drops made of 2 μl of 2-() solution or 2 μl of 2-() (bp A9–T6' exchanged by bp G9–C6') solution and 2–3 μl precipitant solution .
Crystals of 2-() grew in space group P21 under similar conditions when precipitant solution contained 150 mM Na/KPO4 (pH 7.0), 2.4 M Na2-malonate and 2 to 3% 2-methyl-2,4-pentanediol. In all cases crystal quality was improved by micro-seeding.
Data collection, structure determination and refinement
X-ray data were collected at 100 K at the Protein Structure Factory beamline BL14.1 of Free University Berlin at BESSY and processed with HKL2K (12); Table 2.
The structure of 2-() was determined by molecular replacement in PHASER (13) with 202 (PDB code 1IRQ) modeled to 8 bp idealized B-DNA. After manual building of 2-() and restrained refinement in REFMAC5 (14), Fo-Fc maps showed additional electron density for another DNA molecule but not for additional 192. This second DNA molecule was build manually starting from ideal B-form DNA. In the final model the asymmetric unit consists of one 2-() and one free ()-DNA. 19 molecules A' and B could be modeled with all residues 19–71, but A and B' only with residues 23–71 and 25–71, respectively, see Figures 1 and 2 for assignment of A, A', B, B'.
Figure 1 DNA used for cocrystallization and arrangement of (A) 2-() and (B) 2-() and free DNAs in the crystal asymmetric units (large boxes). Subunits sharing the same protein–DNA interactions 19A/B shown in blue and 19A'/B' in white ovals. Solid bars indicate subunits interacting through helices 1, 1'. The pseudo-2-fold axis relating 192 and bound heptads in 2-() indicated by black ellipses in (B). Heptads with sequence 5'-AATCACA/T-3' outlined by arrows, and nucleotides are numbered 0 to 17 and –1' to 16', respectively. Top strands of free DNA in green. Grey nucleotides C17 in both structures and 4 to 15 bp of free DNA in 2-() could not be modeled due to poor electron density.
Figure 2 (A) Structure of 2-(). DNA backbone trace in light grey for top strand and dark grey for bottom strand, see Figure 1. 19 monomers A/A' and B/B' in light and dark green and blue, respectively, helices 1 and 2 labeled with white letters. Helices 1' of A' and 1 of B related by a pseudo-2-fold axis perpendicular to the paper plane (red ellipse) form the 192···192 interface. (B) Superimposition of 19A/A' of 2-() (green) and 19A/A' of 2-() (grey) to show slight positional differences of 19 B/B' associated with palindromic symmetry in (). Pseudo-2-fold axes relating monomers in 192 indicated by dashed lines colored green for 2-() and red for 2-(). Helices 1, 1' involved in 192···192 interactions superimpose well allowing cooperative binding in both complexes. DNA in 2-() is locally kinked by 12° at the centre (bp G11–C4') of heptad A8–T14, and 192 are 0.6 ? closer (vertical distance between pseudo-2-fold axes) than in 2-(), see text. (C) Phosphate backbone of the 14 bp operator regions of operator DNA free () in green and 192-bound () in blue superimposed on ideal B-DNA (grey).
The structure determination of 2-() mutant with bp A9–T6' replaced by bp G9–C6' used difference Fourier technique applied to the isomorphous crystal structure of 2-(), see Table 2.
The structure of 2-() was determined in MOLREP (15) using 2-() as search model. One 2-() was found, and after restrained refinement sparse electron density indicated only four additional bases for free ()-DNA that could not be modeled completely, Figure 1B. Molecule B' could be modeled with all residues, B with residues 25–71 and A, A' with residues 22–71.
For refinement of all three structures in Refmac5 TLS groups were assigned and refined for each polypeptide chain and oligonucleotide, see Table 2 for statistics. No non-crystallographic symmetry was used during refinements. Model quality was examined by Whatcheck and Procheck (16) showing that , torsion angles of most amino acids in all three structures are within the most favored, some are in additionally allowed and none are in forbidden areas of the Ramachandran plot. Figures were generated with MOLSCRIPT (17) and Raster3D (18). Analysis of DNA parameters used program Curves (19).
RESULTS
192 protein
Cocrystallization of wt 2 with operator DNA yielded only crystals with poor X-ray diffraction, but was successful with 19. In vitro, 192 binds specifically to promoter PcopS with 2-fold lower affinity (kD 12 nM) (see Supplementary Figure 10) compared to wt 2 (kD 6 nM) (8), and likewise plasmid-borne 19 gene product represses P utilization in vivo 2-fold weaker compared to wt gene (Table 1). This suggests that even without the N-terminal 19 residues, 2 still binds strongly (only 2-fold weaker) and specifically to DNA heptads in both, in vitro and in vivo gene regulation.
Table 1 Utilization of P:lacZ in the presence of variants in B.subtilis cells
Crystal unit cells contain free and 192 bound operator DNA
The crystal structures of 192 bound to two minimal operators formed by 17 bp DNAs with C, G overhangs and comprising direct () and inverted () heptads (Figures 1 and 2) were determined by molecular replacement at 2.45 and 2.6 ? resolution, respectively (Table 2). The asymmetric units of both complexes contain two 192 bound to operator DNAs (2-() and 2-()) which in turn interact with the ends to free operator DNAs () and (), respectively, to form pseudo-continuous DNA (Figures 1 and 3).
Table 2 Crystallographic data and refinement statistics
Figure 3 2F0-Fc electron density distribution contoured at 1.2 level; view perpendicular to crystallographic b, c planes with the c-axes vertically oriented. (A) Parallel free DNAs (yellow) in 2-()() are at 5 ? shortest distance to each other and stack at both ends with 192-bound DNAs (magenta). Free DNA is not involved in any further crystal contacts. (B) Discontinous electron density for free DNA in 2-()(). Modeled complexes of 2-() and modeled bases of free DNA (see Figure 1B) in yellow. The central parts of free DNA marked by dotted lines (4 to 15 bp, Figure 1B) could not be modeled due to patchy electron density. The shortest distance between free DNAs is 8 ? and they are at an angle of 40°.
On the ‘left’ sides of free DNA (Figure 1), nucleotides C17 of free and 192-bound DNA are not in helical arrangement and not seen in the electron density as they are disordered but were confirmed by MALDI-TOF-spectrometry of dissolved crystals (data not shown). Both, 2-bound and free DNAs, stack with bp G16–C-1' that are related by pseudo-2-fold symmetry.
On the ‘right’ sides of free DNA (Figure 1), the two 3'-G16' overhangs lie in the minor groove of the adjacent duplex and interact with both 5'-G0 to form two consecutive G16'*(G0–C15') base-triplets (Supplementary Figure 11) with similar geometry as reported (20–22).
In the crystals of 192 bound to () in space group C2 and bound to () in space group P21 (Table 2), both 2–DNA complexes interact by protein–protein contacts to form layers parallel to the crystallographic a, b planes. The crystallographic a, b-axes in space group P21 correspond to b, a-axes in C2 thus reflecting space group and lattice packing similarities. The pseudo-continuous DNA helices (Figures 1 and 3) are oriented in c-direction. The doubled c-axis in C2 relative to that in P21 is due to the C-centering, and the DNA helices are parallel to the c-axis at shortest inter-helix distance of 5 ? in 2-()() (base-triplets indicated by ), possibly stabilized by bridging water. By contrast, in 2-()() the pseudo-continuous DNA helices are parallel to the crystallographic a, c plane but inclined at an angle of 40° towards the c-axis and at least 8 ? apart, and poor electron density (Figure 3B) indicates that they are partially disordered as shown by B factors >100 ?2. Consequently, bp 4 to 15 of free ()-DNA in 2-()() could not be modeled (Figures 1B and 3B).
In both crystal unit cells (space groups C2 and P21, Table 2), the DNA-bound 192 show minor structural changes compared to the X-ray structure of free 202 (11). This concerns the loop connecting -helices 1 and 2 and to a minor extent the ?-sheet (residues 27–32, 0.8 ? r.m.s. deviation).
Since gel permeation chromatography of 2-() at 20 mM Tris–HCl (pH 7.5), 150 mM NaCl indicated an apparent molecular mass of 40 kDa (calculated Mw = 33.4 kDa), we assume that under crystallization conditions with 2.4 M Na2-malonate, the complexes associate pairwise through base-triplets G16'*(G0–C15') to form the complex 2-()()-2. If one 2-() of this complex forms layers mediated by protein–protein contacts as in the present structures, it appears that the other 2-() has to release 2 for packing reasons, thereby giving rise to the crystallized 2-()(). The other crystal contact with only DNA–DNA stacking interactions is more flexible than the base-triplets and permits formation of a regular crystal lattice. The same crystallization scenario applies for 2-().
Protein–DNA interactions
Formation of G16'*(G0–C15') base-triplets induces distortions at the 5' ends of heptads A1–A7 as indicated by different, partly water-mediated interactions to 192 compared to heptads A8–A14 in () and A8–T14 in (), Figures 4 and 5. For this reason, we focus here on protein–DNA interactions for the less distorted dimer 19B/B' bound to heptad A8–A14 of 2-() (Figures 2A and 4(left)). In both structures, for each 192-bound heptad the direct (not water-mediated) protein–DNA contacts are comparable.
Figure 4 Schematic representation of interactions between 192 and DNA in 2-() (left) and 2-() (right). The orientation of heptads is indicated by blue arrows and palindromic symmetry by the black ellipse. Yellow bases specifically interact with 192. Residues labeled green for 19 A/B and black for 19 A'/B'. Hydrogen bonds to phosphate oxygens (red) in thin black lines, specific hydrogen bonds to bases in the major groove in black arrows, blue lines indicate water (blue dots) mediated hydrogen bonds. K28 (K28') and Arg31 (Arg31') interact differently in subunits A/B and A'/B'.
Figure 5 Stereo view of specific interactions between 19 B/B' and bound heptad A8–A14 in 2-(). Water molecules; red spheres, hydrogen bonds; dashed lines, phosphorus atoms; magenta, directly contacted bp C11–G4' and base G2' in yellow.
In the major grooves, base pair specific interactions are formed with Thr29 and Arg31 located on the ?-sheet. Thr29O and Thr29'O of 19B and B' bind specifically to the central bp G4'–C11, and Arg31N, hydrogen bond with base G2' (Figure 4(left)). In contrast, the corresponding Arg31' of 19B' hydrogen bonds with N' to Thr29O of 19B and with N1',2' through three water molecules to bases G4', A5' and A9, Figures 4 and 5.
To see whether 192 would bind symmetrically to a palindromic heptad featuring two G to provide both, Arg31 and Arg31', with potential binding partners, we determined the 2.9 ? resolution crystal structure of 2 in complex with a mutated () where bp A9–T6' was replaced by G9–C6' (heptad sequence: 8AGTCACA14, Table 2). However, the same interaction pattern was found as in the original heptad. This agrees with similar (2-fold weaker) binding of wt 2 to mutated () with the same replacement in the first heptad (1AGTCACA7) compared to the original operator (8). A referee suggested to test the binding of 192 to heptads with pseudo-palindromic symmetry, 5'-TGTCACA-3'. In view of the binding geometry of 192 to 8AGTCACA14, we question whether this would provide novel knowledge. This is because the 5'-AAT- or 5'-AGT-ends of the heptads contact 192 in all cases exclusively through unspecific interactions with phosphate groups or are mediated by water molecules (Figure 4A and B). Hence, base pair exchanges in this part of the heptads should not significantly affect binding of 192.
Backbone phosphates of all four heptads contact helices 1 and 2 of 19 subunits with pseudo-2-fold symmetry (Figure 4). In 19B, the 5'-phosphate of A9 caps the N-terminus of 2 by hydrogen bonding to peptide amides of V51 and K52 in a pattern known for RHH2 proteins and other repressors (4,23), and His37N, K41N of 1 bind to the 5'-phosphate of T10. Corresponding residues of 19B' interact with 5'-phosphates of G2' and T3'. However, K28N located on the ?-sheet of 19B forms a salt-bridge with the 5'-phosphate of G4' whereas K28'N on 19B' is >5 ? away from the corresponding 5'-phosphate of C11.
The asymmetry in the binding of each 192 to its particular heptad is reflected by superimposition of C atoms of monomers (A on A' and B on B', Figure 2A) in each 192 showing 0.6 ? r.m.s. deviation partly associated with structural differences in the loops connecting helices 1 and 2 (residues 46–48) and in the ?-strands. When 192 dimers are superimposed on each other in the same orientation (A on B and A' on B'), i.e. the bound heptads point in the same direction, r.m.s. deviation of only 0.3 ? confirms that all 192-heptad interactions are similar.
Thr29 is essential for specific operator binding
To test the importance of Thr29 for heptad sequence recognition in vivo studies were conducted showing that 2Thr29Ala failed to completely repress promoter P utilization (Table 1). Binding of 2Thr29Ala to PcopS operator DNA embedded in 300 bp DNA was tested by electrophoretic mobility shift assays (EMSA) (Supplementary Figure 10). Protein–DNA complexes formed with 192 but Thr29Ala required 100-fold higher concentration (compared to 192) that yielded prominent but unspecific binding to PcopS as confirmed by DNase I footprinting (data not shown).
192···192 interactions
Protein–DNA interactions in both complexes bury 1610 ?2 of solvent accessible surface area. Another 550 ?2 are buried by interaction of pseudo-2-fold axis related 1 helices of adjacent monomers A' and B in 2-() and A' and B' in 2-(), respectively (Figure 1 and Figure 2A and B). In detail, bifurcated hydrogen bonds between His38N and Ala45'O/Lys46'O are augmented by hydrophobic contacts between Ile42, Ala45 of both subunits (Figure 6). Positions of His38 and His38' are identical in both complexes whereas the hydrophobic side-chain of Ile42 adopts different rotamers without affecting the size of the buried interface. These interactions ensure cooperative binding when several 2 associate with multiple heptad repeats as found in natural operators.
Figure 6 Interactions at the 192···192 interface (A) between helices 1' and 1 of subunits A' and B of 2-(), (see Figure 2A). The 192 are related by a pseudo-2-fold axis perpendicular to the paper plane (red ellipse). Residues with contacts <4 ? to neighbor helix are depicted with side-chains, colors as in Figure 2. Bifurcated hydrogen bonds between His38 (His38') and Ala45O/Lys46O to Ala45'O/Lys46'O of the other subunit in dashed lines. (B) Superimposition as in Figure 2B with view on the dimer–dimer interface to show that dimer–dimer interactions are comparable in both structures.
To test the importance of His38, plasmid-borne His38Val was constructed to remove the bifurcated hydrogen bonds but to maintain the hydrophobic character of the interface. 2His38Val repressed P utilization in vivo with 7-fold lower efficiency than wt 2 (Table 1), indicating the important role of His38 for cooperative binding between 2 and multiple consecutive heptads.
Conformation of free and 192-bound DNA
The DNAs in 2-() and 2-() are nearly straight B-form with average helical twist of 36° (range 23° to 43°; Supplementary Figure 12), but show distinct features. The major groove width in 2-() shows strong modulation depending on nucleotide sequence, being 13 to 14 ? except for the 7AAAT10 tract (11 to 12 ?) compared to 11.7 ? for ideal B-DNA (24). In contrast, the major groove width in 2-() is more continuous (13 to 14 ?) due to the palindromic symmetry (Figure 7A). In the 7AAAT10 tract of 2-() and in the 7AAGT10 tract of 2-() the minor groove is narrowed to 2.7 and 4 ?, respectively, compared to 5.7 ? in ideal B-DNA (24) due to an average negative base pair inclination of –8° (Figure 7C), as frequently found for A-tracts (25). This is associated in 7AAAT10 of 2-() with binding of a single spine of four water molecules, one in each base step as reported in (26). Strong opposite buckles of the C–G bp in the 4CAC6 (–15°, 15°) and 11CAC13 (–10°, 15°) segments of 2-() (Figure 7B) widen the major groove to 13 ? to accommodate the ?-sheet, and adjacent (up- and downstream) minor grooves are narrowed (Figure 7A). Similarly, in 2-() the C–G bp in segment 4CAC6 show comparable buckles of –10°, 19° but in 9GTG11 the buckles are reduced to –8°, 5°. Of the four heptads, three feature two subsequent CA steps that are known for their ability to bend B-DNA through positive roll (20,27,28). However, the CA steps do not exhibit unusual structure except for a moderate slide movement of 1 ? at the central C–G bp where 192 binds to both bases (Supplementary Figure 12).
Figure 7 Comparison of operator DNA free and in complex with 192. (A) Minor (m) and major (M) groove widths along the central 11 bp of each operator, as measured by shortest P-P distances (less two phosphate group radii, 5.8 ?) indicated by horizontal dotted black lines for ideal B-DNA (24), 5.7 ? and 11.7 ?, respectively, Note minor groove narrowing in the 7AAAT10 tract of in 192-bound and free () DNA and reduced narrowing in 7AAGT10 tract of (192)2-(). (B) Base pair buckle (°) and (C) inclination angles (°) of each operator.
In the crystal lattice, free ()-DNA contacts 192-bound ()-DNA with both ends to form pseudo-continuous helices (Figure 1) and is not distorted by any further crystal contacts (Figure 3A). Interestingly, free ()-DNA shows comparable structural features as found for 192-bound ()-DNA (Figure 2C). In the 7AAAT10 tract of free ()-DNA, base pair inclination angles are negative, and the minor groove is narrowed to 3.5 ? whereas the major groove is overall widened to an average of 12.5 ?. Additionally, strong opposite buckles are observed for C–G base pair in both CAC segments that are comparable to 192-bound ()-DNA (Figure 7B). These observations suggest that the described significant deviations from ideal B-DNA are a consequence of nucleotide sequence and not induced by binding to 192 (Figures 2C and 7).
However, conformational changes upon repressor binding concern base pair parameters helical twist, roll and slide (Supplementary Figure 12). In free ()-DNA these parameters cluster around values assigned to ideal B-DNA, but show a negative roll of –6° and positive slide of 1 ? for the central G11'–C4 and G4'–C11 bp of both heptads, and helical twist angles for bp A5–T10' and A12–T3' decrease to 25° in 192-bound ()-DNA. Since these deviations are similar for both heptads of 2-() they can be attributed to 192 binding. Such comparison is not possible with free ()-DNA as it could not be modeled completely.
Comparison of 2-() and 2-()
The pseudo-2-fold axis relating the monomers in 192 passes with 0.3 ? downstream (5'3') offset through the central G–C bp of each heptad. Consequently, these symmetry axes are separated by 7 bp in 2-() but they are 0.6 ? closer in 2-(), see Figure 2B. Despite the different dimer–dimer separations helices 1 and 1' forming the 192...192 interfaces superimpose well in both complexes, consistent with similar dissociation constants (kD 20 nM). In contrast, the separation between the two 192 will be 0.6 ? wider in diverging heptads 2-(). Assuming that the interaction pattern between 192 and heptads with () orientation is similar as with () and () heptad orientations, the expected 0.6 ? longer 1···1 contacts are probably less favorable and diminish cooperativity. This agrees with the 6-fold weaker affinity of wt 2 to heptads in () arrangement, the finding that 2-() dissociated during gel filtration, and the drastically reduced binding affinity to heptad repeats spaced by one or more additional base pair (8).
Whereas DNA in 2-() is nearly straight, the superimposition in Figure 2B reveals a kink of 12° at the G11–C4' bp at the centre of the heptad bound to dimer B/B' of 2-(). This kink is associated with 12° rotation of dimer B/B' of 2-(), the rotational pivot point being located in the 192···192 interface. By virtue of these two motions, the position and orientation of helix 1' of subunit B' remain almost as in 2-() and similar 192···192 interactions explain comparable dissociation constants of both complexes.
Structural model of 2-bound to natural operators
Extrapolation of the structures of 2-() and 2-() allowed modeling of 192 in complex with natural P promoter (Figure 8). The model implies that wt 2 binds as left-handed matrix to right-handed, straight B-type operator DNA, each 192 being displaced relative to its neighbor by 7 bp and rotated by 252°. Figure 8B shows that the negatively charged sugar-phosphate backbone of DNA faces positively charged surface of 192.
Figure 8 Model of seven 192-bound heptads () of the natural promoter P. (A) DNA in space filling and 192 as orange/red ribbons, (B) DNA in stick and 192 in surface representation colored according to electrostatic potential (negative and positive charges red and blue, respectively). The model is based on the structures of 2-() and 2-(). Repressors form a left-handed protein-matrix winding around the nearly straight operator.
DISCUSSION
Implications of 192-DNA structures for regulation of transcription
The inc18 family plasmids harbor genes to control their copy number, accurate segregation and stable maintenance during cell division. Since expression of these genes is regulated by the common 2 repressor, a unique mechanism has evolved to fine-tune repressor affinity for the different operators. How this is achieved is shown by the present study. It clearly indicates that the pseudo-symmetric 2 binds with 0.3 ? downstream (5'3') offset relative to the center G–C base pair of the cognate heptad. Since the operators are nearly straight B-DNA, different heptad numbers and orientations lead to different distances between 1 helices of adjacent 2, thereby modulating cooperative interactions between 2 and different operators. The ability to bind to palindromic as well as to non-palindromic operators is a unique feature of 2 and is not shared by other member of the RHH2 family. We associate this with the interactions between 192 that are related by a pseudo-2-fold rotation axis (Figure 2A and Figure 6A and B) so that they interact comparably, no matter what the orientations of adjacent heptads are, (), () or (). It is unlikely that the deleted N-termini would contribute to 192 ···192 interactions as the N-termini of 192 point away from the 192 ···192 interface (Figures 2A and 6A).
When RHH2 bind to DNA, we have to consider two different binding characteristics, the ‘variable’ binding of flexible side-chains of ?-sheets and 1 helices to major groove and phosphate groups, respectively, and ‘stiff’ capping of N-termini of 2 helices by rigid main-chain NH hydrogen bonding to phosphates. Although Arc2, CopG2 and MetJ2 bend minimal cognate operators with two repressor binding sites by 50 to 60°, their cores do not rearrange significantly compared to free repressors and are structurally similar . In contrast, Arc2, CopG2 and MetJ2 superimpose on DNA-bound 192 with higher r.m.s. deviations of 2.3 to 2.8 ? because the ?-sheet of 192 protrudes 3 ? less from the repressor surface and is closer to the N-termini of helices 2 (Supplementary Figure 13). Due to this geometry of 192, the N-termini of helices 2 are in the correct position to clamp the phosphate backbones of straight operator DNA when the ?-sheet is inserted into the major groove (Figures 2 and 4). In contrast, the other three repressors have to bend DNA around the more protruding ?-sheet to place the phosphate backbone in hydrogen bonding distance to N-H groups of the N-termini of their helices 2. It is notable that the distances between phosphates bound by 2 helices is 5 bp in 2 (Figure 4) but 6 bp in the other three RHH2 since DNA has to follow a longer path when bending around the protruding ?-sheet.
The N-termini of RHH2 proteins are associated with different functions
The N-terminus in CopG2 has no obvious function, in MetJ2 it is involved in binding the corepressor S-adenosylmethionine whereas in free Arc2 it is disordered but forms a 310 helix upon and contributes to DNA-binding. The N-terminal residues in 2 do not contribute to DNA-binding in vitro and in vivo as shown in the present work. However, we have recently shown that protein , a homolog to ParA proteins and involved in active plasmid partitioning during cell division, is activated by wt 2 but not by 192 (A. Cicek, F. Pratto, W. Weihofen, J.C. Alonso and W. Saenger, unpublished data). This suggests that 2 is the yet missing ParB protein of inc18 family plasmids in the known ParA/ParB family of plasmid partitioning systems (29).
Direct and indirect readout of DNA sequence
Cocrystallized 192-bound and free ()-DNA show similar phosphate backbone conformation (Figure 2C) with significant deviations from ideal B-DNA (Figure 7), indicating that conformation is predominantly dependent on DNA sequence and not induced by 192 binding. This provides a good example for the ‘indirect’ readout of local DNA conformation by 192 that depends on the particular nucleotide sequence (30). The combination with ‘direct’ readout resulting from interactions of 2 repressor amino acids with heptad bases increases protein–DNA-binding specificity and affinity. This view is consistent with studies in which the heptad nucleotide sequence was mutated base pair by base pair (8). Mutations at heptad positions 2, 3 and 5 (5'-AATCACA-3') resulted in at least 4-fold weaker binding to wt 2, although these base pair are not directly contacted by repressor side-chains contrasting direct read out positions 4 and 6 that were much more sensitive against mutations.
Protein Data Bank accession codes
Atomic coordinates and structure factors have been deposited with the Protein Data Bank under accession codes 2bnw, 2bnz and 2cax.
SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.
ACKNOWLEDGEMENTS
The authors acknowledge beamtime and technical support at the Protein Structure Factory beamlines of Free University Berlin at BESSY/Berlin. This work was partially supported by Fonds der Chemischen Industrie to W.S., BMC2003-00150 from Dirección General de Investigación, Ministerio de Educación y Ciencia to J.C.A., and QLK3-CT-2001-00277 from the European Community to J.C.A. and W.S. Funding to pay the Open Access publication charges for this article was provided by Free University Berlin.
REFERENCES
de la Hoz, A.B., Ayora, S., Sitkiewicz, I., Fernandez, S., Pankiewicz, R., Alonso, J.C., Ceglowski, P. (2000) Plasmid copy-number control and better-than-random segregation genes of pSM19035 share a common regulator Proc. Natl Acad. Sci. USA, 97, 728–733 .
Raumann, B.E., Rould, M.A., Pabo, C.O., Sauer, R.T. (1994) DNA recognition by ?-sheets in the Arc repressor-operator crystal structure Nature, 367, 754–757 .
Gomis-Ruth, F.X., Sola, M., Acebo, P., Parraga, A., Guasch, A., Eritja, R., Gonzalez, A., Espinosa, M., del Solar, G., Coll, M. (1998) The structure of plasmid-encoded transcriptional repressor CopG unliganded and bound to its operator EMBO J, . 17, 7404–7415 .
Somers, W.S. and Phillips, S.E. (1992) Crystal structure of the met repressor-operator complex at 2.8 A resolution reveals DNA recognition by ?-strands Nature, 359, 387–393 .
Brantl, S. (1994) The copR gene product of plasmid pIP501 acts as a transcriptional repressor at the essential repR promoter Mol. Microbiol, . 14, 473–483 .
Pujol, C., Chedin, F., Ehrlich, S.D., Janniere, L. (1998) Inhibition of a naturally occurring rolling-circle replicon in derivatives of the -replicating plasmid pIP501 Mol. Microbiol, . 29, 709–718 .
Ceglowski, P., Boitsov, A., Chai, S., Alonso, J.C. (1993) Analysis of the stabilization system of pSM19035-derived plasmid pBT233 in Bacillus subtilis Gene, 136, 1–12 .
de la Hoz, A.B., Pratto, F., Misselwitz, R., Speck, C., Weihofen, W., Welfle, K., Saenger, W., Welfle, H., Alonso, J.C. (2004) Recognition of DNA by protein from the broad-host range Streptococcus pyogenes plasmid pSM19035: analysis of binding to operator DNA with one to four heptad repeats Nucleic Acids Res, . 32, 3136–3147 .
Glass, C.K. (1994) Differential recognition of target genes by nuclear receptor monomers, dimers, and heterodimers, Endocr. Rev, . 15, 391–407 .
Misselwitz, R., de la Hoz, A.B., Ayora, S., Welfle, K., Behlke, J., Murayama, K., Saenger, W., Alonso, J.C., Welfle, H. (2001) Stability and DNA-binding properties of the regulator protein from the broad-host range Streptococcus pyogenes plasmid pSM19035 FEBS Lett, . 505, 436–440 .
Murayama, K., Orth, P., de la Hoz, A.B., Alonso, J.C., Saenger, W. (2001) Crystal structure of transcriptional repressor encoded by Streptococcus pyogenes plasmid pSM19035 at 1.5 A resolution J. Mol. Biol, . 314, 789–796 .
Otwinowski, Z. and Minor, W. (1997) Processing of X-ray diffraction data collected in oscillation mode In Carter, C.W., Jr and Sweet, R.M. (Eds.). Macromolecular Crystallography Part A, Academic Press pp. 307–326 .
Storoni, L.C., McCoy, A.J., Read, R.J. (2004) Likelihood-enhanced fast rotation functions Acta Crystallogr. D Biol. Crystallogr, . 60, 432–438 .
Murshudov, G.N., Vagin, A.A., Dodson, E.J. (1997) Refinement of macromolecular structures by the maximum- likelihood method Acta. Crystallogr. D Biol. Crystallogr, . 53, 240–255 .
Vagin, A. and Teplyakov, A. (1997) MOLREP: an automated program for molecular replacement J. Appl. Cryst, . 30, 1022–1025 .
Laskowski, R.A., MacArthur, M.W., Moss, D.S., Thornton, J.M. (1993) PROCHECK: a program to check the stereochemical quality of protein structures J. Appl. Cryst, . 26, 283–291 .
Kraulis, P.J. (1991) MOLSCRIPT: a program to produce both detailed and schematic plots of protein structures J. Appl. Cryst, . 24, 946–950 .
Merritt, E.A. and Bacon, D.J. (1997) Raster3D: photorealistic molecular graphics Meth. Enzymol, . 277, 505–524 .
Lavery, R. and Sklenar, H. (1988) The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids J. Biomol. Struct. Dyn, . 6, 63–91 .
Schultz, S.C., Shields, G.C., Steitz, T.A. (1991) Crystal structure of a CAP–DNA complex: the DNA is bent by 90 degrees Science, 253, 1001–1007 .
Spink, N., Nunn, C.M., Vojtechovsky, J., Berman, H.M., Neidle, S. (1995) Crystal structure of a DNA decamer showing a novel pseudo four-way helix–helix junction Proc. Natl Acad. Sci. USA, 92, 10767–10771 .
Nunn, C.M., Garman, E., Neidle, S. (1997) Crystal structure of the DNA decamer d(CGCAATTGCG) complexed with the minor groove binding drug netropsin Biochemistry, 36, 4792–4799 .
Aggarwal, A.K., Rodgers, D.W., Drottar, M., Ptashne, M., Harrison, S.C. (1988) Recognition of a DNA operator by the repressor of phage 434: a view at high resolution Science, 242, 899–907 .
Saenger, W. Principles of Nucleic Acid Structure, (1984) NY Springer-Verlag .
Shatzky-Schwartz, M., Arbuckle, N.D., Eisenstein, M., Rabinovich, D., Bareket-Samish, A., Haran, T.E., Luisi, B.F., Shakked, Z. (1997) X-ray and solution studies of DNA oligomers and implications for the structural basis of A-tract-dependent curvature J Mol. Biol, . 267, 595–623 .
Han, G.W., Kopka, M.L., Cascio, D., Grzeskowiak, K., Dickerson, R.E. (1997) Structure of a DNA analog of the primer for HIV-1 RT second strand synthesis J. Mol. Biol, . 269, 811–826 .
Chen, S., Gunasekera, A., Zhang, X., Kunkel, T.A., Ebright, R.H., Berman, H.M. (2001) Indirect readout of DNA sequence at the primary-kink site in the CAP–DNA complex: alteration of DNA binding specificity through alteration of DNA kinking J. Mol. Biol, . 314, 75–82 .
El Hassan, M.A. and Calladine, C.R. (1998) Two distinct modes of protein-induced bending in DNA J.Mol. Biol, . 282, 331–343 .
Pogliano, J. (2002) Dynamic cellular location of bacterial plasmids Curr. Opin. Microbiol, . 5, 586–590 .
Otwinowski, Z., Schevitz, R.W., Zhang, R.G., Lawson, C.L., Joachimiak, A., Marmorstein, R.Q., Luisi, B.F., Sigler, P.B. (1988) Crystal structure of trp repressor/operator complex at atomic resolution Nature, 335, 321–329 .(Wilhelm Andreas Weihofen, Aslan Cicek, F)