![]()
|
|
||||||||
J. Biol. Chem., Vol. 281, Issue 37, 26893-26903, September 15, 2006
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



1
From the
Instituto Leloir, Patricias Argentinas 435, 1405 Buenos Aires, Argentina and the
Department of Medical Genetics and Microbiology, University of Toronto, Ontario M5S 1A8, Canada
Received for publication, March 6, 2006 , and in revised form, June 28, 2006.
| ABSTRACT |
|---|
|
|
|---|
41 min) through an undissociated intermediate species and the other corresponds to a fast twostep dissociation route (t
0.8 min). In line with this, at least two parallel association events from two populations of protein conformers are observed, with on-rates of 0.251.6 x 108 M1 s1, which occur differentially either in excess protein or DNA molecules. Both parallel complexes undergo subsequent firstorder rearrangements of
2.0 s1 to yield two consolidated complexes. These parallel association and dissociation routes likely allow additional flexible regulatory events for site recognition depending on site availability according to nucleus environmental conditions, which may lock a final recognition event, dissociate and re-bind, or slide along the DNA. | INTRODUCTION |
|---|
|
|
|---|
EBNA1 is the only protein expressed in all types of EBV latent infection (2). It binds to the EBV replication origin oriP where it plays several roles: initiation of DNA replication, segregation of EBV episomes, and transactivation of latent viral gene expression (35). All these functions require the binding of the EBNA1 DNA binding domain to specific 18-bp DNA recognition sites (6, 7). The crystal structure of this domain revealed a particular fold, the dimeric
-barrel, only shared with the papillomavirus E2 DNA binding domain (see Fig. 1) (8, 9). There is no amino acid sequence homology between these two proteins, and they belong to evolutionary unrelated virus families. Comparison of the crystal structures of the DNA-complexed forms suggested a quite different DNA binding recognition mode (10).
The portion of the EBNA1 DNA binding and dimerization domain that resembles the E2 DNA binding domain from papillomavirus is referred to as the core domain (residues 504604). This includes the dimeric barrel interface and two
-helices, one of which is the recognition helix. A noticeable difference in the core domain is the length and position of a loop connecting
-strands 2 and 3. In EBNA1, this loop contains nine extra residues, five of which are proline residues, and is therefore termed the proline loop. Although the helix corresponding to the recognition helix in E2 domain is present in the EBNA1, it is located 6 Å away from the DNA. However, mutagenesis of this helix and the fact that the core alone binds DNA suggest a role of this helix in DNA binding (11). The EBNA1 flanking domain (Fig. 1) is composed of an extended chain (461476) that tunnels through the minor groove of the
-helix (477489), and a connector (490503) linking this domain with the core domain. An acidic tail is located at the C terminus of the protein (608641). Neither acidic nor flanking domains or equivalent domains are present in papillomavirus E2 C-terminal domains. The Epstein-Barr genome includes repetitive contiguous EBNA1 binding elements, and assembly at these sites is cooperative, where the DNA binding domain alone is sufficient for cooperative binding (12). Changes in the DNA, including bending, take place upon binding (10), which result in permanganate sensitivity (13), and DNA unwinding is predicted to accompany assembly of EBNA1 on two adjacent sites (10).
|
-barrel domains described so far that incidentally bind at or near DNA replication origins. In the present work, we analyzed the DNA binding mechanism of EBNA1 to a single specific DNA site in solution, a model for eukaryotic replication origins, using various spectroscopic, biophysical, and calorimetric techniques. We integrated structural thermodynamic and kinetic data and compared them with the E2 DNA binding domain. | EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
DNA SynthesisDouble-stranded 20-bp oligonucleotides containing one EBNA1 recognition sequence (Site 1 of the dyad symmetry element of oriP in the Epstein-Barr B95-8 strain genome) were prepared as follows: single-stranded oligonucleotides were purchased, high-performance liquid chromatography-purified, from Integrated DNA Technologies (Coralville, IA): Site 1A, 5'-CGGGAAGCATATGCTACCCG-3' (recognition sequence is italicized); Site 1B, 5'-6-carboxyfluorescein (FAM) is the complementary strand with a fluorescein molecule attached to the 5'-end via a 6-carbon linker. Singlestranded oligonucleotide concentration was calculated using the molar extinction coefficient at 260 nm obtained from the nucleotide composition. Annealing was performed by mixing equal amounts of the oligonucleotides in 10 mM Bis-Tris HCl buffer, pH 7.0, and 100 mM NaCl, further incubating the mixture for 5 min at 95 °C, and slowly cooling to 25 °C for 16 h. This yielded a double-stranded oligonucleotide termed Site 1 5'-FAM, and no detectable single-stranded oligonucleotide was present as judged by PAGE (not shown). A similar procedure was followed to anneal the Tesi-1-5'-FAM 5'-ATCCCAGCCAGTGTCGAGAG-3'. This oligonucleotide corresponds to the randomized Site 1 sequence keeping the same base composition. The double-stranded E2 oligonucleotide, termed Site 35-5'-FAM and used as nonspecific DNA, was 5'-GTAACCGAAATCGGTTGA-3'.
Protein Expression and PurificationRecombinant fulllength EBNA1 was expressed and purified from baculovirusinfected insect cells as previously described (17). The procedure is much more complex, the yields are lower, and this protein cannot be expressed in bacteria. The recombinant C-terminal domain of EBNA1 (EBNA1452641) was expressed in Escherichia coli BL21(DE3) pLys S, using a modified protocol from the one previously described (18, 19). Bacterial cells expressing EBNA1452641 were grown in a 1.5-liter fermenter (New Brunswick Scientific Co., Inc, Edison, NJ) in Terrific Broth culture medium at 37 °C, to an absorbance of 5 at 600 nm. EBNA1452641 expression was induced with 0.5 mM isopropyl
-D-thiogalactopyranoside, followed by the addition of 150 µg/ml rifampicin 3 h later. The cells were incubated overnight, harvested by centrifugation, and resuspended in 0.1 volume of the lysis buffer (100 mM Tris HCl, pH 6.8, 600 mM NaCl, 1.0 mM EDTA, 1.0 mM phenylmethylsulfonyl fluoride, and 10 mM 2-mercaptoethanol). Cells were frozen at 70 °C and then thawed at room temperature, and subsequently lysed by sonication at 0 °C twice. The lysate was clarified by centrifugation at 20,000 x g for 20 min. The supernatant was then placed in a 75 °C water bath until the temperature of the protein solution reached the temperature of the bath. Heating was followed by a 10-min incubation on ice. The supernatant was clarified by centrifugation at 20,000 x g for 40 min and then loaded onto a heparin HyperD (BioSepra, Villeneuve la Garenne, France) affinity column equilibrated with buffer A (50 mM Tris-HCl, pH 7.2, 0.2 M NaCl, and 5 mM 2-mercaptoethanol), washed with 5 column volumes of the same buffer, and eluted with 0.151.0 M NaCl linear gradient. The fractions that were over 90% pure were pooled, dialyzed against buffer A without 2-mercaptoethanol, and digested with thrombin (1.5 units/mg of protein) at 37 °C for 7 h to remove the His6 tag. Digestion was stopped with 0.2 mM phenylmethylsulfonyl fluoride, and the protein was loaded onto a Mono Q ion exchange column (Amersham Biosciences), equilibrated with buffer A, and washed with 5 column volumes of the same buffer. The protein eluted in a 15-ml linear gradient from 0.2 to 1.0 M NaCl in buffer A. The fractions that were greater than 95% pure as judged by SDS-PAGE were pooled, concentrated using Centriprep-10 (Amicon, Bedford, MA), and loaded onto a Superdex 75 gel-filtration column (Amersham Biosciences). This procedure yielded around 20 mg/liter of >98% pure EBNA1452641. The purified protein was dialyzed against 25 mM Tris HCl, pH 7.2, 500 mM NaCl, 0.2 mM EDTA, 10% glycerol, and 5 mM 2-mercaptoethanol, and stored at
150 µM fractions at 70 °C after snap freezing in liquid nitrogen. Protein concentration was determined using an extinction coefficient of 4.19 x 104 M1 cm1.
DNA BindingFluorescence measurements were recorded in an Aminco Bowman series 2 luminescence spectrometer assembled in "L" geometry. For fluorescein anisotropy measurements excitation was set to 495 nm with a 4 nm slit, and emission was recorded at 520 nm. When fluorescein concentration was lower than 20 nM, the excitation slit was set to 8 nm. The temperature was kept constant at 25 ± 0.1 °C through all experiments. All titrations were performed adding small amounts of a concentrated solution of the variable ligand to fixed amounts of a concentrated solution of the other and allowed to equilibrate for 2 min. In all cases, maximal dilution was 20%, and the data were corrected accordingly.
Dissociation constants for the EBNA1·DNA complex were performed in 25 mM Bis-Tris-HCl (pH 7.0), 200 mM NaCl, and 1mM dithiothreitol by measuring the steady-state fluorescence anisotropy of the 5'-FAM·DNA as a function of added EBNA1. Data were fitted using nonlinear least squares to Equation 1,
![]() |
F is the difference in the signal between the EBNA1·DNA complex and free DNA; [DNA] and [EBNA1] are the oligonucleotide and protein concentrations, respectively; and KD is the dissociation constant for the interaction. No computational corrections for emission intensity were required, because the quantum yield did not change significantly upon binding. Electrophoretic Mobility Shift AssayStandard equilibrium EBNA1·Site 1 binding was carried out in identical conditions to equilibrium binding by fluorescence spectroscopy (25 mM BisTris-HCl pH 7.0, 0.2 M NaCl, and 1 mM dithiothreitol) at 2 µM Site 1 5'-FAM and different molar ratios of EBNA1452641, followed by 60-min incubation at room temperature in a final volume of 50 µl. For the EBNA1·Site 1 dissociation experiment, reaction mixtures in the same buffer containing 1 µM EBNA1452641·Site 1 5'-FAM complex and large excess of unmodified Site 1 were incubated 60 min at room temperature in a final volume of 15 µl. Finally, for EBNA1·Site 1 off-rate kinetics, the procedure was as follows: reaction mixtures in the same buffer containing a 5 µM EBNA1452641·Site 1 5'-FAM complex were incubated 60 min at room temperature and then a 30-fold excess of unmodified Site 1 was added to a final volume of 100 µl to displace the Site 1 5'-FAM from the complex. Mixtures were loaded continuously into running 10% non-denaturing polyacrylamide gels containing TBE (0.1 M Tris-HCl, pH 8.0, 0.15 M sodium borate, 4 mM EDTA) 0.5x. The gels were resolved at 4 V/cm, 48 °C for 28 h. Fluorescein bands were detected by UV transilumination. After visualization and documentation, the gels were silver stained.
Circular DichroismCD spectra were monitored in the near UV region using Jasco J-810 equipment. Ten scans were averaged for each measurement at 25.0 ± 0.1 °C controlled by a peltier, and contribution of the protein was subtracted. It is expected that in the 250 to 320 nm range, the ellipticity of the protein does not change substantially (i.e. the contribution is linear), and thus we can assume we are evaluating mainly changes in DNA conformation (20).
Isothermal Titration CalorimetryAll experiments were conducted using a VP ITC (MicroCal Inc., Northampton, MA). Titrations were performed as described elsewhere (21, 22). In a typical experiment, 8-µl injections of EBNA1452641 were made into DNA solution in the cell. The heat of dilution of the protein into the buffer was determined in separate experiments and subtracted from the titration prior to data analysis. The data were analyzed using the ORIGIN software supplied with the calorimeter. The dissociation binding constants in the temperature range of 293308 K were too high to be determined from the ITC results, because of the high concentrations of protein that were required to detect the heats of binding and the tight binding involved.
Stopped-flow KineticsAll fluorescence stopped-flow kinetic experiments were performed using an SX.18MV stopped-flow apparatus (Applied Photophysics, Leatherhead, UK). Concentrations reported are those resulting from mixing equal volumes of each syringe at 25 ± 0.1 °C, unless otherwise stated. The reactions were monitored using the fluorescence of the 5'-FAM moiety of the modified oligonucleotides or the intrinsic Trp fluorescence of the protein. Excitation was set to 490 nm or 280 nm and the total fluorescence recorded through cut-off filters (Schott, PA). Five to ten kinetics traces were collected and averaged for each concentration point. The data were analyzed to extract the rates and amplitudes using non-linear least squares fitting software provided by the manufacturer, using single or double exponential equations. The errors for the constants were between 2 and 10%. The Double-jump association-dissociation experiments were performed mixing 0.5:0.5 volumes of protein and 5'-FAM·DNA, and, after the "delay time," the reaction sample was subsequently mixed with one volume of nonfluoresceinated DNA, and the fluorescent decay was measured.
| RESULTS |
|---|
|
|
|---|
|
As a probe for the binding from the DNA end, we used near-UV CD, which allows discrimination of changes exclusively taking place on the DNA. Upon addition of EBNA1452641 to an unmodified Site 1 duplex, there was a change in the spectra, indicative of a substantial change on its conformation upon binding of the domain (Fig. 2C). An increase in the bands at
280 nm could be indicative of DNA unwinding, but because the contributions in this region are so complex, it is virtually impossible to assign a particular structure (20). The inset of the figure shows that the stoichiometry was also 1:1, as judged by the ellipticity change in the DNA followed at two different wavelengths (Fig. 2C, inset).
To determine the dissociation constant, we carried out binding isotherms under near dissociation conditions, following fluorescein anisotropy at different DNA and protein concentrations (Fig. 3A). The data were fitted to a quadratic binding equation, and the residuals are indicated in the inset. The KD was determined to be 0.5 ± 0.1 nM, in excellent agreement with previous EMSA results (11). Using the same procedure, the binding of a nonspecific sequence with the same base content yield as the specific E2 site showed a KD > 10 µM (Fig. 3, inset), and the discrimination capacity (KDnon-specific/KDspecific) is at least 20,000, in the same range of what was observed for HPV16 E2C (14).
As a test for the biological relevance of studying the DNA binding domain rather than the difficult to express full-length form, we determined the dissociation constant of the baculovirus-expressed EBNA1 (see "Experimental Procedures") for comparison. We carried out the comparison in 400 mM NaCl in the same buffer, because it weakens the extremely tight interaction, and allows an accurate determination with the best signal-to-noise ratio possible. In these conditions, the KD of EBNA1 was 1.0 ± 0.1 nM, compared with 1.2 ± 0.2 nM for the EBNA1452641 domain to the same DNA site, which validates our analysis.
Thermodynamics of EBNA1452641-DNA BindingTo evaluate the thermodynamic parameters governing the interaction of EBNA1452641 with its specific DNA site, we carried out ITC experiments. Upon addition of EBNA1452641 to a solution containing the EBNA1 Site 1 duplex, a strong enthalpic change was observed, where the raw data indicated a tight and saturable interaction, as expected (Fig. 4A). The data were fitted to a standard binding equation, showing a 1:1 stoichiometry, and the dissociation could not be determined accurately, due to the limits imposed by the technique. The KD used was that determined from fluorescence titration experiments: 0.5 x 109 ± 0.1 x 109 M (free energy for the binding process
G = 12.7 ± 0.1 kcal mol1); the
Hobs was 11.9 ± 0.14 kcal mol1 (Fig. 4A). Using these data, the entropic change,
S was calculated to be only +2.7 ± 0.9 cal mol1 K1, indicating that the interaction is, to a great extent, enthalpically driven.
|
Cp), which is related to the amount of hydration of the polar or apolar molecular surfaces (24). ITC experiments were carried out at different temperatures, and the
Hobs was plotted against the temperature, and the slope corresponds to the
Cp of the interaction. The formation of the EBNA1·DNA complex was accompanied by a
Cp of 1.22 ± 0.08 kcal mol1 K1.
Kinetic Dissociation of EBNA1452641·DNATo determine the overall binding mechanism, we started analyzing the kinetic dissociation of EBNA1452641 to DNA. For this, we incubated an EBNA1452641·5'-FAM Site 1 complex with an excess of unmodified EBNA1 Site 1 DNA duplex, and determined the anisotropy change after 4 h. A 25-fold excess of unmodified DNA ensures dissociation of the complex as the anisotropy change indicated (not shown). The formation of the complex was very fast as the change in anisotropy of the 5'-FAM Site 1 was immediately increased upon addition of EBNA1452641 protein (not shown). After a 30-min stabilization period, addition of a 25-fold excess of unmodified Site 1 DNA decreased the anisotropy, as expected for the dissociation (Fig. 5A). The anisotropy change that accompanies dissociation was fitted to two exponential phases, and the residuals are indicated (Fig. 5A, bottom panel). The phases contribute each with 50% of the amplitude, and the first order rates were 2.8 x 104 s1 for kr3, and 2.7 x 103 s1 for kr2, corresponding to t
values of 41 and 4.3 min, respectively. A similar experiment was carried out in the stopped-flow equipment monitoring fluorescence intensity. The advantage of this technique is that faster reactions can be observed, but there is a limitation for slow reactions that can be observed in a standard fluorometer. Fig. 5B shows the data in the 0 to 1000-s range where a fast phase corresponding to 20% of the amplitude and a kFAM1 rate of 0.011 s1 was observed, and a slower phase with a rate of 0.001 s1 was observed, which we assigned by coincidence to the kr2 from anisotropy (Fig. 5A). We confirmed the complete dissociation of the DNA by an EMSA experiment. The EBNA1452641·5'-FAM Site 1 complex (Fig. 5B, inset, lane 3) was displaced into the free 5'-FAM Site 1 by excess unmodified Site 1 DNA (Fig. 5B, inset, lanes 48).
|
|
of
100 min, in very good agreement with the anisotropy experiment (41 min), considering the large difference in the accuracy of the experiments (Fig. 6B).
|
Two-step Association-Rearrangement MechanismThe association kinetics was investigated by a pseudo first order experiment, where a fixed concentration of 5'-FAM Site 1 DNA oligonucleotide was mixed with increasing concentrations of EBNA1452641 in a stopped-flow fluorometer. There was a fast fluorescence increase involving at least two phases, as the data fit to two exponential components (Fig. 7A, and insets). The observed pseudo first order rates were plotted against protein concentration, and the data are shown in Fig. 7B. The major phase (k1, Fig. 7B), accounting for 85% of the amplitude, is concentration-dependent and thus corresponds to the association event, with a kFAMon1 of 1.6 x 108 M1 s1, in the range of a diffusion controlled reaction. The minor phase, accounting for 15% of the amplitude, did not depend on the concentration, indicating that it corresponds to a conformational rearrangement, most likely after the collision event, with a first order rate of 2.0 s1 (k2, Fig. 7B). The extrapolation to zero protein provides an estimate of the koff, but it was only possible to confirm it as
0.5 s1.
We carried out a similar stopped-flow experiment in pseudofirst order conditions but following the change in the intrinsic tryptophan fluorescence of EBNA1452641 upon binding to DNA. In this experiment two binding events were observed: a major phase (k1, Fig. 7C) accounting for 85% of the amplitude was also observed with a kWon1 of 1.0 x 108 M1 s1, confirming the association rate with a different probe (Fig. 7C). A minor phase (k3, Fig. 7C) corresponding to 15% of the amplitude and a second-order rate kWon2 of 0.25 x 108 M1 s1 indicated two parallel binding events. The sensitivity of tryptophan fluorescence was comparatively much lower than the 5'-FAM fluorescence, so we were not able to measure the minor phase as accurately. In any case, the minor phase was slightly but significantly lower.
Double-jump ExperimentsA double-jump association-dissociation experiment consists of mixing EBNA1452641 protein with 5'-FAM Site 1 DNA, allowing the binding reaction to proceed, and at different "delay" times ranging from 0.01 to 100 s, it dissociates the complex with an excess of unmodified Site 1 DNA and monitors the fluorescence change accompanying the dissociation. When the complex was formed in excess of protein (800 nM protein: 100 nM 5'-FAM DNA), and displaced with excess of unlabeled DNA, a single phase was observed in a 200-s time frame with a rate of 0.014 s1 (kDJ1, Fig. 8C), coincident with the faster rate observed for dissociation of the complex followed by stopped flow (kFAM1, Fig. 5C). Under these experimental conditions the kobs for complex formation was
100 s1 (see Fig. 7C), which was faster than the dead-time of the double-jump experiment (10-ms delay time). Therefore, its disappearing amplitude was not observable. The progressive change of the amplitudes with the delay time is expected to be equivalent to the formation of the species in the "forward" direction. The increase in the single observed amplitude fitted to a single exponential (Fig. 8A) showed a rate of 0.8 s1 in agreement with the rate of the slow phase from pseudo-first order (2 s1, Fig. 7B), which we now identify as the rearrangement of the encounter complex to yield the consolidated complex.
|
|
41 min, Figs. 5A and 6A) was too slow to fit in the experimental window of a stopped-flow experiment. The rates obtained from the amplitude change with delay time were 1.8 and 1.4 s1, k3 and k4, respectively (Fig. 8B). These rates were identical to the slow unimolecular rate observed in pseudo-first order, so they could not be discriminated in the forward direction but were only dissected in a double-jump experiment from their reverse reactions. The existence of both amplitudes from the earliest delay times indicates that they correspond to parallel rearrangement channels, and the fact that the percent amplitudes changed with time strongly suggests that they are connected, i.e. they can exchange. As the reaction progressed, the amplitudes stabilized at 65 and 35%, respectively (Fig. 8B, inset); however, we do not consider there should be a linear relationship between the fluorescence change and the actual percentage of the population.
Finally, under the experimental conditions for the association silent phase of the association-dissociation double-jump experiment (100 nM, 1:1 EBNA1·DNA), the association kobs was
20 s1 (t
35 ms). Considering a 10-ms dead-time for the delay mixing, one would have expected the amplitude corresponding to the encountered complex to disappear progressively, thus yielding the subsequent rearrangement described above. However, this phase did not yield any fluorescence change and could not, therefore, be followed due to experimental restrictions.
| DISCUSSION |
|---|
|
|
|---|
The binding affinity of EBNA1452641 to DNA in solution was 0.5 nM, which is in excellent agreement with that determined from EMSA, particularly considering the differences in the experimental conditions (11). The discrimination capacity (KDnon-specific/KDspecific) was determined to be
20,000, which is comparable to that of HPV16 E2C·DNA interaction with one of its specific sites: 10,000 (14).
Thermodynamic analysis of the interaction indicates that it is entirely enthalpically driven, which appears evident from the comparison of
G (measured by fluorescence spectroscopy) and
H, and a negligible
S. This also suggests that most of the enthalpic change comes from binding as opposed to local folding events, suggesting that the molecules behave as rigid bodies compared with other systems (26). The
Cp estimated from established algorithms using burial of surface area (
ASA) values (27) is 0.54 kcal mol1 K1, compared with 1.22 kcal mol1 K1 for that determined experimentally in this work. This type of discrepancy is often found in protein-DNA interactions, because it is not as predictable as in protein folding reactions or in small molecules (28). On the other hand, the measured
Cp for the structurally homologous and functionally related E2 DNA binding domain from human papillomavirus is 0.37 kcal mol1 K1 (29), which shows little discrepancy with the predicted value, 0.46 kcal mol1 K1. The EBNA1452641·DNA complex buries 5,850 Å2 against 3,544 Å2 for the E2C·DNA counterpart. In addition, whereas 50% of the buried surface area in EBNA1452641·DNA complex formation is polar; this figure is 64% in the E2C·DNA complex, suggesting a less polar interface for EBNA1·DNA.
Although it is generally accepted that a large fraction of the
Cp is accounted by burial of surface area, there are several possible sources of
Cp, and therefore, for the discrepancy between the "surface area only" empirical approach and the experimental data (30). The most obvious is a local or global folding process coupled to binding (24). The EBNA1 domain is very stable to chemical denaturation (11),3 and this evidence limits a significant contribution by a coupled local, let alone global, folding process. Other possible sources are protonation and ion binding events (31), but the EBNA1·DNA complex presents little effect of pH and salt concentration on the binding.4 A recently described alternative explanation is the role of the network of water molecules at the interface, related in part to changes in soft vibrational modes of the trapped interface water molecules (30). With the present limited information, we are in favor of the latter interpretation, i.e. local changes in the cooperative network of interactions, possibly involving structural interfacial water molecules. This view agrees very well with a less polar interface in the EBNA1·DNA complex. Perhaps more polar interfaces provide better estimations for the algorithm.
The association reaction displays at least two major phases, one of them corresponding to an encounter complex taking place near the diffusion limit, comparable to the association rate of the HPV16 E2C domain and other systems (15, 32, 33). This species undergoes a conformational rearrangement to yield the final consolidated complex (15). An identical association phase is observed when monitoring changes in intrinsic tryptophan fluorescence, which is carried out in excess of DNA instead of protein, confirming there is one major association phase. Our interpretation is that the difference in the observed phases in either excess of protein or DNA comes from the presence of not unexpected protein conformers in solution, where in excess of protein, the fastest reacting protein conformer reacts preferentially, and a single phase was observed (Fig. 7B). On the other hand, in excess of DNA, there were plenty of these molecules to react with all EBNA1452641 species (at least two populations of conformers), and two phases were observed. The fact that the on-rates of these phases are separable supports conformationally distinct populations (Fig. 7C). The microheterogeneity of the protein may be potentially more extensive, but only two rates are distinguishable. In support for this, we observed heterogeneity in the EBNA1·DNA complex in EMSA experiments (Fig. 6A). Based on the unusual positioning of the major "binding helix" in the crystal structure of EBNA1·DNA and the presence of a flanking domain that travels through the DNA helix making base specific contacts, consecutive binding and rearrangement steps are the most likely explanation for this reaction pathway (11).
Three phases were observed for the dissociation reaction, in different time ranges, using different spectroscopic probes, experimental set-ups, and even a kinetic EMSA experiment. The best estimation of these rates are kDJ1 = 0.014 s1, kDJ2 = 0.0014 s1, and kr3 = 0.00028 s1, where the latter (t
41 min) corresponds to the rate-limiting dissociation, because it agrees with that obtained in the EMSA experiment, a direct physical evidence of a slow dissociation, considering the large difference in experimental conditions. However, although the EMSA experiment cannot yield a large number of points to fit complex reactions, it does provide an accurate indication that, at the experimental dead time (
1 min), 23% of the molecules have already dissociated (Fig. 6B). This is consistent with the fastest rate (kDJ1, 0.014 s1).
When the double-jump experiments were carried out with excess of protein in the association stage (first mixing), a single amplitude in the forward direction (Fig. 8A) corresponded to the rearrangement observed in the pseudo-first order experiment (
2 s1, Fig. 7B). An additional phase in the forward direction appeared in the analysis of the double-jump amplitudes when the complex was formed in equimolar concentrations (Fig. 8B), where both protein conformers could bind in parallel to DNA molecules. This additional forward phase, judging by the rate (k4, 1.4 s1), would not be separable from the observed rearrangement in the forward direction even if it were accompanied by a strong fluorescence signal (k2, 2 s1). Because it is not formed at the expense of the disappearance of k2, we propose that this corresponds to the post-collision rearrangement of the parallel pathway, which not surprisingly appears similar to the single rearrangement directly observed in the forward reaction. Fig. 9 provides the simplest model for association and dissociation events that is compatible with the results presented here but in no way rules out other interpretations (Fig. 9, A and B).
|
kWon1 = kon1, 1.6 x 108 s1 M1). However, the existence of a fast dissociation route in these conditions (kFAM1 = 0.011 s1) coincident with the early fast dissociating species from the double-jump experiment (kDJ1 = 0.014 s1 in excess of either protein or DNA) strongly suggests the presence of a parallel fast route, such as the one we described for the HPV E2·DNA complex (15). The collision events of direct and rearrangement mediated pathways are diffusion-controlled reactions and likely to display similar values. Therefore the kDJ1/kFAMon1 ratio yields a value of 0.090.14 nM for the KD, in excellent agreement with that determined from equilibrium experiments (0.5 nM, Fig. 3), supporting the existence of this direct fast route. We can hypothesize that the fast reacting EBNA1452641 conformer involved in the direct route is that containing an ordered "arm" domain (Fig. 9C, left). The route progressing through rearrangement will involve a population of EBNA1452641 conformers that bind through the core domain only (Fig. 9C, right) and require rearrangement at the interface to yield a consolidated complex. Overall, we pay more attention to the protein conformers, because the conformational heterogeneity of the DNA is likely to be greatly diminished in the context of the genome, as opposed to a small duplex, and is expected to be in fast exchange. Because the encounter need not discriminate base sequence, the complex formed in the two-state route, with a faster off-rate, may therefore need to undergo dissociation to re-bind to the specific target sequence while the intermediate encounter complex may slide along the DNA and find the correct sequence. A possible mechanism, based on the crystal structure and mutagenesis, involves the interaction of the core domain followed by the interaction of the flanking domain with the DNA site (11). We have shown that such a mechanism is consistent with the results we present now, but further mutagenesis of the EBNA1452641 domain, in particular at the major core helices and the DNA-contacting arm, will be required to confirm the species involved in each step. In addition, NMR techniques in solution will be required for a detailed structural and dynamic investigation of the protein·DNA complex.
| FOOTNOTES |
|---|
1 A Career Investigator from Consejo Nacional de Investigaciones Científicasy Técnicas. To whom correspondence should be addressed. Tel.: 54-0115238-7500; Fax: 54-011-5238-7501; E-mail: gpratgay{at}leloir.org.ar.
2 The abbreviations used are: OBP, origin binding protein; EBNA1, Epstein-Barr nuclear antigen 1; EBV, Epstein-Barr virus; Bis-Tris, 2-[bis(2-hydroxyethyl)amino]-2-(hydroxymethyl)propane-1,3-diol; EMSA, electrophoretic mobility shift assay; ITC, isothermal titration calorimetry; HPV, human papillomavirus. ![]()
3 C. Oddo, E. Freire, L. Frappier, and G. de Prat-Gay, unpublished results. ![]()
4 C. Oddo and G. de Prat-Gay, unpublished results. ![]()
| REFERENCES |
|---|
|
|
|---|