Cre Induces an Asymmetric DNA Bend in Its Target loxP Site*,

Linda Lee {ddagger}, Linda C. H. Chu and Paul D. Sadowski §

From the Department of Medical Genetics and Microbiology, University of Toronto, Toronto M5S 1A8, Canada

Received for publication, March 5, 2003 , and in revised form, April 8, 2003.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Cre initiates recombination by preferentially exchanging the bottom strands of the loxP site to form a Holliday intermediate, which is then resolved on the top strands. We previously found that the scissile AT and GC base pairs immediately 5' to the scissile phosphodiester bonds are critical in determining this order of strand exchange. We report here that the scissile base pairs also influence the Cre-induced DNA bends, the position of which correlates with the initial site of strand exchange. The binding of one Cre molecule to a loxP site induces a ~35° asymmetric bend adjacent to the scissile GC base pair. The binding of two Cre molecules to a loxP site induces a ~55° asymmetric bend near the center of the spacer region with a slight bias toward the scissile A. Lys-86, which contacts the scissile nucleotides, is important for establishing the bend near the scissile GC base pair when one Cre molecule is bound but has little role in positioning the bend when two Cre molecules are bound to a loxP site. We present a model relating the position of the Cre-induced bends to the order of strand exchange in the Cre-catalyzed recombination reaction.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
The Cre recombinase of bacteriophage P1 assists in the efficient segregation of the low copy P1 plasmid by resolving dimeric lysogenic P1 plasmids into monomeric units (1). Cre is a member of the {lambda} integrase or tyrosine recombinase family of conservative site-specific recombinases (25). The tyrosine recombinases share a common mechanism of catalysis. A conserved tyrosine (Tyr-324 in Cre) serves as the catalytic nucleophile that cleaves a specific phosphodiester bond in the DNA target sequence and attaches the recombinase to the DNA via a 3'-phosphotyrosine bond (see Fig. 1a below) (39). Recombination proceeds via two sequential strand exchanges, forming a four-armed Holliday structure as an intermediate (1017).



View larger version (28K):
[in this window]
[in a new window]
 
FIG. 1.
Cre-mediated site-specific recombination. a, schematic diagram of the mechanism of Cre-loxP recombination. Two Cre molecules bind to and induce an asymmetric bend in each of the loxP sites (25). We refer to the Cre subunits and the DNA strands using the nomenclature of Van Duyne and coworkers (24, 25). One of the two Cre molecules bound to the lox site is poised to cleave the adjacent scissile phosphate ("cleaving" subunit, dark oval), whereas the other "non-cleaving" monomer (gray) is in the inactive conformation (21, 23, 25). The "crossing" strand (thin lines in i–iii) is defined as the strand containing the activated scissile phosphate and is pointed toward the center cavity of the synapse ready for strand exchange after the first strand cleavage (25). The "continuous" strand (thick lines in i–iii) contains the inactive scissile phosphate and is on the concave side of the DNA bend (25). Cleavage initiates on the crossing strands (bottom strand of loxP, thin lines) and results in covalent attachment of the cleaving Cre molecule to the 3'-end of the DNA via a phosphotyrosine bond. Exchange of the cleaved strands leads to the formation of a Holliday junction intermediate. The Holliday intermediate is then resolved on the top strand (thick lines) to yield reciprocal recombinant products. b, the wild type loxP site. The loxP site is composed of two 13-bp inverted repeats or symmetry elements a and b (normal type and represented by the outer arrows) flanking an asymmetric 8-bp spacer region (boxed boldface type). The orientation of the loxP site is designated by the middle arrow. We define the top and bottom strands of the lox site as those illustrated here. The scissile A and G nucleotides are defined as the nucleotide immediately 5' to the top and bottom cleavage sites (small vertical arrows) at positions 4' and 4, respectively. The numbers adjacent to the cleavage sites indicate the order of strand cleavage and exchange. The cI and cII bends result from the binding of one (cI) and two (cII) Cre molecules to a loxP site, respectively (see Fig. S2b). The location of the cI and cII bend centers in loxP induced by wild type (Wt) and K86A Cre proteins are indicated by vertical arrows above and below the sequence, respectively (Table II). The dashed arrows denote the cI bend center within a DNA containing only one symmetry element (se a or se b; Table I). The error in estimating the position of the bend center was about ±4 bp (see Supplementary Data, Fig. S7). The direction of the Cre-induced bend (either toward the major or minor groove of the DNA) at the bend center is given in parentheses. The direction of the cI bend at position 8' in loxP and se a is ambiguous, but we speculate that it is probably toward the minor groove like the other cI bends (see "Results"). c, the lox4 site. In the mutated lox4 site, the scissile base pairs at positions 4' and 4 have been interchanged relative to the wild type loxP site. The schematic of the lox4 site and the Cre-induced bends within the lox4 site are illustrated in a similar manner to loxP in b.

 


View this table:
[in this window]
[in a new window]
 
TABLE II
Summary of the Cre-induced DNA bends in various lox sites

 

View this table:
[in this window]
[in a new window]
 
TABLE I
The sequence of the lox sites studied in this report

 
The DNA target sequence for the Cre protein is called loxP (Fig. 1b) and consists of two identical 13-bp inverted symmetry elements surrounding an 8-bp asymmetric sequence (18). This 8-bp sequence defines the orientation of the loxP site, and we refer to it as the "spacer" region. Each symmetry element serves as a site for sequence-specific binding of a Cre monomer (1921). The scissile phosphodiester bonds are 6 bp apart, and the Cre protein attaches covalently to the 3'-phosphoryl A residue on the top strand and the 3'-phosphoryl G residue on the bottom strand (22). We refer to these nucleotides as the "scissile A" and "scissile G," respectively.

Several crystal structures of the Cre-lox complexes have provided remarkable insights into the conformations of the various intermediates in the Cre-lox reaction (21, 2327). Two Cre-bound loxP sites are brought together by a cyclic network of protein-protein interactions to form an approximately square-planar synaptic complex (Fig. 1a) (21, 2330). The two Cre molecules bound to a lox site are conformationally and functionally different: one is poised to cleave the DNA ("cleaving" subunit), and the other is in an inactive conformation ("non-cleaving" subunit; see Fig. 1a) (21, 2327).

Cre catalyzes recombination with a defined order of strand exchange (see Fig. 1, a and b): it first cleaves and exchanges the two bottom strands adjacent to the scissile G nucleotide to form the Holliday intermediate (11, 2932) that is then resolved preferentially on the top strand adjacent to the scissile A nucleotide to generate reciprocal recombinant products (11, 3133). We found previously that the order of strand exchange was dictated primarily by the scissile base pairs at positions 4' and 4 in the loxP site (32, 33). The order of strand exchange was reversed when the scissile base pairs were interchanged in the mutated lox4 site (see Fig. 1c). Furthermore, we found that Lys-86, which contacts the scissile nucleotides (21, 2327), is important for establishing the strand selectivity in the resolution of the loxP-Holliday intermediate but not in the initiation of recombination between the loxP sites (32, 33).

The crystal structure of the Cre-lox synaptic complex reveals the presence of an asymmetric DNA bend in the lox spacer region, positioned 5 bp away from the activated cleavage site (25). Guo et al. (25) proposed that the position and/or direction of the DNA bend dictate the site of initial strand exchange. In this report we have further characterized the roles of the scissile base pairs and the Lys-86 residue in the Cre-induced DNA bending to better understand the basis for the order of strand exchange. We find that Cre induces asymmetric bends in the loxP site, and the position of the bends is dictated by the scissile base pairs: the binding of one Cre molecule (complex I, cI)1 induces a bend near the margin of the spacer region adjacent to the scissile G nucleotide, whereas the binding of two Cre molecules (complex II, cII) induces a bend near the center of the spacer region with a slight bias toward the scissile A (see Fig. 1, b and c). Lys-86 has a role in positioning the Cre-induced bend in the loxP cI, but not in the loxP cII. Changes in the Cre-induced DNA bends within cII correlate with the site of initial strand exchange as originally proposed by Guo et al. (25). We present a model relating the position of the Cre-induced bends to the order of strand exchange in the Cre-catalyzed recombination reaction.


    EXPERIMENTAL PROCEDURES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Oligonucleotides, DNA Constructs, and Proteins—The S82 substrates used in the chemical footprinting experiments were constructed by annealing two 82-nucleotide complementary oligonucleotides, 5'-GATCCAGACTGCAGCG(lox)GCTAATCGATCAGAATCTGGCTATGACGATCC-3'/3'-CTAGGTCTGACGTCGC(lox)CGATTAGCTAGTCTTAGACCGATACTGCTAGG-5', where (lox) is the lox sequence of interest (see Table I).

