Published online 1 August 2005
Article |
Determination of proteinDNA binding constants and specificities from statistical analyses of single molecules: MutSDNA interactions
1Department of Chemistry, University of North Carolina at Chapel Hill Chapel Hill, NC 27599-3290, USA 2Curriculum in Applied and Materials Sciences, University of North Carolina at Chapel Hill Chapel Hill, NC 27599-3290, USA 3Genetics and Biochemistry Branch, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health Bethesda, MD 20892, USA
*To whom correspondence should be addressed. Tel: +1 919 962 6370; Fax: +1 919 966 3675; Email: derie{at}unc.edu
Received May 4, 2005. Revised June 13, 2005. Accepted June 29, 2005.
| ABSTRACT |
|---|
|
|
|---|
Atomic force microscopy (AFM) is a powerful technique for examining the conformations of proteinDNA complexes and determining the stoichiometries and affinities of proteinprotein complexes. We extend the capabilities of AFM to the determination of proteinDNA binding constants and specificities. The distribution of positions of the protein on the DNA fragments provides a direct measure of specificity and requires no knowledge of the absolute binding constants. The fractional occupancies of the protein at a given position in conjunction with the protein and DNA concentrations permit the determination of the absolute binding constants. We present the theoretical basis for this analysis and demonstrate its utility by characterizing the interaction of MutS with DNA fragments containing either no mismatch or a single mismatch. We show that MutS has significantly higher specificities for mismatches than was previously suggested from bulk studies and that the apparent low specificities are the result of high affinity binding to DNA ends. These results resolve the puzzle of the apparent low binding specificity of MutS with the expected high repair specificities. In conclusion, from a single set of AFM experiments, it is possible to determine the binding affinity, specificity and stoichiometry, as well as the conformational properties of the proteinDNA complexes.
| INTRODUCTION |
|---|
|
|
|---|
Understanding proteinDNA interactions is fundamentally important for dissecting the molecular mechanisms underlying many biological processes. Association constants and specificities of protein binding to DNA are the primary thermodynamic properties for understanding proteinDNA interactions. Many methods, such as electrophoretic mobility shift assays (EMSA), filter binding assays, surface plasmon resonance (SPR) and calorimetric assays are used to investigate the thermodynamic equilibrium constants of proteinDNA interactions (15). Although these methods are very powerful, they all have two significant limitations. First, all are bulk measurements; therefore, the observed affinities are the weighted sum of all interactions occurring between the protein and the DNA (Figure 1a) (6). For example, if a protein has a significant binding affinity for the ends of the DNA, the apparent binding constant may represent this preference, especially for nonspecific binding. Second, in all of these assays, the measurement of binding is indirect, and it is generally assumed that the signal, such as heat in calorimetry or refractive index in SPR, is linearly proportional to the binding (Figure 1a). While this situation is often the case, there are many cases when this assumption is not valid (2).
|
A single molecule method to determine proteinDNA binding constants can overcome these limitations. Accordingly, we have developed a single molecule method using atomic force microscopy (AFM) to determine proteinDNA binding constants and specificities directly at the level of DNA-binding sites (DNAi, Figure 1b). Using AFM, it is possible not only to determine the extent of binding of a protein to a DNA fragment but also to determine where on the DNA the protein is bound, (i.e. fractional occupancies at the specific site, the ends, or nonspecific sites, Figure 1b). In the sections that follow, we present the theory that demonstrates how binding affinities, as well as binding specificities, can be determined from the analysis of AFM images of proteinDNA complexes. We then demonstrate its utility and accuracy by analyzing complexes of the mismatch repair protein MutS with DNA and comparing our results with those of bulk measurements.
The MutS family of proteins is highly conserved and the sole known factor for the recognition of DNA mismatches and small insertion/deletion loops in the DNA mismatch repair (MMR) pathway (7,8). MutSDNA interactions are crucial in the regulation of the initiation of MMR (9,10). Studies on Escherichia coli MutS and eukaryotic MutS homologs using traditional bulk techniques show that the binding specificities to various mismatches are very low (
30 or less) (11,12). This relatively low binding specificity to mismatches versus much higher expected MMR specificity is one of central puzzles in MMR (13,14). Interestingly, EMSA studies of Taq MutS binding to the single T-bulge, however, suggest a much higher binding specificity (>1000), although the specificities for other mismatches are similarly low (11). In this paper, we present a detailed analysis of MutSDNA interactions using AFM. Our results indicate that the binding specificities of MutS are greatly underestimated in the previous studies and suggest that this underestimation is due, in part, to a high affinity of MutS to DNA ends.
| THEORY |
|---|
|
|
|---|
Site-specific binding constant
From the lattice binding model of proteinDNA interactions (5,15), a protein interacts at another DNA-binding site whenever it moves 1 bp or more away from the current binding position. In other words, the number of binding sites (N) on a linear DNA fragment is Nbp Nbp,P + 1, where Nbp is the length of the fragment in base pairs and Nbp,P is the number of base pairs occupied by protein (5,15). (If end binding is a different mode than nonspecific binding, then the total number binding sites on the DNA fragment will increase by 2, i.e. N = Nbp Nbp,P + 3.) The number of binding sites on a circular DNA fragment is equal to the number of base pairs (i.e. N = Nbp). Assuming that all sites are independent, the binding constant, Ki, to a given site i is:
![]() | (1) |
), and [P] and [D] are the total concentrations of protein and DNA fragments, respectively. The right hand side of Equation 1 is derived by expressing the concentrations in terms of occupancies. Specifically, [ProteinDNAi] = Oi x [DNAi]Total, [Protein] = [P] [D] x OFragment (i.e. the free protein concentration equals the total protein concentration minus the concentration of protein bound to the DNA), and
![]() |
The term
is included because the protein-binding site size (Nbp,P) is >1 bp, and therefore, the concentration of free site i will depend not only on those proteins bound at i, but also on those bound at sites that are within Nbp,P 1 bp of i (assuming that Nbp,P is independent of position).
Specific and nonspecific binding constants
Normally, only one or a few binding sites on a DNA fragment are specific and all other binding sites are nonspecific with similar protein binding affinities. Consequently, the binding constant for the specific site (KSP) and the average binding constant for nonspecific binding sites (KNSP) are of interest. In addition, for a linear DNA fragment, proteins may have significant binding affinities for DNA ends (KE). A linear fragment containing NSP specific sites and a total of N binding sites has N NSP 2 nonspecific sites and two end binding sites. Using Equation 1, the binding constant for a given specific site (KSP), and the average binding constants for a DNA end (KE) and for a nonspecific site (KNSP) can be expressed in terms of fractional occupancies of the sites. Specifically, for a linear DNA fragment containing NSP specific sites that are separated from each other and from the DNA ends by at least Nbp,P nonspecific sites:
![]() | (2) |
![]() | (3) |
![]() | (4) |
and
are included in Equation 4 to account for the occlusion of the nonspecific sites by protein binding at or near the specific sites and the DNA ends, respectively.
Under conditions of low occupancy [ONSP << 1 and (OSP + OE)/NNSP << 1, i.e. DNA is sufficiently long], the binding constants for a linear fragment containing a single specific site can be approximated within experimental error as:
![]() | (5) |
![]() | (6) |
![]() | (7) |
Binding specificity
Specificity is the relative affinity of a protein binding to a specific site versus a nonspecific site, i.e. S = KSP/KNSP. It is also the relative occurrence probability of protein binding at different sites on the same DNA fragment. Consequently, it is not necessary to know the absolute binding constants to determine the specificity, and AFM provides a straightforward method for estimating binding specificity (1618). If a piece of DNA contains no specific sites, a uniform distribution of protein is seen on the DNA; however, if the DNA contains a single specific site, a Gaussian distribution of protein will be observed centered at the specific site. The Gaussian distribution is a result of the error in the determination of the position of the protein. All of the proteins in the Gaussian distribution (Asp, Figure 1c) result from the presence of the specific site. Accordingly, the specificity, S, can be determined by integrating the areas under the Gaussian curve of the protein distribution (Asp and Ansp, Figure 1c). The binding specificity can be given by the following mathematically different but physically consistent forms (see derivation in Supplementary Material):
![]() | (8) |
Theoretically, the distribution of protein-binding positions is definite if the number of binding sites (N) and the binding specificity (S) are given. Consequently, another way to determine specificity or to compare specificities of two different binding sites is given by the following correlation of protein binding on any two DNA fragments that contain up to one specific site in each fragment (see derivation in Supplementary Material):
![]() | (9) |
), determining the specificity in the second fragment simply requires determining the probability of protein binding at nonspecific sites (
), because
. Consequently, this correlation provides a simpler way to determine the binding specificity as long as the binding probability on nonspecific sites (Pmin) can be determined. This method is especially useful under the conditions where determining the integration of Asp is difficult, e.g. if the position distribution is not fit well by a Gaussian distribution.
Error of measurements
The error induced by the uncertainty in Oi and OFragment can be determined from multiple AFM experiments and the relative errors of the binding constant and specificity can be calculated by applying the theory of error propagation (19) for Equations 5 7 and S = KSP/KNSP. The relative errors are:
![]() | (10) |
![]() | (11) |
![]() | (12) |
represents the standard error of physical variables and the upper bar represents the mean of physical variables in multiple AFM experiments. Inspection of these relations shows that the relative errors of binding constants and specificities are close to the relative errors of the DNA fractional occupancy, given conditions of low occupancy (OSP, OE < 0.5 and ONSP << 1) and
(i.e. at least 50% proteins are unbound), which are fulfilled in this study of MutS.
Relationship between site-specific and fragment-specific binding constants and specificities
The macroscopic DNA-binding constant determined by bulk solution measurements is a measure of the protein-binding affinity to the entire DNA fragment (Figure 1a). Consequently, to compare results from AFM, which provides a measure of binding affinity to individual sites (Figure 1b), it is necessary to calculate the binding constant to the entire DNA fragment from AFM microscopic constants. For a linear DNA fragment, it can be shown, based on the lattice binding model of proteinDNA interactions, that the binding constants to DNA fragments with N binding sites that contain either a single specific site (KSP,Macro) or no specific sites (KNSP,Macro) are:
![]() | (13) |
![]() | (14) |
Similarly, the bulk specificities are determined from the ratio of binding affinities to two DNA fragments: one with and one without the specific site (Figure 1a); whereas, AFM yields a direct measure of the relative affinities for different sites (Equation 8). The macroscopic specificity can be calculated from the ratio of KSP,Macro/KNSP,Macro (Equations 13 and 14). Inspection of Equations 13 and 14 shows that both specificities and binding constants determined by bulk measurements may be skewed if specificities are low such that nonspecific binding makes a significant contribution to the overall binding or if there is significant affinity for DNA ends. On the other hand, for proteins that have high specificities and low end binding affinities, the macroscopic and microscopic binding constants and specificities should be similar. EcoRI, which has a very high specificity (
106), provides an example of the latter case. Specifically, a recent single-molecule study, which used an optical trap to unwind a DNA double helix with EcoRI bound, estimated the binding constant of EcoRI to its specific site and found it to be very similar to the bulk measurements (20). Finally, in comparing bulk measurements with AFM, it is important to keep in mind the limitations of the methods for determining different types of binding. For example, proteinDNA interactions that are very dynamic (rapid on and off rates) are often not detected in EMSA or filter binding assays (21). In contrast, with AFM, such interactions should be detected approximately as well as those with slow dissociation rates, because deposition of both proteins and DNA on mica is diffusion-limited, and both bind to mica irreversibly over the time scale of the deposition (seconds to minutes) (2225). Consequently, AFM detection of proteinDNA interactions is normally not affected by the non-linear response; whereas it is a common problem for bulk methods (2).
Limitations of AFM
The primary assumption in using AFM to determine proteinDNA association constants is that the populations of bound and free DNA on the surface are the same as those in solution, i.e. that deposition on the surface does not alter the populations (23). If the free DNA deposits more efficiently onto the surface than the proteinDNA complexes, or if the surface causes the protein to dissociate from the DNA, the apparent binding constant determined by AFM will be less than the actual constant. One case where we have encountered this problem is for proteins that induce a 3D topology in the DNA such that the proteinDNA complex must be distorted to lie flat on the surface (M. Guthold, O.K. Wong, J. Gelles and D.A. Erie, unpublished data). Another case is that rinsing the AFM surface during sample preparations may cause the dissociation of proteinDNA complexes, especially for small proteins that have limited interaction with the surface when bound to DNA (M. Guthold, O.K. Wong, J. Gelles and D.A. Erie, unpublished data). In this latter case, the protein is often seen on the surface near the DNA (M. Guthold, O.K. Wong, J. Gelles and D.A. Erie, unpublished data), because the protein generally binds to the surface after dissociating from the DNA. In contrast, it is possible that the proteinDNA complex may deposit more efficiently than the free DNA, in which case the apparent binding constant would be overestimated. This latter case may be a problem if there is a very high occupancy of protein on the DNA; however, for long DNAs with low occupancy of protein, it is unlikely that the protein will change the efficiency of deposition of the DNA. For example, DNA fragments with streptavidinhorseradish peroxidase fusion protein bound to both ends exhibit the same rate and efficiency of deposition as the unbound fragments and no significant change in conformation (23,24). The binding constants could also be overestimated if the amount of the protein deposited on the surface is so high that there is a high probability that the protein coincidentally lands on the DNA (see below). This latter case may become a problem if the binding interaction is weak and a high concentration of protein is required to observe binding. In general, AFM should be a good method for determining binding constants as long as the occupancy of protein on the DNA is low and the protein does not fold the DNA into a 3D structure. Notably, our method does not depend on the relative populations of free protein and proteinDNA complexes, but only on that of free DNA and proteinDNA complexes. In addition, although both selective and systemic alterations in the DNA occupancy by the surface can affect absolute constants, only biased alterations between the occupancy on the specific site and that on nonspecific sites will affect specificities. Such bias could occur if the conformations of the specific and nonspecific proteinDNA complexes are significantly different.
Estimate of maximal contribution of random landing events to observed occupancies
The probability that a protein lands close enough to the DNA to be counted as a complex can be estimated from the size of the protein and the surface area covered by the DNA. To estimate the maximum possible contribution of a protein randomly landing on the DNA to the total observed complexes, we use a simple model. If we model the DNA as a cylinder and the protein as a sphere, we can define an area on the surface around the DNA in which a protein would be defined as bound. For a cylinder of width w and length L and a sphere of diameter d, the area is approximately (d + w) x (d + L). Accordingly, the probability of a single protein randomly landing within this area is
![]() | (15) |
| MATERIALS AND METHODS |
|---|
|
|
|---|
MutS protein and DNA substrates
Taq MutS was expressed and purified as described previously (11,26). Three linear DNA fragments with blunt DNA ends, named 782Homo, 783TBulge and 982GT, were generated and purified as described previously (13), where the number represents the length of DNA in base pairs followed by the text representing homoduplex DNA or the single mismatch at a dedicated position on the DNA. For 783TBulge, an extra T is 213 bp (27%) away from its closest DNA terminus. For 982GT, a GT mismatch is 412 bp (42%) away from its closest terminus. Two other linear homoduplex DNA fragments with 3'-overhang DNA ends, named 817Puc18 and 1869Puc18, were obtained by digesting Puc18 circular plasmid with the restriction enzyme, DrdI, where 817 and 1869 are the length of the resulted linear fragments. Two fragments were purified together with GFX DNA purification kit (Amersham Pharmacia Biotech), and incubated together with MutS, but they were identified and analyzed separately by their different lengths in AFM images.
AFM imaging and analysis
The MutSDNA reactions were carried out by incubating 1225 nM Taq MutS (dimer) with 0.43 nM DNA substrates (double strand) for 13 min at room temperature (23°C) or at 65°C in a binding buffer of 20 mM HEPES, pH 7.8, 50 mM NaCl and 5 mM MgCl2. The different incubation times did not affect the measured DNA occupancies, suggesting that equilibrium was established prior to deposition. The reaction mixture was deposited onto freshly cleaved mica surface (Spruce Pine Mica Company) at 23 or 65°C, incubated for <1 min on the mica, rinsed with deionized water, dried under a gentle flow of nitrogen and imaged as described previously (13). The extent of rinsing did not have significant effect on the population or conformation of the complexes on the surface. Specifically, comparison of images where the surface has been rinsed only a couple of times with those that have been rinsed dozens of times does not result in any significant differences in coverage. It is possible that some complexes were lost in the initial rinse; however, analyses of protein and DNA binding to mica indicate that they bind irreversibly over the time scale of our depositions (2225). All images were captured with a Nanoscope IIIa microscope (Digital Instruments) in oscillating mode. Pointprobe® oscillating mode silicon probes (Molecular Imaging Cooperation) with spring constants of
50 Nm1 and resonance frequencies of
170 kHz were used. All images were collected at the scan size of 1 µm x 1 µm, scan speed of 3 Hz and resolution of 512 x 512 pixels.
Generation of position distributions and determination of binding specificities
The MutSDNA complexes were selected for the position measurement independent of whether or not multiple proteins were bound on a single DNA fragment. Only MutS molecules that completely overlapped with the DNA were counted as complexes. Only complexes in which DNA contour lengths were within the standard deviation of the DNA length were used in the position histograms. The distance from the center of complex i to its closest DNA terminus (di) and the contour length of DNA i (Li) were measured in the program Nanoscope 5.12r3 (Digital Instruments). The reproducibility of this measurement was confirmed by conducting measurements in triplicate, which indicated that the error in determining the length of a given fragment is
1%. The position of the complex i is defined by Xi = di/Li. From a large number of complexes, the position histograms were plotted between the positions from 10 to 50% away from the closest end, because two DNA ends are not distinguishable in AFM images and the shortest distance to ends was used to define the position. The histograms were presented as occurrence probability
versus position, where i represents the position of the individual bins,
is the total number of binding occurrences observed within the position range (1050% DNA full length away from the closest DNA end), and Nbp,bin is the number of DNA base pairs in each position bin. Only complexes that were
10% from the end were counted because of the larger error in determining the absolute positions for the complexes close to the ends. This procedure makes the assay of specificities more efficient and has little effect on the characterization of specificities (see the discussion under Equation 8). The program Kaleidagraph (Synergy Software) was used to fit the position histograms into the position distributions. For DNA fragments containing a mismatch site (the specific site), the equation
was used for the fitting, where m1m4 were the fitting variables. This equation represents a weighted sum of one uniform distribution (P = m1) and one Gaussian distribution
. This statistical analysis is reasonable because MutS binding to mismatch-containing DNA fragments can be viewed as the sum of the binding to the homoduplex DNA and to a single mismatch. From position distributions, binding specificities were determined using Equation 8.
Determination of DNA fractional occupancies and binding constants
The occurrence probability in position histograms was not directly used for deducing the actual occupancies of individual DNA sites, because the occurrence probability at a given site will be affected by the occurrence probability at nearby sites due to the error in distance measurements. Instead, we have used the following counting of proteinDNA complexes and DNA fragments. The total number of DNA fragments (nFragment), the total number of MutSDNA complexes at internal DNA contours (nComplex,Int) and the total number of DNA termini bound by MutS (nComplex,Ter) were counted from a set of AFM images for each MutSDNA reaction. The DNA fragments on the edge of images, partial DNA fragments and overlapping DNA fragments were excluded from the counting. The fractional occupancy of the DNA fragment (OFragment), the apparent fractional occupancy of DNA termini (OTer) and the average fractional occupancy of internal binding sites (
) were determined by
and
. The number of binding sites at internal DNA contours (Nint) was assumed to be Nbp50, because 25 bp at the vicinity of each DNA terminus were allocated as the range of end binding based on the size of MutS in AFM images, due to the resolution limitation of AFM. To deduce the fractional occupancies of specific, nonspecific and end binding sites from OTer and
, DNA-binding sites of MutS were categorized into the internal mismatch site (the specific site), internal homoduplex sites (nonspecific sites) and DNA ends. For MutS interacting with homoduplex DNA, the nonspecific fractional occupancy is
, whereas for MutS interacting with the mismatched DNA, the nonspecific and specific fractional occupancies can be partitioned from
based on the specificity and the number internal binding sites (Nint):
![]() |
![]() |
. It is necessary to subtract 25 x ONSP from OTer to calculate the occupancies at the DNA ends (OE) because nonspecific binding of MutS at up to 25 sites away from the ends were likely counted as end binding, due to the size of MutS in AFM images. The relative standard errors of the occupancies were determined by at least three independent measurements.
To estimate the maximum possible contribution of proteins randomly landing on DNA to the observed occupancies on the DNA fragments, we used Equation 15, the diameter of DNA (w = 2.5 nm) and the average diameter of MutS based on the crystal structure (d = 7 nm) and the length (L) of the DNA used. For the highest concentration of protein used for the 23°C depositions (20 nM), an average of
100 proteins are seen in a 1 µm x 1 µm image (nP = 100 and ATOT = 1 µm2). For the GT-containing DNA, L is 320 nm. Accordingly, nPPP/OFragment = 0.22 (OFragment = 1.4, Table A in Supplementary Material).
As can be seen from the inspection of Equations 57, to calculate site-specific binding affinities, it is necessary to know the binding site size of the protein (Nbp,P) and the total concentrations of protein, [P], and DNA, [D], in addition to the site occupancies. Nbp,P was set to be 15 based on the inspection of the crystal structures of MutSDNA complexes (2729). [P] was calculated as the dimer concentration because Taq MutS exists primarily as a dimer in the absence and presence of DNA at the concentration used in this study (13,26,30). The relative errors of the binding constant and specificity were calculated based on the error of occupancies using Equations 1012, from which the standard errors were obtained by multiplying the relative error with the mean. The binding constants determined at different concentrations were within error of one another and the reported constants and standard deviations are the average of the constants determined from different concentrations.
DNA substrates and procedures for fluorescence measurements
DNA substrates for fluorescence measurements were purchased high-performance liquid chromatography (HPLC)-purified from MWG Biotech, Inc. and Integrated DNA Technologies, Inc. TAMRA-labeled ssDNA (5'-TACCTCATCTCGAGCGTGCCGATA-TAMRA-3') was annealed with complementary strands to create T-bulge dsDNA (5'-TATCGGCACGTCTCGAGATGAGGTA-3'), GT mismatch dsDNA (5'-TATCGGCACGTTCGAGATGAGGTA-3') and homoduplex dsDNA (5'-TATCGGCACGCTCGAGATGAGGTA-3'). The oligonucleotides were annealed in buffer containing 50 mM HEPES, pH 7.8, 100 mM NaCl and 5 mM MgCl2 in a 1:1 ratio at 55°C for 20 min then slowly cooled to room temperature. Binding reactions were performed at 23°C in the same binding buffer as used for the AFM experiments. DNA concentrations used were between 5 and 100 nM. Taq MutS was incubated with DNA for 5 min prior to measurement acquisition.
Fluorescence anisotropy was measured using a Jobin Yvon Horiba Fluorolog-3 fluorometer in T-format equipped with a Wavelength Electronics temperature control box. TAMRA-labeled dsDNA substrates were excited at 535 nm and emission was measured at 582 nm. Excitation and emission slit widths were set between 5.0 and 7.0 nm. Intensity measurements were corrected for the dark photon count, and fluorescence anisotropy was calculated using the software provided by the instrument. The fluorescence anisotropy was measured as a function of MutS concentration.
Fluorescence binding data analysis
The fluorescence anisotropy is plotted as a function of the MutS (dimer) concentration. The binding curves for the T-bulge DNA and the 10 nM GT-DNA were fit by a weighted nonlinear regression to a binding isotherm using
![]() | (16) |
versus
, where
is the fraction of DNA bound by MutS) (data not shown). This observation indicates that there are two binding events on the GT-DNA. Fitting the binding curves to two binding constants directly was not possible because the fits did not converge. Consequently, we used a both plots of
versus [MutSFree] and Scatchard plots and to estimate the high affinity binding constant. For the Scatchard plots, the constant was estimated by fitting the linear portion of the plot to a straight line. For
versus [MutSFree], the curves were fit to the sum of two binding curves [
1[MutSFree]/([MutSFree] + Kd1) + (
2[MutSFree]/([MutSFree] + Kd2))]. The reported constant is the average from both of these fits. The high affinity constants determined from these fits are consistent with the binding constants from the fits of the 10 nM data to Equation 15. For homoduplex DNA, the change in anisotropy of the DNA upon addition of MutS was too small to be able to obtain binding constants. | RESULTS AND DISCUSSION |
|---|
|
|
|---|
Binding specificities and contributions of microscopic constants to macroscopic constants
Representative AFM images of free DNA and DNA in the presence of Taq MutS are shown in Figure 2. MutS can be seen bound to the DNA in the deposition in the presence of protein. From a large number of such images, position distributions for Taq MutS bound to several different DNA fragments were obtained (Figure 3). As expected, the distributions of MutS bound to DNAs that do not contain mismatches (nonspecific DNA) are uniform, whereas those for MutS bound to DNAs that contain a single mismatch (specific DNA) are represented by the sum of a uniform and a Gaussian distribution. From these data, we calculated the specificities of Taq MutS binding to a T-bulge and a GT mismatch using Equation 8 (Table 1). In addition, we determined the standard errors of the specificity using Equation 12, and the relative accuracy is
10% (Table 1). The binding specificity of MutS can also be calculated from the binding probability on the nonspecific sites (Pmin) using Equation 9 or from directly estimating the number of specific and nonspecific complexes (16,17). The obtained specificities using these two alterative methods (data not shown) are consistent with those obtained using Equation 8, which integrates the areas under the Gaussian distribution.
|
|
|
Inspection of Table 1 reveals that Taq MutS has a larger binding specificity to a T-bulge than to a GT mismatch (1660 versus 300), which is consistent with data from gel shift assays (11). The specificity for a T-bulge is the same for AFM and gel shift assays (1660 for AFM versus 1700 for EMSA). The consistency between AFM and EMSA for a T-bulge suggests that both methods are suitable for determining the specificity of high affinity proteinDNA interactions. We have also used our method to analyze two other proteinDNA interactions for which position distributions have been published. In one study, binding of the human DNA damage recognition complex, XPC-HR23B, to an 800 bp DNA containing a single cholesterol moiety was investigated (31). Analyzing their raw position histogram yields a specificity of
2600 to a cholesterol moiety, which is consistent with the biochemical studies in a similar system (32). In another recent study, human 8-oxoguanine DNA glycosylase (hOGG1) binding to a 1024 bp linear DNA containing a single oxoG was investigated (17). Analysis of their raw position histogram yields a binding specificity of 390 for hOGG1 to the single oxoG, which is the same as their estimation of 400 and is consistent with bulk measurements (17). Taken together, these results indicate that AFM is a good method for directly determining the specificities of proteinDNA interactions. In addition, these results support our suggestion that proteinDNA interactions with different dynamics are detected with similar efficiencies using AFM, because the nonspecific complexes are more dynamic than the specific ones. Interestingly, the specificity for Taq MutS binding to a GT mismatch determined by AFM is significantly higher than that determined by EMSA (300 for AFM versus 12 for EMSA). AFM provides a direct measure of specificity, in that the relative probability of the protein bound to different sites on a single DNA fragment are compared, and the relative probabilities are a direct measure of the relative affinities for the different sites. In bulk studies, specificity is determined by measuring the binding affinity to two different DNA fragments: one containing a specific site and one without a specific site. Consequently, to understand the difference between the specificities determined by AFM and EMSA, it is necessary to inspect the contributions of different types of binding to macroscopic binding constants (Equations 13 and 14).
From fractional occupancies of MutS bound to different DNA sites (Supplementary Table A), we calculated the binding affinities of Taq MutS to specific and nonspecific sites and to DNA ends using Equations 57, as well as the apparent macroscopic binding affinities using Equations 13 and 14 (Table 1). Inspection of Table 1 reveals that Taq MutS has a weak binding affinity for nonspecific DNA sites (2035 µM; the agreement between Puc18 and 782Homo suggests that the engineered DNA fragments are similar to plasmid DNA fragments) and that nonspecific binding makes only a small contribution to the calculated macroscopic binding constants (KDNA,AFM) for a DNA fragment containing a T-bulge or a GT mismatch (Table 1). In contrast, Taq MutS has a high affinity for DNA ends (
50 nM), similar to that for a GT mismatch (77 nM). Consistent with this result, analysis of the position distribution of E.coli MutS on a DNA fragment containing a GT mismatch indicates that E.coli MutS binds to DNA ends with an affinity that is only approximately five times less than that of a GT mismatch [Supplementary Material in ref. (13); H. Wang, P. Hsieh and D.A. Erie, unpublished data]. End binding makes an increasingly important contribution to KDNA,AFM as specificity decreases, with KDNA,AFM for nonspecific DNA being dominated by the contribution from end binding (Table 1). Consequently, macroscopic binding constants determined using bulk methods may vary significantly depending on whether or not end binding is detected in the assay. As discussed in the following sections, the apparent differences in specificities of Taq and E.coli MutS determined from bulk measurements likely result from differences in the extent to which end binding is being detected in bulk assays.
Comparison of AFM and bulk studies
To compare the AFM results with the existing data from bulk studies, we have calculated the macroscopic binding constants (KDNA,AFM) of Taq MutS to 60 bp DNA fragments (the length of the DNA fragment in the bulk EMSA study) using Equations 13 and 14 (Table 1). Comparison of the calculated macroscopic AFM binding constants, KDNA,AFM, with the apparent binding constants determined from gel shift assays, KDNA,EMSA, reveals that the constants for Taq MutS binding to DNA containing a T-bulge, which is the highest affinity mismatch, are very similar (12 nM for AFM versus 2 nM for EMSA). Interestingly, AFM, however, yields higher binding affinities than the EMSA studies for DNA containing a GT mismatch (18 nM for AFM versus 310 nM for EMSA) and for homoduplex DNA (23 nM for AFM versus 3800 nM for EMSA). The consistency between the results from AFM and EMSA for MutS binding to a T-bulge indicates that both AFM and gel shift assays are good methods for determining this tight binding constant. The differences for homoduplex DNA and GT-containing DNA, however, indicate that either AFM is overestimating the binding constants or that the EMSA experiments are underestimating them. As discussed below, it is more likely that the gel shift assays underestimate the binding constants for the GT-mismatch and homoduplex DNA fragments.
For AFM to yield an apparent binding constant that is greater than the actual binding constant, the free DNA would have to be less efficiently deposited on the surface than the DNA with protein bound. The agreement between the AFM and EMSA data for MutS binding to the DNA containing a T-bulge suggests that there is no significant difference between the deposition efficiency of free and bound DNA for this fragment. Consequently, it is unlikely that there would be a difference for free and bound DNA containing a GT mismatch or no mismatch, because both the occupancies of MutS on the DNA (Supplementary Table A) and the overall conformations of the proteinDNA complexes are similar (13). In addition, as mentioned above the rates and efficiency of deposition of DNA fragments with protein bound on both ends are the same as those of free DNA (23,24). Another way in which AFM could overestimate the binding affinity is if the coverage of protein on the surface is too high, such that the protein coincidently lands on DNA. Because the protein coverage on the surface is low in these studies (Figure 2), this potential problem is minimized. Specifically, using a statistical analysis (see Theory), we have estimated the maximal contribution of proteins randomly landing on the DNA to the observed binding constants. This analysis indicates that for the highest protein and DNA concentrations used, random landing could result in no more than
20% increase in the observed binding constant (see Materials and Methods), which is insufficient to explain the observed discrepancies between AFM and EMSA. Furthermore, the binding constants of MutS for homoduplex sites and for DNA ends determined from different DNA and protein concentrations are consistent with one other (within 0.5 kcal/mol) and show no trend with protein or DNA concentration (Table 1 and Supplementary Table A), indicating that the contribution from proteins randomly landing on the DNA is negligible. As discussed in Theory, it is more likely that this AFM method would underestimate a binding constant than overestimate it.
The above discussion leads to the suggestion that the gel shift assays may be underestimating the binding affinities of Taq MutS to GT-containing and homoduplex DNA fragments. One significant limitation to using gel shift assays for determining binding constants is its insensitivity to weak or dynamic proteinDNA interactions (21). Notably, the AFM and EMSA results agree for the T-bulge, which has the highest binding affinity. As discussed above, MutS binding to DNA ends makes an increasingly important contribution to the calculated macroscopic binding constants, KDNA,AFM, as specificity decreases (Table 1). Consequently, if binding of Taq MutS to DNA ends is not detected in the gel shift assays, the apparent binding constants for nonspecific and lower specificity DNA fragments determined by EMSA would be significantly less than those calculated from AFM. To assess this possibility further, we calculated apparent binding constants from AFM ignoring the contribution from end binding (Table 1, italic values). Although the affinities are still higher than that from EMSA (68 nM for AFM versus 310 nM for EMSA for GT-containing fragment; 495 nM for AFM versus 3800 nM for EMSA for homoduplex DNA fragment), they are similar given the different techniques and different DNA fragments that were used in the two studies.
To further resolve this discrepancy, we have measured the binding affinity, using fluorescence anisotropy (33,34), of MutS to short (24 bp) DNA fragments that have a fluorescent label on one end. Figure 4 shows binding curves for T-bulge and GT-containing DNA. The curves for T-bulge DNA are fit well to a single binding isotherm (Figure 4). The binding constant for T-bulge DNA, determined from the average of seven measurements, is 5 ± 4.9 nM, which is in good agreement with both the AFM and EMSA results. This result indicates that both AFM and EMSA yield accurate binding constants for MutS binding to DNA containing a T-bulge. The binding curves for the GT-DNA at 10 nM are fit well to a single binding isotherm (Figure 4); however, at 30 and 50 nM GT-DNA, they are biphasic, which is clearly revealed by non-linearity in Scatchard plots (data not shown). This result indicates that there are multiple binding events on the GT-DNA, which is consistent with the AFM data, which yields similar binding constants for DNA ends and a GT mismatch. The binding constant for GT-DNA, determined from the average of four independent measurements (see Materials and Methods), is 40 ± 25 nM. This binding affinity is in good agreement with the binding constant determined by AFM and is significantly tighter than that determined by EMSA. For homoduplex DNA, the change in anisotropy upon the addition of MutS was too small to determine a binding constant (data not shown).
|
Taken together, these analyses strongly suggest that end binding of Taq MutS was not detected by EMSA and that the binding to the weaker nonspecific sites and to a GT mismatch was only partially detected. On the other hand, AFM is expected to be significantly less sensitive to the dynamics of the proteinDNA interactions because the deposition of the complexes is rapid and irreversible over the time scale of the deposition (23,24). This suggestion is supported both by the agreement between specificities calculated from biochemical and AFM data (see previous section) and by the observation that the binding constant for Taq MutS to nonspecific DNA is similar at 23 and 65°C (Table 1). In summary, these results strongly suggest that the differences between AFM and EMSA are a result of the dynamic limitations of EMSA and that AFM provides an accurate measure of the binding constants.
Binding specificities of E.coli and Taq MutS
E.coli MutS presents an interesting counter example to Taq MutS, in that end binding by E.coli MutS appears to be detected by EMSA. Consistent with our AFM studies (H. Wang, P. Hsieh and D.A. Erie, unpublished data), recent biochemical studies indicate that E.coli MutS has a strong binding affinity to DNA ends and suggest that this end binding is stable to electrophoresis (35). This suggestion is further supported by inspection of the gels from the EMSA study, which reveals a super-shifted band on increasing E.coli MutS concentration, indicating that two MutS proteins are bound to the fragments at higher concentrations (11). In contrast, no super-shifted bands are seen with Taq MutS at a concentration of 100 nM (13), supporting our assertion that end binding of Taq MutS is not stable to electrophoresis. Detection by EMSA of end binding by E.coli MutS but not by Taq MutS likely results from the interactions between E.coli MutS and DNA ends being less dynamic than those between Taq MutS and DNA ends. Consistent with this suggestion, kinetic studies indicate that the dissociation rate of Taq MutS from a T-bulge (substrate with the highest affinity) is four times faster than that of E.coli MutS from a T-bulge (36).
The differences in the detection of end binding provide an explanation for the apparent differences in the specificities of E.coli and Taq MutS determined from EMSA studies. The bulk specificity for a T-bulge (KT-bulge,EMSA/KHomo,EMSA) is
1700 for Taq MutS, whereas it is only
30 for E.coli MutS (Table 1). The differences in these specificities result primarily from differences in the apparent binding affinities to homoduplex DNA (30 nM for E.coli MutS versus 3000 nM for Taq MutS). The higher apparent binding affinity of E.coli MutS for homoduplex DNA results from stable end binding and is predicted by our analysis of the microscopic binding constants and the observations of high binding affinity of E.coli MutS for DNA ends. Consequently, for E.coli MutS, the ratio KT-bulge,EMSA/KHomo,EMSA is not a good measure of specificity. It is likely that E.coli and Taq MutS have similar specificities for a T-bulge, given that their EMSA binding affinities are similar (Table 1). Published AFM data of the position distribution of E.coli MutS bound to DNA fragments containing a GT mismatch support this idea. Analysis of these data [see Supplementary Material in ref. (13)] reveals that E.coli MutS has a specificity of
1000 for a GT mismatch, which is similar to that for Taq MutS (
300) but significantly higher than the specificity (KT-bulge,EMSA/KHomo,EMSA) determined by EMSA (
7).
Our results indicate that MutS has significantly higher specificities for mismatches than those previously determined from bulk measurements. Taq MutS maintains a high specificity (
200) for a T-bulge even at its physiological temperature of 65°C (Table 1). These higher specificities resolve, in part, the previous conundrum of apparent low MutS binding specificities, but highly efficient repair (13). AFM is powerful method for determining specificities, because it provides a direct measure of the relative binding affinities to different sites on the same DNA fragment and knowledge of the absolute binding affinities is not required. The only way that specificities determined from AFM studies can be skewed is if there is a difference in deposition efficiency of specific and nonspecific complexes, which is unlikely unless there is a significant difference in the overall conformation of the specific and nonspecific complexes.
| CONCLUSION |
|---|
|
|
|---|
Intrinsic limitations of bulk techniques on the quantification of proteinDNA binding constants and specificities have been long recognized, especially for proteins with low binding specificity and/or a second specific site, such as DNA ends. Only apparent binding constants to the entire DNA fragment can be determined from these methods. Furthermore, many of these bulk techniques are sensitive to the dynamics of the interactions. We have established a single-molecule method for determining proteinDNA binding constants and specificities in a site-specific fashion. This method provides direct measure of the binding affinity to a particular site because the occurrence of a protein binding to each site on the DNA is directly observed. Similarly, the binding specificities from AFM are site-to-site (microscopic) comparisons of binding constants instead of fragment-to-fragment comparisons. In addition, this method is straightforward because it relies only on counting and 1D-distance measurements of hundreds of complexes, which are sufficient for the assay. Finally, from a single set of AFM experiments, it is possible to determine the binding affinity, specificity and stoichiometry, as well as the conformational properties of the proteinDNA complexes (23,37). Application of this method to MutSDNA interactions has revealed that the apparent differences between Taq and E.coli MutS DNA binding properties are due to the strong end-binding affinity of MutS, as well as to a problem in detecting end binding and nonspecific DNA binding by Taq MutS in the EMSA measurements. The significantly higher DNA-binding specificities of MutS obtained in this study are obviously important for achieving efficient DNA repair in the cell. The biological significance of strong binding affinity for DNA ends is unclear; however, it may be functionally important in other biological processes in which MutS is involved, such as DNA double-strand break repair and recombination, because DNA ends are critical intermediates in these pathways.
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at NAR Online.
| ACKNOWLEDGEMENTS |
|---|
The authors are grateful for helpful discussions with Xiaolei Zhou. The authors thank Dr Ashutosh Tripathy, Director, UNC Macromolecular Interactions Facility, for use of fluorescence instrumentation. This work was supported by the National Institutes of Health grant GM R01-54316 and the American Cancer Society (DAE).
Conflict of interest statement. None declared.
| REFERENCES |
|---|
|
|
|---|
- Carey, J. (1991) Gel retardation Methods Enzymol., 208, 103117[Medline] .
- Lohman, T.M. and Bujalowski, W. (1991) Thermodynamic methods for model-independent determination of equilibrium binding isotherms for proteinDNA interactions: spectroscopic approaches to monitor binding Methods Enzymol., 208, 258290[Medline] .
- Myszka, D.G. (2000) Kinetic, equilibrium, and thermodynamic analysis of macromolecular interactions with BIACORE Methods Enzymol., 323, 325340[CrossRef][Web of Science][Medline] .
- Oda, M., Furukawa, K., Ogata, K., Sarai, A., Nakamura, H. (1998) Thermodynamics of specific and non-specific DNA binding by the c-Myb DNA-binding domain J. Mol. Biol., 276, 571590[CrossRef][Web of Science][Medline] .
- Winzor, D.J. and Sawyer, W.H. Quantitative Characterization of Ligand Binding, (1995) New York, NY Wiley-Liss, Inc .
- Record, M.T., Jr, Ha, J.H., Fisher, M.A. (1991) Analysis of equilibrium and kinetic measurements to determine thermodynamic origins of stability and specificity and mechanism of formation of site-specific complexes between proteins and helical DNA Methods Enzymol., 208, 291343[Medline] .
- Modrich, P. (1989) Methyl-directed DNA mismatch correction J. Biol. Chem., 264, 65976600
[Abstract/Free Full Text] . - Modrich, P. and Lahue, R. (1996) Mismatch repair in replication fidelity, genetic recombination, and cancer biology Annu. Rev. Biochem., 65, 101133[CrossRef][Web of Science][Medline] .
- Hsieh, P. (2001) Molecular mechanisms of DNA mismatch repair Mutat. Res., 486, 7187[Web of Science][Medline] .
- Schofield, M.J. and Hsieh, P. (2003) DNA mismatch repair: molecular mechanisms and biological function Annu. Rev. Microbiol., 57, 579608[CrossRef][Web of Science][Medline] .
- Schofield, M.J., Brownewell, F.E., Nayak, S., Du, C., Kool, E.T., Hsieh, P. (2001) The Phe-X-Glu DNA binding motif of MutS. The role of hydrogen bonding in mismatch recognition J. Biol. Chem., 276, 4550545508
[Abstract/Free Full Text] . - Gradia, S., Acharya, S., Fishel, R. (2000) The role of mismatched nucleotides in activating the hMSH2hMSH6 molecular switch J. Biol. Chem., 275, 39223930
[Abstract/Free Full Text] . - Wang, H., Yang, Y., Schofield, M.J., Du, C., Fridman, Y., Lee, S.D., Larson, E.D., Drummond, J.T., Alani, E., Hsieh, P., et al. (2003) DNA bending and unbending by MutS govern mismatch recognition and specificity Proc. Natl Acad. Sci. USA, 100, 1482214827
[Abstract/Free Full Text] . - Wang, H. and Hays, J.B. (2002) Mismatch repair in human nuclear extracts. Quantitative analyses of excision of nicked circular mismatched DNA substrates, constructed by a new technique employing synthetic oligonucleotides J. Biol. Chem., 277, 2613626142
[Abstract/Free Full Text] . - McGhee, J.D. and von Hippel, P.H. (1974) Theoretical aspects of DNAprotein interactions: co-operative and non-co-operative binding of large ligands to a one-dimensional homogeneous lattice J. Mol. Biol., 86, 469489[CrossRef][Web of Science][Medline] .
- de Jager, M., Wyman, C., van Gent, D.C., Kanaar, R. (2002) DNA end-binding specificity of human Rad50/Mre11 is influenced by ATP Nucleic Acids Res., 30, 44254431
[Abstract/Free Full Text] . - Chen, L., Haushalter, K.A., Lieber, C.M., Verdine, G.L. (2002) Direct visualization of a DNA glycosylase searching for damage Chem. Biol., 9, 345350[CrossRef][Web of Science][Medline] .
- Solis, F.J., Bash, R., Yodh, J., Lindsay, S.M., Lohr, D. (2004) A statistical thermodynamic model applied to experimental AFM population and location data is able to quantify DNAhistone binding strength and internucleosomal interaction differences between acetylated and unacetylated nucleosomal arrays Biophys. J., 87, 33723387[CrossRef][Web of Science][Medline] .
- Taylor, J.R. An Introduction to Error Analysis: The Study of Uncertainties in Physical Measurements, (1997) 2nd edn Sausalito, CA University Science Books .
- Koch, S.J., Shundrovsky, A., Jantzen, B.C., Wang, M.D. (2002) Probing proteinDNA interactions by unzipping a single DNA double helix Biophys. J., 83, 10981105[Web of Science][Medline] .
- Lim, W.A., Sauer, R.T., Lander, A.D. (1991) Analysis of DNAprotein interactions by affinity coelectrophoresis Methods Enzymol., 208, 196210[Medline] .
- Lee, C.S. and Belfort, G. (1989) Changing activity of ribonuclease A during adsorption: a molecular explanation Proc. Natl Acad. Sci. USA, 86, 83928396
[Abstract/Free Full Text] . - Ratcliff, G.C. and Erie, D.A. (2001) A novel single-molecule study to determine proteinprotein association constants J. Am. Chem. Soc., 123, 56325635[CrossRef][Web of Science][Medline] .
- Rivetti, C., Guthold, M., Bustamante, C. (1996) Scanning force microscopy of DNA deposited onto mica: equilibration versus kinetic trapping studied by statistical polymer chain analysis J. Mol. Biol., 264, 919932[CrossRef][Web of Science][Medline] .
- Gettens, R.T., Bai, Z., Gilbert, J.L. (2005) Quantification of the kinetics and thermodynamics of protein adsorption using atomic force microscopy J. Biomed. Mater. Res. A, 72A, 246257[Medline] .
- Biswas, I., Ban, C., Fleming, K.G., Qin, J., Lary, J.W., Yphantis, D.A., Yang, W., Hsieh, P. (1999) Oligomerization of a MutS mismatch repair protein from Thermus aquaticus J. Biol. Chem., 274, 2367323678
[Abstract/Free Full Text] . - Obmolova, G., Ban, C., Hsieh, P., Yang, W. (2000) Crystal structures of mismatch repair protein MutS and its complex with a substrate DNA Nature, 407, 703710[CrossRef][Medline] .
- Lamers, M.H., Perrakis, A., Enzlin, J.H., Winterwerp, H.H., de Wind, N., Sixma, T.K. (2000) The crystal structure of DNA mismatch repair protein MutS binding to a G x T mismatch Nature, 407, 711717[CrossRef][Medline] .
- Natrajan, G., Lamers, M.H., Enzlin, J.H., Winterwerp, H.H., Perrakis, A., Sixma, T.K. (2003) Structures of Escherichia coli DNA mismatch repair enzyme MutS in complex with different mismatches: a common recognition mode for diverse substrates Nucleic Acids Res., 31, 48144821
[Abstract/Free Full Text] . - Bjornson, K.P., Allen, D.J., Modrich, P. (2000) Modulation of MutS ATP hydrolysis by DNA cofactors Biochemistry, 39, 31763183[CrossRef][Medline] .
- Janicijevic, A., Sugasawa, K., Shimizu, Y., Hanaoka, F., Wijgers, N., Djurica, M., Hoeijmakers, J.H., Wyman, C. (2003) DNA bending by the human damage recognition complex XPC-HR23B DNA Repair (Amst.), 2, 325336[CrossRef][Medline] .
- Hey, T., Lipps, G., Sugasawa, K., Iwai, S., Hanaoka, F., Krauss, G. (2002) The XPCHR23B complex displays high affinity and specificity for damaged DNA in a true-equilibrium fluorescence assay Biochemistry, 41, 65836587[CrossRef][Medline] .
- Jabonski, A. Bull Acad Pol Sci Ser A, (1960) 8, 259264 .
- Lakowicz, J.R. Topics in Fluorescence Spectroscopy, (1999) New York Kluwer Academic/Plenum Publishers pp. 291318 Chapter 10 .
- Acharya, S., Foster, P.L., Brooks, P., Fishel, R. (2003) The coordinated functions of the E. coli MutS and MutL proteins in mismatch repair Mol. Cell, 12, 233246[CrossRef][Web of Science][Medline] .
- Schofield, M.J., Nayak, S., Scott, T.H., Du, C., Hsieh, P. (2001) Interaction of Escherichia coli MutS and MutL at a DNA mismatch J. Biol. Chem., 276, 2829128299
[Abstract/Free Full Text] . - Yang, Y., Wang, H., Erie, D.A. (2003) Quantitative characterization of biomolecular assemblies and interactions using atomic force microscopy Methods, 29, 175187[CrossRef][Web of Science][Medline]
.
This article has been cited by other articles:
![]() |
M. Labazi, L. Jaafar, and H. Flores-Rozas Modulation of the DNA-binding activity of Saccharomyces cerevisiae MSH2-MSH6 complex by the high-mobility group protein NHP6A, in vitro Nucleic Acids Res., October 20, 2009; (2009) gkp649v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Wagner, G. Moolenaar, J. van Noort, and N. Goosen Single-molecule analysis reveals two separate DNA-binding domains in the Escherichia coli UvrA dimer Nucleic Acids Res., April 1, 2009; 37(6): 1962 - 1972. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Tessmer, Y. Yang, J. Zhai, C. Du, P. Hsieh, M. M. Hingorani, and D. A. Erie Mechanism of MutS Searching for DNA Mismatches and Signaling Repair J. Biol. Chem., December 26, 2008; 283(52): 36646 - 36654. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
