The circularly permuted plasmids were constructed by ligating the 40 bp of annealed, complementary oligonucleotides 5'-tcgac(lox)g-3'/3'-g(lox)cagct-5', containing the lox sequence of interest (Table I) flanked by SalI ends, into SalI-digested pBend2 (34). This results in the insertion of the lox site between two tandemly repeated sets of restriction enzyme sites. The circularly permuted DNA substrates were obtained by digestion of the plasmid with EcoRV, SmaI, NruI, SspI, BamHI, MluI, BglII, NheI, and SpeI, and isolation of the 162-bp digested fragments.

The phasing plasmids were constructed by excising the 156-bp BamHI-BamHI DNA fragment from the plasmids pK10, pK12, pK14, pK16, pK18, and pK20 (35) and ligating in the 90-bp annealed, complementary oligonucleotides: 5'-GATCCACGATCAGACTGCAGCCATGGCACG(lox)GCTAAGATCTCAGAATCTGGCTATCG-3'/3'-GTGCTAGTCTGACGTCGGTACCGTGC(lox)CGATTCTAGAGTCTTAGACCGATAGCCTAG-5' containing the lox sequence of interest (Table I). Constructs containing the lox site inserted in the forward and reverse orientations were isolated. Phasing substrates were obtained by isolating the 347- to 357-bp RsaI-PvuII fragments. The loxSA site introduced a novel RsaI site in the symmetric spacer region and so the loxSA phasing substrates were obtained by PCR using the primers RsaI2959F (5'-ACATATTGTCGTTAGAACGCG-3') and PvuII270R (5'-CTGGCTTATCGAAATTAATAC-3'). ~0.2 µg of the loxSA phasing plasmid was amplified in 50 µl of 3 mM MgCl2, 0.4 mM of each dNTPs, 1 mM of each primers, 2.5 units of Taq polymerase (Invitrogen) for 35 PCR cycles of heating at 95 °C for 45 s, 55 °C for 1 min, and 72 °C for 1 min. The 347- to 357-bp PCR products were purified on a preparative 4% native PAGE. Note that the S182 substrates used in the recombination assays previously described (32) were 182-bp SacI-PvuII fragments isolated from the phasing plasmids pK10LPA and pK10L4A, which are pK10-derived plasmids with the loxP and lox4 sites, respectively, inserted in the forward orientation.

All oligonucleotides were synthesized by Invitrogen. The procedures for the purification of the oligonucleotides, annealing, and 5'-32P labeling the DNA substrates were as described previously (32, 33). The restriction enzymes and T4 polynucleotide kinase were from New England BioLabs. The Cre proteins were purified from an induced Escherichia coli BL21(DE3 pLysS) culture as described previously (33, 36).

DMS Methylation Protection Analysis—2 nM S82 substrate (5'-32P-labeled on either the top or the bottom strand) was incubated with 0.5 µM Cre in 50 µl of 50 mM sodium cacodylate, pH 8, 30 mM NaCl, and 0.05 mg/ml denatured calf thymus DNA at room temperature for 15 min. The reaction was then treated with 3.3 µl of 10% DMS (v/v, diluted in ethanol, Aldrich) at room temperature for 6 min (37). Analysis of a 5-µl aliquot of the DMS-treated reaction on a native 6% polyacrylamide gel showed that cII was the predominant species with some unbound substrate, cI, and higher order complexes present (data not shown) (32). The DMS reaction was quenched as described previously (38). The modified DNA was cleaved at methylated G > A residues in alkali as described by Craig and Nash (39) and concentrated by ethanol precipitation. The cleaved DNA was analyzed by 8% denaturing PAGE (Diamed SequaGel).

Dried gels were scanned using a Amersham Biosciences Phosphor-Imager and quantified using the ImageQuaNT 5.0 software program. The peak intensity of each band was expressed relative to the total intensity in the lane and then normalized to the relative intensity of the band corresponding to nucleotide "25 G" on the top strand or "19' G" on the bottom strand (Fig. 2, a and b). The extent of protection or enhancement was calculated as a ratio of the normalized relative band intensities in the presence or absence of Cre (+ to – Cre ratio).



View larger version (76K):
[in this window]
[in a new window]
 
FIG. 2.
DMS methylation protection of loxP and lox4 by wild type and K86A Cre proteins. A 82-bp DNA (S82) containing either the loxP site (lanes 1–4)orthe lox4 site (lanes 5–8) was 5'-radioactively labeled on the top strand (a) or the bottom strand (b).2nM DNA was incubated without (lanes 1 and 5) or with 0.5 µM of wild type (Wt), His-tagged wild type (HisWt) or HisK86A Cre protein, then treated with DMS as described under "Experimental Procedures." The DMS-treated DNA was cleaved at methylated G and A residues with alkali and analyzed on an 8% denaturing polyacrylamide gel. Representative autoradiograms are shown. The loxP sequence corresponding to the cleaved bands are indicated to the left of the audioradiogram. The residues mutated to create lox4 are indicated in parentheses. The symmetry elements (normal type) are represented by vertical arrows, and the spacer region (boldface type)is boxed with the small arrow denoting the cleavage site. The band intensities were quantified as described under "Experimental Procedures." The average ratios (from at least three experiments) of the normalized band intensities in the presence to absence of Cre (+ to – Cre ratio) are represented graphically on a log scale to the right of the autoradiogram. The nucleotides are numbered from the center of the lox site as shown in Fig. 1b, with the nucleotides in the spacer in boldface type and the flanking non-lox sequences in lowercase. The dotted vertical lines indicate 1.5-fold protection or enhancement (hypermethylation). Solid blue bars, loxP plus Wt Cre (non-tagged and His-tagged); striped blue bars, loxP plus K86A; solid red bars, lox4 plus Wt Cre; striped red bars, lox4 plus K86A. c, summary of the DMS footprints of the loxP and lox4 spacer region and proximal nucleotides by wild type (Wt) and K86A. Sites of protection (green diamonds) and enhancement (yellow triangles) greater than 1.5-fold are indicated. Although the scissile G at position 4' in lox4 was protected by wild type Cre by slightly less than 1.5-fold, the protection was reproducible.

 

OP-Cu Footprinting—2nM S82 substrate (5'-32P labeled on either the top or the bottom strand) was incubated with 1 µM Cre in 20 µl of Cre reaction buffer (32) for 30 min at room temperature. The OP-Cu footprinting (Supplementary Data, Fig. S1) was performed as described by Sigman et al. (40). The chemicals used in the OP-Cu reactions (1,10-phenanthroline and 3-mercaptopropionic acid) were from Sigma. The OP-Cu reaction was quenched with 1 µl of 28 mM neocuproine (Sigma) and 100 µl of 10 mM Tris-HCl (pH 8.0), 1 mM EDTA, 0.3 M sodium acetate, 20 µg/ml yeast tRNA, followed by phenol/chloroform extraction and ethanol precipitation. The DNA was analyzed by 8% denaturing PAGE (Diamed SequaGel).

Circular Permutation Analysis—The binding assay was performed by electrophoretic mobility shift assay (EMSA) as described previously (32) with the following modifications. 2 nM 5'-32P-labeled circular permutation substrates (obtained by digestion with various restriction enzymes as described above) was incubated with 0.05 µM Cre, and the reaction was analyzed on a 5% polyacrylamide gel. Determination of the bend centers and angles by circular permutation is described in detail in the Supplementary Data (Fig. S2).

Phasing Analysis—The binding assay was performed by EMSA as described previously (32) with the following modifications. 2 nM 5'-32P-labeled phasing substrates was incubated with 0.05 µM Cre, and the reaction was analyzed on a 4% polyacrylamide gel. Dried gels were scanned using an Amersham Biosciences PhosphorImager and analyzed using the ImageQuaNT 5.0 software program. The marker M bands in each lane were aligned. The phasing analysis was conducted according to the method of Zinkel and Crothers (35) as follows. The relative electrophoretic mobility (µ) of a protein-DNA complex was calculated as the mobility of the complex divided by the mobility of the unbound substrate to correct for small variations in the mobility of the unbound substrate. The relative mobility was then normalized to the average relative mobility (µave) of the particular complex from all six phasing substrates (µ/µave). The normalized relative mobilities were plotted as a function of the linker length. The linker length is defined as the distance (bp) between the center of the kinetoplast DNA that contains the A-tracts (35) and the middle of the lox site. The curve was plotted to the best fit polynomial curve (Microsoft Excel).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Guo and colleagues (25) have proposed that the position and/or direction of the asymmetric DNA bend in the lox spacer region revealed in the crystal structure of the Cre-lox synaptic complex determine the site of initial strand exchange. We therefore investigated whether the scissile base pairs and Lys-86 dictate the order of strand exchange by affecting the Cre-induced DNA bending. We have probed the DNA conformation of the loxP and lox4 sites (Fig. 1) bound to Cre using chemical footprinting, circular permutation, and phase-sensitive analyses.

Cre Increases the Sensitivity to DMS Methylation in the loxP and lox4 Spacer Region on the Strand Containing the Scissile G—DNA footprinting provides information about protein-DNA contacts and distortions in the DNA conformation such as bending and unwinding (4044). To detect subtle changes in the DNA conformation, we probed the Cre-bound DNA using small chemical compounds. We first analyzed the sensitivity of guanine and adenine residues in the loxP and lox4 sites to methylation by dimethyl sulfate (DMS). DMS methylates the N7 group of guanine in the major groove of DNA and, to a lesser extent, the N3 group of adenine in the minor groove (41).

An 82-bp DNA substrate (S82) containing either the wild type loxP or the mutated lox4 site was treated with DMS in the absence and presence of the wild type Cre protein (Fig. 2). The DMS footprints for the top and bottom strands are shown in Fig. 2 (a and b) and summarized in Fig. 2c. The two strands of the loxP and lox4 sites differed in the sensitivity of their spacer region to DMS modification upon binding to Cre, and the differences were dependent on the location of the scissile G. The scissile G at position 4 on the bottom strand of loxP (Fig. 2b, lane 2) was protected by the wild type Cre protein, and this protection was dependent on Lys-86 (see below). Bases 3' to the scissile G on the bottom strand of loxP (most notably position 3'A) exhibited enhanced reactivity to DMS upon binding to Cre. These enhanced sensitivities suggested that Cre causes structural alterations in the lox spacer region, possibly due to DNA bending (see "Discussion").

Cre also protected the scissile G and enhanced the DMS methylation of downstream bases (2'G, 1A, 3G, and 5T) on the top strand of lox4 (Fig. 2a, lane 6). The increased sensitivity to DMS was stronger on the top strand of lox4 than on the bottom strand of loxP. Because T residues are normally not methylated by DMS, the enhanced cleavage at position 5T on the top strand of lox4 suggests a significant conformational alteration near this position. Cre also caused hypermethylation of the scissile 4A on the bottom strand of lox4, whereas the corresponding scissile 4'A on the top strand of loxP did not exhibit increased reactivity to DMS. Using 1,10-phenanthroline-copper (OP-Cu), we also detected subtle differences in the spacer regions of loxP and lox4 when incubated with Cre (see Supplementary Data, Fig. S1). The Cre-induced DMS and OP-Cu footprinting patterns in the symmetry elements were similar for loxP and lox4 (Figs. 2 and S1) and generally agree with the contacts seen in the Cre-lox crystal structures (21, 2327).

Lys-86 Contacts the Scissile G in Both loxP and lox4 —Because Lys-86 influences the order of strand exchange (32, 33), we investigated the role of Lys-86 in the conformation of the Cre-lox complexes by mutating Lys-86 to Ala (K86A). Although the Cre K86A mutant protein is fused to a N-terminal His10 tag, the His-tagged wild type (HisWt) Cre protein exhibited similar properties to the non-tagged wild type Cre protein (Figs. 2 and S1) (32, 33).

We found that the K86A mutation relieved the protection of the scissile G from DMS methylation in both loxP (Fig. 2b, lane 4) and lox4 (Fig. 2a, lane 8). This supports the observation from several Cre-lox crystal structures that Lys-86 contacts the scissile G nucleotide (2325). The K86A mutation not only relieved the protection but strongly increased the reactivity of the scissile G to DMS methylation compared with the substrate without any protein added. The hyper-reactivity of the scissile G in the presence of Cre K86A required DMS treatment (data not shown), indicating that the hyper-reactivity was not due to Cre-catalyzed cleavage. The Lys-86-dependent protection of the scissile G was also observed in the DMS footprint of the substrate "se b" that contains only symmetry element b (Table I; data not shown). This suggests that the interaction between Lys-86 and the scissile G nucleotide does not require the binding of a Cre dimer. Bases 3' to the scissile G on the bottom strand of loxP and on the top strand of lox4 remained sensitive to DMS methylation when bound to Cre K86A. Furthermore, the Cre K86A mutant protein generally increased the sensitivity of the central region of the loxP and lox4 spacer to OP-Cu compared with the wild type protein (Fig. S1). The Cre K86A-induced DMS and OP-Cu footprinting patterns in the symmetry elements were similar to the wild type Cre protein.

Cre Bends loxP Toward the Major Groove Near the Middle of the Spacer Region—Because the chemical footprinting analyses showed evidence for Cre-induced conformational changes in the spacer region of the lox site, we used circular permutation and phasing analyses to detect the Cre-induced DNA bending in the lox site. Bent DNA migrates more slowly than unbent DNA in a polyacrylamide gel due to its shorter end-to-end distance (46).

We first used circular permutation analysis to map the position and angle of the Cre-induced DNA bends (Table II, Fig. 1(b and c), and Supplementary Data, Figs. S2–S6) (46, 47). The binding of a Cre monomer (cI) induces a bend of 35 ± 5° in the symmetry element generally near the margin of the spacer region ("cI bend"). The binding of two Cre molecules (cII) induces a bend of 55 ± 5° in the lox site near the center of the spacer region ("cII bend"). However, the small bend angles made it difficult to position the bend center accurately (see Supplementary Data, Fig. S7) and to interpret the significance of small differences between loxP and lox4. In addition, circular permutation cannot distinguish between a directed bend and a non-directed protein-induced flexure of the DNA (35, 47, 48).



View larger version (35K):
[in this window]
[in a new window]
 
FIG. 3.
Phasing analysis of the Cre-induced bends in loxP and lox4. a, schematic diagram of the phasing substrates. A 90-bp DNA containing the lox site (blue box) flanked by non-lox sequences (striped boxes) was inserted into the BamHI site of a set of phasing vectors (pK10, 12, 14, 16, 18, and 20) constructed by Zinkel and Crothers (35) (see "Experimental Procedures"). The phasing vectors were named such that the number refers to the length of the variable region (10–20 bp). Hence, the linker length (the distance between the center of the kinetoplast DNA that contains the A-tracts (35) and the middle of the lox site) was varied in 2-bp increments (111–121 bp) by one helical turn. The small arrow within the blue box indicates the orientation of the inserted 90-bp fragment as well as the lox site: the forward (F) orientation (<-) corresponds to the orientation of the lox site as shown in Fig. 1b and in Table I. Rs, RsaI; S, SacI; B, BamHI; N, NcoI; P, PvuII. b, gel mobility of the Cre-DNA complexes bound to the phasing substrates containing loxP (lanes 1–6) and lox4 (lanes 7–12) inserted in the forward orientation. 2 nM 32P-labeled DNA substrate was incubated with 0.05 µM wild type Cre protein and analyzed on a 4% native polyacrylamide gel as described under "Experimental Procedures." cI and cII, one or two Cre molecules bound to loxP, respectively (45); M, a 464-bp RsaI vector fragment that was contaminated with the 347- to 357-bp RsaI/PvuII phasing substrates (S). The M fragment was not bound by Cre (data not shown) and served as a useful non-phased size marker. Note that the 347- to 357-bp phasing substrates S migrated more slowly than the larger M fragment on the polyacrylamide gel likely as a consequence of its intrinsic DNA bend. This mobility anomaly was not observed on an agarose gel (data not shown). The cis isomer (the A-tract bend and the Cre-induced bend are in phase) and the trans isomer (the two bends are out of phase) are diagrammed to the left of the gel. The Cre molecule is represented as an oval. The asterisk marks the lane with the slowest relative mobility (cis isomer) for cII (see d). c, graph of the relative mobility of Wt Cre-loxP cI (open triangles and dashed line) and cII (filled squares and solid line) as a function of the linker length. The loxP site is inserted in the forward orientation. The vertical arrows denote the linker length at which the Cre-induced bend and A-tract bend are in cis (minimum relative mobility). The error bars represent the experimental errors from at least four experiments. The graphs were plotted to the best fit polynomial curve (Microsoft Excel). d, comparison of the relative mobility of cIIs formed by Wt Cre bound to loxP and lox4 in the forward (F) and reverse (R) orientations. loxP-F, filled squares and solid blue line; loxP-R, open diamonds and dashed green line; lox4-F, filled triangles and solid red line; lox4-R, open circles and dashed magenta line. e, graph of the relative mobility of Wt Cre-lox4 cI (open circles and dashed line) and cII (filled triangles and solid line) as a function of the linker length. The lox4 site is inserted in the forward orientation. The vertical arrows denote the linker length at which the Cre-induced bend and A-tract bend are in cis.

 
We therefore turned to phasing analysis to magnify the differences in electrophoretic mobility as well as to determine the direction of the Cre-induced bends. The Cre-induced bend was analyzed relative to a sequence-directed A-tract bend (Fig. 3) (35). We varied the helical phasing of the two bends (and hence the end-to-end distance of the DNA fragment) by varying the linker length between the lox site and the A-tract in 2-bp increments through one helical turn (Fig. 3a). When the two bends are in-phase (in the same direction), the end-to-end distance would be the shortest in the cis isomer, and so the protein-DNA complex would have the minimum mobility (Fig. 3b) (35, 46). In contrast, maximum mobility is obtained in the trans isomer when the two bends are out of phase (in the opposite direction). We define the linker length as the number of base pairs between the center of the kinetoplast DNA that contains the A-tracts (35) and the middle of the lox site. We chose the middle of the lox site as an arbitrary reference center for comparing the various lox sites. The sequences of the various lox sites studied are listed in Table I, and the results of the phasing analysis are summarized in Table II.

We found that Cre induced a directed bend in the loxP site. The electropherogram of wild type Cre bound to the loxP-containing phasing fragments is shown in Fig. 3b (lanes 1–6). The unbound substrates themselves exhibited a slight phasing in their mobility, suggesting that the lox-containing DNA possessed an intrinsic bend. Both cI and cII exhibited phase-dependent variations in gel mobility with a periodicity of about 10 bp indicative of a directed bend (Fig. 3c). The amplitude for the cII phasing curve was about 2-fold greater than that obtained for cI (cI, 0.04 versus cII, 0.08), suggesting that cII exhibited a larger bend angle than cI and supporting the relative bend angles measured by circular permutation analysis (Fig. S2c and Table II). The relative mobility of the Cre-loxP cI and cII reached a minimum mobility at a linker length (Lmin) of 116 and 115 bp, respectively (Fig. 3c), both of which correspond to about 11 helical turns with a helical repeat of 10.5 bp/turn. Because the A-tract is bent toward the major groove at the center of the kinetoplast DNA and the two in-phase bends are separated by an integral number of turns, we concluded that the Cre-induced bends in cI and cII were also toward the major groove in the middle of the loxP site (the arbitrarily designated reference center). This is further confirmed by analysis of the trans isomer in which the two out-of-phase bends were separated by ~121 (cI) or 120 bp (cII), corresponding to 11.5 turns (Fig. 3c). From the circular permutation analysis (Fig. S2c), we estimated that the Cre-induced bend center in cI is close to position 4, about half a helical turn from the reference phasing center, and so, Cre bends the loxP in cI toward the minor groove at position 4. In contrast, the bend center in cII is near position 1' (close to the reference phasing center), and so, the loxP in cII is bent toward the major groove near the middle of the spacer region.

Cre Induces an Asymmetric Bend in loxP—An asymmetrically positioned DNA bend was observed in the crystal structure of the Cre-lox synaptic complex (25). We assessed the asymmetry of the Cre-induced bend in the loxP site by comparing the relative phasing of the loxP site in the forward and reverse orientations. The 90-bp test fragment containing the loxP site precisely in its middle was inserted into the phasing vectors in either the forward (F) or the reverse (R) orientation (Fig. 3a). Therefore, the linker length between the A-tract and the middle of the lox site was the same in both orientations. Because the position and magnitude of the Cre-induced bends should be identical regardless of the orientation of the fragment, differences in the relative mobility between the two orientations would reflect a difference in the position of the actual Cre-induced bend center relative to the A-tracts. If the center of the Cre-induced bend is exactly in the middle of the lox site (i.e. at the reference center) and the DNA is bent toward one of the grooves, then the phasing curves should be the same for both orientations (46). On the other hand, if Cre induces an asymmetric DNA bend, then the actual distance between the two bends would be different in the two orientations and hence, the phasing curves would be altered (unless the difference is a full helical turn).

We found that wild type Cre does indeed bend the loxP site asymmetrically and that the asymmetric bend is dictated by the asymmetric spacer sequence. The gel mobility shown in Fig. 3 (b and c) was for the forward orientation, "loxP-F," which corresponds to the orientation of the loxP site illustrated in Fig. 1b. In Fig. 3d, we compare the relative mobility of the Cre-loxP cII bound to the "loxP-F" substrates and the reverse "loxP-R" substrates. We found that the two phasing curves were slightly different: the "loxP-R" curve was consistently shifted by about +1 bp relative to the "loxP-F" curve. The "loxP-F" and "loxP-R" cII had Lmin values of 115 and 116 bp, respectively, both of which still correspond to a bend toward the major groove in the middle of the loxP sequence (11 turns, 10.5 bp/turn). Therefore, reversing the orientation of the loxP site causes a small but highly reproducible shift in the phasing curves, suggesting that the Cre-induced bend in the cII was asymmetrically positioned in the loxP site. Asymmetry in the cI bend was also seen (Table II) and is discussed in more detail in a later section.

If the apparent asymmetry in the Cre-induced DNA bend is indeed due to the asymmetric loxP site and not to differences in the flanking non-lox sequences, then a symmetric lox site should abolish this asymmetry and the phasing curves would be the same in both the forward and reverse orientations. In fact, we found that the Cre-induced bend was indeed positioned symmetrically in both the symmetric loxSA and loxSB substrates (Table I), with an Lmin of around 115.5 bp in both orientations for cI and cII (Table II and Fig. S8). Therefore, the asymmetric Cre-induced bend in the loxP site is dictated by the asymmetric spacer sequence.

The Scissile Base Pairs Determine the Cre-induced DNA Bends: the Cre-induced Bends in lox4 Resemble Those for loxP in the Reverse Orientation—The DMS footprint induced by Cre in the lox4 spacer region was almost the reverse of that for loxP (Fig. 2). To investigate whether the Cre-induced bends in loxP and lox4 exhibit a similar inverse relationship, we also constructed two sets of lox4 phasing substrates with the lox4-containing DNA inserted in the forward ("lox4-F") and reverse ("lox4-R") orientations. The electrophoretic mobility shift assay of the "lox4-F" phasing substrates is shown in Fig. 3b (lanes 7–12) and the phasing curves are plotted in Fig. 3e. We observed differences between the phasing of the Cre-lox4 and the Cre-loxP complexes with the most notable difference being in cII (Fig. 3d). We discuss the differences in cI at a later section. Although the cII bend centers in loxP and lox4 differed by only 1 bp by circular permutation (Figs. 1b, 1c, S2, and S3), we confirmed by phasing analysis that they were indeed different and asymmetric (Fig. 3d). Recall that the relative mobility of the "loxP-F" cII had a Lmin of 115 bp (Fig. 3b, lane 3). In contrast, the "lox4-F" cII reached an apparent minimum gel mobility at a linker length of 117 bp (Fig. 3b, lane 10). After we corrected for the intrinsic bend in the substrate, the phasing curve for the "lox4-F" cII reproducibly had a Lmin of 116 bp, a +1 bp shift relative to "loxP-F" (Fig. 3d). Therefore, although Cre appears to bend the lox4 site in the same general direction as loxP (i.e. toward the major groove at the center of the lox4 site), the 1-bp difference in phasing may reflect the 1-bp difference in the position of the bend center measured by circular permutation (Fig. 1, b and c). The Cre-induced DNA bend in lox4, like that in loxP, was asymmetrically positioned, because the cII phasing curve for "lox4-R" was shifted by almost –1.5 bp relative to "lox4-F" (Fig. 3d). Interestingly, the phasing curve for "lox4-F" cII coincided with that for loxP in the reverse orientation ("loxP-R") and correspondingly, the "lox4-R" curve resembled the "loxP-F" curve. Therefore, the effect of interchanging the 4' and 4 bp is analogous to reversing the entire loxP site.

Single Symmetry Elements and cI: Wild Type Cre Bends loxP and lox4 Near the Scissile G Base Pair in cI—We wished to examine whether the asymmetric lox spacer sequence also influences the Cre-induced bend in cI. Because cI formed on a full lox site may consist of a mixture of complexes of a Cre molecule bound to either the left or the right symmetry element, we constructed lox half-sites in which we replaced one or the other symmetry element with a random sequence to disrupt Cre binding (Table I and Fig. 4a). The "se a" substrates contain an intact left symmetry element a and spacer region, but the right symmetry element had been disrupted. Similarly, the "se b" substrates have a mutated left element and an intact right symmetry element b. Note that symmetry elements a and b have identical sequence and are defined by their position relative to the asymmetric spacer sequence (Fig. 1b).



View larger version (31K):
[in this window]
[in a new window]
 
FIG. 4.
Phasing analysis of the Cre-induced bend in the single symmetry element substrates se a and se b. a, a schematic of the phasing substrates se a and se b in the forward (F) and reverse (R) orientations. The se a and se b sequences are listed in Table I. The disrupted symmetry element (X) was generated by replacement with a random sequence to allow analysis of cI bound to the intact symmetry element. The functional symmetry element in the "se a-F " and the "se b-R" substrates is proximal to the A-tracts than is the functional symmetry element in "se a-R" and "se b-F " (distal). b, EMSA of wild type Cre bound to the "se a-F " (lanes 1–6) and the "se b-F " (lanes 7–12) phasing substrates. The linker length (bp) is indicated above each lane. The asterisk marks the lane with the slowest relative mobility (cis isomer) for cI (see c). c, comparison of the relative mobility of cI bound to the se a and se b phasing substrates in the forward and reverse orientations. Note that the phasing profile of the "se a-F " cI (filled triangles and red line) is similar to that of "se b-R" cI (open diamonds and green dashed line), and, conversely, the distal pairs "se a-R" (open circles and dashed magenta line) and "se b-F " (filled squares and blue line) are similar to each other. The vertical arrows denote the linker length at which the Cre-induced bend and A-tract bend are in cis. d, comparison of the phasing curves of wild type Cre cI bound to the loxP, lox4, se a, and se b in the forward orientation. loxP, open diamonds and green dashed line; lox4, open circles and dashed magenta line; se a, filled triangles and red line; se b, filled squares and blue line.

 

We detected a difference in the phasing of the cI bound to "se a-F " and "se b-F " (forward orientation): the Lmin for "se a-F " was at 115.1 bp, whereas that for "se b-F " was at 116.4 bp (Fig. 4b, lanes 3 and 10, and Fig. 4c). The distance of Lmin was calculated using the middle of the spacer region as the reference center, which does not correspond to the bend centers determined by circular permutation: position 8' for se a and position 6 for se b (Fig. S4). For example, the bend center in se b is separated by (115.1 + 5.5) or ~121 bp from the A-tracts in the cis isomer, and 121 bp corresponds to 11.5 turns with a helical repeat of ~10.5 bp/turn, implying that Cre bends se b toward the minor groove at position 6.

The se a-F and se b-F phasing substrates differ in the distance of the functional symmetry element relative to the A-tracts: symmetry element a is closer to the A-tracts than the more distal symmetry element b in the forward orientation (Fig. 4a). The 1-bp difference in the Lmin for the phasing curve of se a-F and se b-F may simply reflect this difference in relative distance. To better evaluate the two symmetry elements, we compared substrates in which the distance of the functional symmetry element relative to the A-tracts was the same, but the orientation of the spacer sequence was reversed. For example, the functional symmetry element in the "se a-F " and the "se b-R" substrates are proximal (closer) to the A-tracts than the functional symmetry element in the "se a-R" and "se b-F " substrates, which are more distal to the A-tracts (Fig. 4a). The relative mobility profile of the "se a-F " substrates was similar to the proximal "se b-R" substrates, and conversely, the distal pairs "se a-R" and "se b-F " were similar to each other (Fig. 4d). Therefore, Cre appears to bend se a and se b toward the minor groove near the margin of the spacer region independently of the spacer sequence.

However, in a full lox site with two competing symmetry elements, the spacer sequence influences the predominant cI species. The cI bends induced by wild type Cre in loxP and lox4 mapped to opposite ends of the spacer region by circular permutation (Fig. 1, b and c, and Supplemental Figs. S2 and S3) and phasing analyses (Fig. 4d). The cI bend in the loxP site was at the scissile G adjacent to symmetry element b, and its phasing curve resembled that of the se b substrate more closely than that of the se a substrate. In contrast, the cI bend in the full lox4 site was in symmetry element a (also adjacent to the scissile G), and its phasing curve resembled that of the se a substrate. The differences in the cI phasing between the forward and reverse orientations were smaller (~0.5 bp) for the full lox substrates than the individual single symmetry elements (≥1 bp; Table II). This may have been due to the shuffling of the Cre monomer between the left and the right symmetry elements during electrophoresis (37, 49). In conclusion, the wild type Cre bends loxP and lox4 near the scissile G base pair in cI, and this suggests preferential binding of the first Cre molecule to the adjacent symmetry element (see "Discussion").

Lys-86 Is Important for the Positioning of the Cre-induced DNA Bend in loxP and lox4 cI—Because Lys-86 contacts the scissile bases (21, 2327) and affects the order of strand exchange (32, 33), we investigated its role in the position and direction of the Cre-induced bends. Mutation of Lys-86 to Ala (K86A) dramatically altered the position of the DNA bend in cI, shifting the bend center (determined by circular permutation) by 11 bp from position 4 to position 8' in loxP (Figs. 1b and S5b) and shifting the Lmin in the phasing analysis by –1 bp to a linker length of about 115 bp (Fig. 5b). Note that the phasing curve for Cre K86A-loxP cI resembled that of the substrate se a (Table II), positioning the Cre K86A-induced cI bend in loxP in symmetry element a adjacent to the scissile AT base pair. The K86A mutation did not significantly affect the position of the cI bend in the single symmetry element substrates se a and se b (Fig. 1b and Table II), suggesting that Lys-86 is not involved in actually bending the DNA but rather in selecting the symmetry element to be occupied by the first Cre monomer within a full lox site (see "Discussion"). The K86A-loxP cII exhibited a similar phasing profile to wild type Cre-loxP cII with a Lmin of about 115 bp (Fig. 5c), suggesting that both these proteins induced a similar cII bend in loxP. Circular permutation analysis revealed that Cre K86A also bends loxP at similar position to the wild type Cre protein in cII (Figs. 1b and S5c). Therefore, Lys-86 is important for positioning the cI bend in loxP but has little if any role in directing the cII bend.



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 5.
Phasing analysis of the Cre K86A-induced bends in loxP. a, EMSA of His-tagged wild type Cre (HisWt; lanes 1–6) or HisK86A (lanes 7–12) Cre protein bound to loxP-F phasing substrates. b, comparison of the phasing curves of loxP-F cI bound by HisWt and HisK86A. The arrows indicate the linker length at which the minimum mobility is reached. c, comparison of the phasing curves of loxP cII bound by HisWt and HisK86A. For b and c: HisWt, open squares and dashed line; HisK86A, filled circles and solid line. Note that the phasing curves for HisWt are similar to those of the untagged Wt Cre protein (Fig. 3c). The His-tagged proteins gave two cI bands, although HisWt gave predominantly the lower cI band (see a). The two cI bands may be due to the His tag for the following reasons. (i) They were also observed using substrates containing only a single symmetry element (se a and se b, data not shown) and so did not appear to correspond to the binding of a Cre monomer to different symmetry element within the lox site. (ii) The two cIs exhibit similar circular permutation and phasing profiles (data not shown), suggesting that their bend centers are located at similar positions. (iii) Analysis of the complexes by two-dimensional electrophoresis (first, native PAGE; second, SDS-PAGE) did not detect a covalently attached Cre in either cIs (data not shown). We show in b the phasing curves for the lower cI, but the upper cI band exhibits similar phasing profiles (data not shown).

 

In contrast to loxP, Lys-86 influences both the cI and cII bends in lox4 (Figs. 6 and S6). As in loxP, the cI bend centers for the wild type and K86A Cre protein were located at opposite ends of the lox4 spacer region (Wt Cre, at position 5' in the left symmetry element; Cre K86A, position 4 near the right margin of the spacer; Figs. 1c and S6b). The Lmin for the phasing curve of the Cre K86A-lox4 cI was 116 bp, a difference of ~1 bp from that for the wild type Cre-lox4 cI (Fig. 6b) and close to that for the substrate se b (Table II). Hence, Cre K86A bends both the loxP and lox4 sites near the scissile AT base pair, whereas wild type Cre bends the DNA near the scissile GC base pair in cI.



View larger version (31K):
[in this window]
[in a new window]
 
FIG. 6.
Phasing analysis of the Cre K86A-induced bends in lox4. a, EMSA of HisK86A Cre protein bound to lox4-F phasing substrates. b, comparison of the phasing curves of lox4-F cI bound by HisWt and HisK86A. The arrows indicate the linker length at which the minimum mobility is reached. c, comparison of the phasing curves of lox4-F cII bound by HisWt and HisK86A. For b and c: HisWt, open squares and dashed line; HisK86A, filled circles and solid line. Note that the phasing curves for HisWt are similar to those of the untagged Wt Cre protein (Fig. 3e). d, comparison of the phasing curves of HisK86A-lox4 cII in the forward (F) and reverse (R) orientations. F, filled circles and solid line; R, open triangles and dashed line.

 

Unlike loxP, the K86A mutation also altered the cII bend in lox4, shifting the phasing curve by –0.4 bp relative to wild type Cre with the lox4-F substrates (Fig. 4c). Although the difference in the cII bend between the wild type and K86A Cre proteins is small, this difference was also evident and slightly larger (~0.7 bp) in the lox4-R substrates (Table II). The K86A-lox4 cII appeared to be bent symmetrically around the middle of the lox4 site, analogous to the bends in the symmetrical loxSA and loxSB sites (Fig. S7). The difference in Lmin for the Cre K86A cII bound to the lox4-F and lox4-R substrates was only ~0.3 bp (Fig. 6d), which is smaller than the 1.4-bp difference in Lmin between the two orientations for the wild type Cre-lox4 cII (Fig. 3d). This symmetric cII bend induced by Cre K86A correlates with its loss of strand preference in the initiation of strand exchange in lox4 (see "Discussion") (32).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
The Scissile Base Pairs Determine the Cre-induced DNA Bends in cI and cII—The Cre recombinase has long been known to initiate recombination by preferentially exchanging the bottom strands of the loxP site and resolving the Holliday intermediate on the top strands (11, 2932). We previously found that the scissile base pairs at positions 4' and 4 were critical in determining this order of strand exchange (31, 32). The evidence we report here links the scissile base pairs to their influence on the Cre-induced DNA bends within cI and cII.

In the Cre-lox crystal structures, the two DNA strands adopt different conformations and are contacted differently in the spacer region by the cleaving and non-cleaving Cre subunits (21, 2327). The "continuous" strand contains the inactive scissile phosphate and is on the concave side of the DNA bend (Fig. 1a) (25). It makes numerous direct contacts with the N-terminal domain of the cleaving Cre subunit directly opposite the activated scissile phosphate. The more extended "crossing" strand contains the activated scissile phosphate and is contacted directly by the Cre molecules only at the activated scissile phosphate and scissile nucleotide but not throughout the rest of the spacer region (25). The crossing strand points toward the center cavity of the synapse ready for strand exchange after cleavage (Fig. 1a).

The major differences in the Cre-induced footprints of loxP and lox4 were in the spacer region (Fig. 2). The footprints within the lox4 site were similar to those for the loxP site in the reverse orientation. Cre increased the sensitivity to DMS methylation of the bases in the spacer region on the strand containing the scissile G nucleotide in both loxP and lox4. We propose that the DMS-sensitive, scissile G-containing strand corresponds to the more exposed "crossing" strand in both loxP and lox4 complexes (Fig. 7, see below).



View larger version (25K):
[in this window]
[in a new window]
 
FIG. 7.
Model relating the Cre-induced DNA bending and the order of strand exchange. The scissile base pairs dictate the order of strand exchange by influencing the position of the Cre-induced DNA bends in the loxP (a) and lox4 sites (b). We propose that the first Cre monomer is directed via Lys-86 to the symmetry element adjacent to the scissile G and induces a bend near the margin of the spacer region. The binding of the second Cre monomer shifts the bend toward the center of the spacer region with a slight bias toward the scissile A. The cleaving strand is on the more exposed convex side of the DNA bend (i.e. incipient "crossing" strand). One of the two Cre molecules bound to the lox site ("cleaving" subunit, dark oval) is poised to cleave the crossing strand at the scissile G (arrow), whereas the other "non-cleaving" monomer (gray) is in the inactive conformation (see Fig. 1a) (21, 23, 25). The scissile G-containing crossing strand is the bottom strand (thin line) of loxP (a) or the top strand (thick line) of lox4 (b). Interchanging the scissile base pairs to create the mutated lox4 site changes the position of the asymmetric cI and cII DNA bends and, as a consequence, alters the site of initial strand exchange (arrow). Note that the DNA bends in the Cre-lox4 complexes resemble those in which the loxP site is inverted by 180°.

 

The conformational alterations seen in the chemical footprints are supported by changes in the Cre-induced DNA bending measured by circular permutation and phasing analyses (Fig. 1, b and c, and Table II). The cI bends induced by the wild type Cre protein in loxP and lox4 were both situated near the scissile G base pair where recombination initiates (11, 2932). Likewise we demonstrated that the cII bends in the loxP and lox4 sites were different and asymmetric (Fig. 3d). The cII bends induced by the wild type Cre protein in loxP and lox4 were near the middle of the spacer region with a slight bias toward the scissile A. In fact, interchanging the scissile base pairs as in the lox4 site resulted in changes similar to those seen by inverting the entire loxP site by 180°. Therefore, the scissile base pairs dictate the Cre-induced DNA bending and the order of strand exchange (32).

Model for DNA Bending and the Order of Strand Exchange—We have demonstrated that the Cre-induced DNA bends in loxP and lox4 correlate with our previously observed alteration in the site of initial strand exchange (32), confirming the original proposal by Guo et al. (25). We present in Fig. 7 a revision of the model proposed by Guo and colleagues (25). We propose that the first Cre molecule binds preferentially to the symmetry element adjacent to the scissile G and induces a cI bend near the margin of the spacer. We speculate that because this Cre molecule is bound adjacent to the initial cleavage site, it may be in a "pre-activation/cleaving" conformation even prior to the binding of the second Cre monomer. The binding of a second Cre molecule redirects the bend asymmetrically toward the middle of the spacer region away from the scissile G. Within the cII and/or synaptic complex, the strand containing the scissile G adopts the "crossing" configuration (as suggested by the footprinting analysis) and is cleaved adjacent to the scissile G. Interchanging the scissile base pairs as in the lox4 site reverses the DNA bend induced by Cre from those in loxP and, as a consequence, alters the order of strand exchange (32).

The proposed preferential binding of Cre to a particular symmetry element in cI remains to be verified. Kinetic and equilibrium studies indicate that Cre binds the individual symmetry elements a and b with similar affinity (50). In addition, we did not detect preferential protection of one symmetry element relative to the other element on a full lox site by footprint analyses of the isolated cI (data not shown).

We found that Cre induced a 35° cI bend in both loxP and lox4; this small cI bend may correspond to the ~25° smooth bend observed within each symmetry element in the crystal structures of the Cre-lox synaptic complex (21, 2327). We measured a cII bend of ~55° toward the major groove at position 1' in loxP and at position 1 in lox4. Our cII bends resemble most closely the DNA bend observed in the Cre-lox crystal structures for the covalent cleavage intermediate and the Holliday intermediate in which the DNA was bent near the center of the spacer region toward the major groove (23, 24). However, our results differ to some extent from the asymmetric DNA bend observed in the crystal structure of the Cre-loxS synaptic complex in which the DNA was bent 78° toward the minor groove between positions 2/3 (or 2'/3') (25). It is possible that the DNA bend angles within a cII in solution and a crystallized synaptic complex are different. The latter consists of two lox sites and four Cre molecules. The act of synapsis and/or crystallization may alter the bend angle and position. Furthermore, because the bend angle measured by circular permutation analysis was small, the positioning of the bend was subject to appreciable error (see Supplementary Data, Fig. S7). The circular permutation analysis may also underestimate the magnitude of the Cre-induced DNA bend (48, 51). Other possible reasons for the discrepancies include the type of Cre proteins and lox sites used. The crystal structures were solved using mutant Cre R173K and Y324F proteins to block cleavage (25), whereas we used wild type and K86A Cre proteins. In addition, Guo et al. (25) used a symmetric loxS, which corresponds to our loxSB sequence (Table I). Although we did not determine the position of the cII bend center in our symmetrical lox substrates, we observed symmetric bends in these substrates by phasing analysis. It is possible that the symmetric phasing curves represent an average of rapidly alternating, asymmetric bends as seen in the crystal structures (25). Correspondingly, the Cre-loxP and Cre-lox4 cII may also consist of rapidly alternating isomers whose bend center and direction may not correspond precisely to the overall average measured. Nevertheless, the predominant cII species is different for loxP and lox4.

Lys-86 and DNA Bending—The scissile base pairs of the lox spacer region are critical in determining the position of the Cre-induced bends and the order of strand exchange, but how does Cre distinguish between the scissile base pairs? We confirmed by DMS methylation that Lys-86 makes major groove contacts with the scissile G nucleotide in both loxP and lox4 as observed in several of the Cre-lox crystal structures (2325). We observed this contact even when only symmetry element b was present (data not shown), indicating that the Lys-86-scissile G interaction can occur within cI. In a recent crystal structure of the wild type Cre protein bound to a loxP-Holliday junction, Lys-86 of the cleaving Cre subunit contacts the scissile A at N7 and the adjacent scissile phosphate was activated for cleavage (27). However, Lys-86 of the non-cleaving Cre subunit did not contact the adjacent scissile G. Because resolution occurs preferentially adjacent to the scissile A (11, 31, 32), this structure may represent the conformation of the Cre-lox synaptic complex just prior to resolution, and this may differ from the initial cII or the synaptic complex containing two linear lox sites.

Mutation of Lys-86 to Ala significantly altered the position of the cI bend, shifting the bend center by 8–11 bp to position 8' in loxP and to position 4 in lox4 (Fig. 1, b and c, and Table II). The cI bends induced by Cre K86A in both loxP and lox4 were near the scissile AT base pair and were at the opposite half of the lox site relative to the cI bends induced by the wild type Cre protein. Mutation of Lys-86 to Ala did not significantly alter the position of the cI bend in the single symmetry element substrates se a and se b. We suggest that Lys-86 is not involved in actually bending the DNA, but is responsible for directing the Cre-induced cI bend and possibly the first Cre molecule toward the scissile GC base pair on a full lox site. The interaction between Lys-86 and the scissile G may help stabilize the cI in which the wild type Cre protein is bound adjacent to the scissile G.

Although important for positioning the cI bend, Lys-86 was not critical for establishing the cII bend in loxP (Table II). Despite the fact that the wild type and K86A Cre proteins bent loxP at opposite ends of spacer region in cI, they induced similar cII bends. This implies that, unlike wild type Cre, the first or predominant Cre K86A molecule to bind to the loxP site in cI does not necessarily become the cleaving subunit in cII and/or the synaptic complex. The similarities in the circular permutation and phasing analyses of the loxP cII bends induced by the wild type and K86A Cre proteins further complement the footprinting results. Like the wild type Cre protein, Cre K86A also increased the sensitivity of the spacer region on the bottom strand of loxP to DMS and OP-Cu modifications (Figs. 2 and S1). This could imply that the bottom strand is the incipient crossing strand in the wild type and K86A Cre-loxP complexes as discussed above (Fig. 7). The K86A mutation may not have dramatically altered the DNA conformation in cII but rather simply revealed the already exposed crossing strand, and this would account for the hypermethylation of the scissile G. Therefore, although the scissile base pairs are still critical for directing the cII bend in loxP, it appears that their interaction with Lys-86 is not essential. This is consistent with the fact that Lys-86 is not important for determining the strand preference in loxP during the first strand exchange event (32).

Although Lys-86 did not determine the cII bend in loxP,itdid contribute to the cII bend in lox4. Whereas the wild type Cre protein induced an asymmetric cII bend in lox4, the cII bend induced by the Cre K86A protein was more symmetrical. The symmetric K86A-lox4 cII bend is consistent with our previous observation that the K86A mutation abolished the strand bias for the initiation of strand exchange in lox4 (32). Therefore, Lys-86 may be responsible for establishing the asymmetric cII bend and the site of initial strand exchange in lox4.

In summary, the location of the scissile base pairs is important for dictating the asymmetric Cre-induced cI and cII bends in loxP and lox4 as well as the site of the first strand exchange event. Lys-86 is critical for directing the cI bend in loxP and lox4 toward the scissile GC base pairs as well as the position of the cII bend and the site of initiation in lox4. In contrast, the cII bend and the site of initiation in loxP are not dependent on Lys-86. Additional bases in the lox spacer region as well as other residues in Cre may also contribute to the Cre-induced DNA bending and the order of strand exchange.


    FOOTNOTES
 
* This work was supported in part by grants from the Medical Research Council (MRC) of Canada and the Canadian Institutes for Health Research (CIHR) (to P. D. S.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Back

The on-line version of this article (available at http://www.jbc.org) contains Supplemental Figs. S1–S8. Back

{ddagger} Supported by Studentships from the MRC and CIHR and by a University of Toronto Open Fellowship. Back

§ To whom correspondence should be addressed. Tel.: 416-978-6061; Fax: 416-978-6885; E-mail: p.sadowski{at}utoronto.ca.

1 The abbreviations used are: cI, complex I; cII, complex II; DMS, dimethyl sulfate; EMSA, electrophoretic mobility shift assay; OP-Cu, 1,10-phenanthroline-copper; Wt, wild type. Back


    ACKNOWLEDGMENTS
 
We thank Barbara Funnell, Brigitte Lavoie, and Linda Beatty for helpful comments.



    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Austin, S., Ziese, M., and Sternberg, N. (1981) Cell 25, 729–736[Medline] [Order article via Infotrieve]
  2. Hoess, R. H., and Abremski, K. (1990) in Nucleic Acids and Molecular Biology (Eckstein, F., and Lilley, D. M. J., eds) Vol. 4, pp. 99–109, Springer-Verlag, Berlin
  3. Sadowski, P. D. (1993) FASEB J. 7, 760–767[Abstract/Free Full Text]
  4. Landy, A. (1993) Curr. Opin. Genet. Dev. 3, 699–707[Medline] [Order article via Infotrieve]
  5. Jayaram, M., Tribble, G., and Grainge, I. (2002) in Mobile DNA II (Craig, N. L., Craigie, R., Gellert, M., and Lambowitz, A. M., eds) pp. 192–218, ASM Press, Washington, D. C.
  6. Argos, P., Landy, A., Abremski, K., Egan, J. B., Haggard, L. E., Hoess, R. H., Kahn, M. L., Kalionis, B., Narayana, S. V., Pierson, L. d., Sternberg, N., and Leong, J. M. (1986) EMBO J. 5, 433–440[Abstract]
  7. Abremski, K. E., and Hoess, R. H. (1992) Protein Engineering 5, 87–91[Abstract]
  8. Nash, H. (1996) in Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology (Neidhardt, F., Curtiss, R., Ingraham, J., Lin, E., and Low, K., eds) pp. 2363–2376, ASM Press, Washington, D. C.
  9. Nunes-Düby, S. E., Kwon, H. J., Tirumalai, R. S., Ellenberger, T., and Landy, A. (1998) Nucleic Acids Res. 26, 391–406[Abstract/Free Full Text]
  10. Holliday, R. (1964) Genet. Res. 5, 282–304
  11. Hoess, R., Wierzbicki, A., and Abremski, K. (1987) Proc. Natl. Acad. Sci. U. S. A. 84, 6840–6844[Abstract]
  12. Nunes-Düby, S., Matsumoto, L., and Landy, A. (1987) Cell 50, 779–788[Medline] [Order article via Infotrieve]
  13. Kitts, P. A., and Nash, H. A. (1987) Nature 329, 346–348[CrossRef][Medline] [Order article via Infotrieve]
  14. Kitts, P. A., and Nash, H. A. (1988) Nucleic Acids Res. 16, 6839–6856[Medline] [Order article via Infotrieve]
  15. Meyer-Leon, L., Huang, L. C., Umlauf, S. W., Cox, M. M., and Inman, R. B. (1988) Mol. Cell. Biol. 8, 3784–3796[Medline] [Order article via Infotrieve]
  16. Jayaram, M., Crain, K. L., Parsons, R. L., and Harshey, R. M. (1988) Proc. Natl. Acad. Sci. U. S. A. 85, 7902–7906[Abstract]
  17. Arciszewska, L. K., and Sherratt, D. J. (1995) EMBO J. 14, 2112–2120[Abstract]
  18. Hoess, R. H., Ziese, M., and Sternberg, N. (1982) Proc. Natl. Acad. Sci. U. S. A. 79, 3398–3402[Abstract]
  19. Hoess, R., Abremski, K., and Sternberg, N. (1984) Cold Spring Harbor Symp. Quant. Biol. 49, 761–768[Medline] [Order article via Infotrieve]
  20. Hoess, R., Abremski, K., Irwin, S., Kendall, M., and Mack, A. (1990) J. Mol. Biol. 216, 873–882[Medline] [Order article via Infotrieve]
  21. Van Duyne, G. D. (2001) Annu. Rev. Biophys. Biomol. Struct. 30, 87–104[CrossRef][Medline] [Order article via Infotrieve]
  22. Hoess, R. H., and Abremski, K. (1985) J. Mol. Biol. 181, 351–362[Medline] [Order article via Infotrieve]
  23. Guo, F., Gopaul, D. N., and van Duyne, G. D. (1997) Nature 389, 40–46[CrossRef][Medline] [Order article via Infotrieve]
  24. Gopaul, D. N., Guo, F., and Van Duyne, G. D. (1998) EMBO J. 17, 4175–4187[Abstract/Free Full Text]
  25. Guo, F., Gopaul, D. N., and Van Duyne, G. D. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 7143–7148[Abstract/Free Full Text]
  26. Gopaul, D. N., and Van Duyne, G. D. (1999) Curr. Opin. Struct. Biol. 9, 14–20[CrossRef][Medline] [Order article via Infotrieve]
  27. Martin, S. S., Pulido, E., Chu, V. C., Lechner, T. S., and Baldwin, E. P. (2002) J. Mol. Biol. 319, 107–127[CrossRef][Medline] [Order article via Infotrieve]
  28. Hamilton, D. L., and Abremski, K. (1984) J. Mol. Biol. 178, 481–486[Medline] [Order article via Infotrieve]
  29. Hoess, R. H., Wierzbicki, A., and Abremski, K. (1990) in Structure & Methods Vol. 1: Human Genome Initiative & DNA Recombination (Sarma, R. H., and Sarma, M. H., eds) pp. 203–213, Adenine Press, Guilderland, NY
  30. Shaikh, A. C., and Sadowski, P. D. (2000) J. Biol. Chem. 275, 30186–30195[Abstract/Free Full Text]
  31. Lee, G., and Saito, I. (1998) Gene (Amst.) 216, 55–65[CrossRef][Medline] [Order article via Infotrieve]
  32. Lee, L., and Sadowski, P. D. (2003) J. Mol. Biol. 326, 397–412[CrossRef][Medline] [Order article via Infotrieve]
  33. Lee, L., and Sadowski, P. D. (2001) J. Biol. Chem. 276, 31092–31098[Abstract/Free Full Text]
  34. Zwieb, C., Kim, J., and Adhya, S. (1989) Genes & Dev. 3, 606–611[Abstract]
  35. Zinkel, S. S., and Crothers, D. M. (1987) Nature 328, 178–181[CrossRef][Medline] [Order article via Infotrieve]
  36. Shaikh, A. C., and Sadowski, P. D. (1997) J. Biol. Chem. 272, 5695–5702[Abstract/Free Full Text]
  37. Beatty, L. G., and Sadowski, P. D. (1988) J. Mol. Biol. 204, 283–294[Medline] [Order article via Infotrieve]
  38. Funnell, B. E., and Gagnier, L. (1993) J. Biol. Chem. 268, 3616–3624[Abstract/Free Full Text]
  39. Craig, N. L., and Nash, H. A. (1984) Cell 39, 707–716[Medline] [Order article via Infotrieve]
  40. Sigman, D. S., Kuwabara, M. D., Chen, C. B., and Bruice, T. W. (1991) Methods Enzymol. 208, 414–433[Medline] [Order article via Infotrieve]
  41. Wissmann, A., and Hillen, W. (1991) Methods Enzymol. 208, 365–379[Medline] [Order article via Infotrieve]
  42. Tullius, T. D., Dombroski, B. A., Churchill, M. E. A., and Kam, L. (1987) Methods Enzymol. 155, 537–558[Medline] [Order article via Infotrieve]
  43. Dixon, W. J., Hayes, J. J., Levin, J. R., Weidner, M. F., Dombroski, B. A., and Tullius, T. D. (1991) Methods Enzymol. 208, 380–413[Medline] [Order article via Infotrieve]
  44. Price, M. A., and Tullius, T. D. (1992) Methods Enzymol. 212, 194–219[Medline] [Order article via Infotrieve]
  45. Mack, A., Sauer, B., Abremski, K., and Hoess, R. (1992) Nucleic Acids Res. 20, 4451–4455[Abstract]
  46. Crothers, D. M., and Drak, J. (1992) Methods Enzymol. 212, 46–71[Medline] [Order article via Infotrieve]
  47. Wu, H. M., and Crothers, D. M. (1984) Nature 308, 509–513[Medline] [Order article via Infotrieve]
  48. Kerppola, T. K., and Curran, T. (1991) Science 254, 1210–1214[Medline] [Order article via Infotrieve]
  49. Luetke, K. H., and Sadowski, P. D. (1998) Nucleic Acids Res. 26, 1401–1407[Abstract/Free Full Text]
  50. Ringrose, L., Lounnas, V., Ehrlich, L., Buchholz, F., Wade, R., and Stewart, A. F. (1998) J. Mol. Biol. 284, 363–384[CrossRef][Medline] [Order article via Infotrieve]
  51. Hagerman, P. J. (1992) Biochim. Biophys. Acta 1131, 125–132[Medline] [Order article via Infotrieve]