| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
New York University, Department of Chemistry, New York, New York 10003
Correspondence: Address reprint requests to Alexander Vologodskii, E-mail: alex.vologodskii{at}nyu.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The active form of the Escherichia coli enzyme consists of two A and two B subunits (reviewed by Reece and Maxwell, 1991
). For many years it was accepted that the tetramer binds a segment of double-stranded DNA
140 bp in length that is supposed to wrap around the enzyme forming approximately one turn of the right-handed helix (Fisher et al., 1981
; Kirkegaard and Wang, 1981
; Liu and Wang, 1978a
). The bound DNA segment, named the gate (G) segment, is cleaved approximately in the middle, and the 5' ends of the broken strands are covalently attached to the protein (Fig. 1 A). Another DNA segment, the so-called transporting (T) segment, then passes through the double-stranded break and the break is resealed (Reece and Maxwell, 1991
; Wang, 1998
). This model of the G-segment conformation is supported by the results from nuclease protection experiments (Liu and Wang, 1978a
; Morrison and Cozzarelli, 1981
) as well as by DNase I (Fisher et al., 1981
; Kirkegaard and Wang, 1981
) and hydroxyl radical footprinting (Orphanides and Maxwell, 1994
). It has been shown that binding DNA gyrase to nicked DNA in the absence of ATP increases the linking number of the DNA strands, measured after ligation of the nicks (Kampranis et al., 1999
; Liu and Wang, 1978b
; Peng and Marians, 1995
). Quantitative analysis of these experiments showed that each complex between the enzyme and DNA increases the linking number in the ligated molecules by 0.50.8. These data are considered to provide additional support for the model. If both the G- and T-segments belong to the same closed circular DNA, each strand-passing reaction results in the change of the linking number (Lk) of DNA complementary strands by two (Brown and Cozzarelli, 1979
). If the T-segment belongs to second circular DNA, the reaction can link or unlink the molecules.
|
140 bp are wrapped around the enzyme, there are only
40 bp available to make the external loops. Recent data obtained by atomic force microscopy are also in certain disagreement with the model (Heddle et al., 2004
Although DNA gyrase can catenate, decatenate, and unknot circular DNA molecules (Kreuzer and Cozzarelli, 1980
; Mizuuchi et al., 1980
), the efficiency of these reactions is very low (Marians, 1987
; Ullsperger and Cozzarelli, 1996
). The major catalytic activity of the enzyme is selective introduction of negative supercoils into circular DNA. The right-handed helix turn, formed by G-segment, has been assumed to provide this selectivity (Fig. 1 A). However, DNA molecules of a few thousand base pairs in length adopt many different conformations in solution, and it is hard to believe that a small chiral loop can provide the required selectivity. Indeed, other types of conformations (with a potential T-segment juxtaposed with the enzyme gate) would result in different topological changes. In particular, the outcome from two conformations shown in Fig. 1, B and C, would be the introduction of two positive supercoils and formation of a positive trefoil. None of these products has been observed as a reaction product under conditions optimal for supercoiling. Thus, we suggested that the model might not account for the experimentally observed selectivity of the enzyme action. In this article we analyze this question quantitatively by computer simulation.
There is convincing evidence that computer simulations provide accurate quantitative description of the large-scale conformational properties of DNA. Simulations reproduce experimental data on hydrodynamic properties of DNA molecules (Hagerman, 1981
; Hagerman and Zimm, 1981
; Rybenkov et al., 1997d
), DNA cyclization (Hagerman, 1990
; Levene and Crothers, 1986
; Taylor and Hagerman, 1990
; Vologodskaia and Vologodskii, 2002
), equilibrium distributions of topological states (Klenin et al., 1988
, 1989
; Rybenkov et al., 1993
, 1997c
; Shaw and Wang, 1993
; Vologodskii and Cozzarelli, 1993
), elasticity of single molecules (Vologodskii, 1994
; Vologodskii and Marko, 1997
), and light and neutron scattering data on supercoiled DNA (Gebe et al., 1996
; Hammermann et al., 1998
, 1997
; Klenin et al., 2000
). The simulations are based on the statistical-mechanical treatment of a well-established model of the double helix. There are only three parameters in the model and all of them have been reliably determined for various solution conditions. Thus, the simulations are able to provide reliable quantitative information on many DNA properties that are hard to measure experimentally. The simulations are very useful for the analysis of enzymatic reactions that involve formation of DNA loops (Grainge et al., 2002
; Klenin et al., 2002
; Vologodskii et al., 2001
). Here we apply the simulation to estimate the distribution of the reaction products corresponding to the model of DNA gyrase shown in Fig. 1. The results leave no doubt that the model is inconsistent with selective formation of negative supercoils by DNA gyrase.
The simulation results allow us to formulate a general requirement for a correct model of the enzyme. Models that satisfy this requirement were suggested by Kampranis et al. (1999)
and recently by Corbett et al. (2004)
. We discuss and elaborate important details of such a model.
| METHODS OF COMPUTATION |
|---|
|
|
|---|
rigid cylinders of equal length with harmonic potentials with respect to the angles between adjacent segments and the value of torsional deformation,
Tw (Klenin et al., 1991
Lk value. To preserve the chain topology during a simulation run we used a topological invariant, the Alexander polynomial,
A topological invariant has the same value over all conformations with a particular topology. We calculated
for
and
after each move of the Metropolis procedure (Vologodskii et al., 1974
or
were different for the current conformation and trial conformations, the trial conformation was rejected. A rigid subchain of five straight segments that form nearly one turn of the right-handed helix corresponded to the G-segment bounded with DNA gyrase (Fig. 2 A). The geometry of the subchain did not change during the simulation.
|
Parameters
All computations were performed for a DNA persistence length of 50 nm (Hagerman, 1988
), a torsional rigidity, C, of 3.0 x 1019 erg·cm (Horowitz and Wang, 1984
; Klenin et al., 1989
), and a DNA effective diameter, d, of 5 nm. The later value of d corresponds to a 0.2 M solution of NaCl (Brian et al., 1981
; Rybenkov et al., 1993
, 1997b
; Stigter, 1977
). The length of the straight segments was equal to 10 nm.
Analysis of the simulated conformations
A segment of the chain was considered to be juxtaposed with the G-segment if all of the following conditions were satisfied:
(Fig. 2 A);
(0°
90°) between the tested segment and the axis of the helix, which specifies the orientation of the G-segment, did not exceed
(Fig. 2 A);
for
(Vologodskii et al., 1975
If the link did not appear, we considered that there is no juxtaposition. Such cases were very rare, however.
Although the probabilities of juxtaposition depend on the values chosen for
and
for sufficiently small
and
the product distributions do not depend on their choice. In most of the calculations we used
nm and
If a particular conformation of the chain passed all tests for the juxtaposition, we determined the potential reaction product by calculating the change of chain writhe,
Wr, and the values of
or
after the extension of the G-segment shown in Fig. 2 C. The Alexander polynomial, however, does not distinguish between a knot and its mirror image. To solve this problem, we used
Wr (Kleninet al., 2002
). It can be proven that, if the substrate is an unknotted conformation and the product is a trefoil, then the "sign" of the product coincides with the sign of
Wr. Similarly, if the substrate is a trefoil and the product is an unknot, then the "sign" of the substrate is opposite to the sign of
Wr (see, e.g., Murasugi, 1996
; in particular Example 6.4.2 and Theorem 6.4.7).
| RESULTS |
|---|
|
|
|---|
Lk, and by the topology of DNA axis, which can be unknotted or form various knots. Averaging the results of such analysis over a large volume of sampled conformations gives us the desired distribution of the reaction outcomes.
|
Lk = 0) unknotted chains. Fig. 4 A shows the results obtained for DNA molecules of different lengths. One can see from the figure that introduction of two () supercoils would be the dominant reaction outcome only for DNA molecules shorter than 2000 bp. For molecules of 3000 bp or longer, formation of positive trefoils should be comparable with the introduction of negative supercoiling. Clearly, for the relaxed molecules of these lengths the chirality of the G-segment does not provide the selectivity of reaction outcomes comparable with what is observed experimentally.
|
Fig. 5 shows the product distributions obtained from negatively supercoiled molecules. The results were obtained for a circular DNA 3500 kb in length. Although for this DNA length negative supercoiling is the major reaction outcome for all tested initial values of supercoiling, the total fraction of other products remains close to 30%. The distribution of the other products changes, however, with the increase of () supercoiling. Surprisingly, the fraction of positive trefoils, which was the second largest fraction for the relaxed chains (Fig. 4 A), declines rapidly as () supercoiling increases. Instead, an increase of Lk by (+2) becomes the second major outcome of the reaction.
|
| DISCUSSION |
|---|
|
|
|---|
What reliable experimental data have to be explained by a realistic model of the enzyme action? Wrapping the G-segment around the enzyme seems to be an established feature of the complex. A model also has to account for the fact that the double-stranded break is introduced nearly in the middle of DNA segment bound to the enzyme (reviewed in Reece and Maxwell, 1991
). Additional strong restriction for a model imposes the finding that the enzyme is able to introduce supercoiling into very small DNA circles, 174 bp in length (Bates and Maxwell, 1989
). We believe that the simplest and perhaps the only way to fit all of the above requirements for DNA minicircles corresponds to the diagram shown in Fig. 6. Here the T-segment enters the gate from inside rather than from outside the loop formed by the G-segment. This allows the G- and T-segments to be combined into a single stretch that interacts with the enzyme along its entire length. If the enzyme has twofold symmetry, as is usually assumed, the pathway of the short DNA circle keeps (or nearly keeps) this symmetry during the strand-passing reaction in the model. To generalize the model for longer DNA molecules we extend one of the two short internal loops formed by the G- and T-segments outside the enzyme surface (Fig. 7). Conformations of the DNA stretch interacting with the protein during the strand-passing reaction are completely specified by this interaction, and this provides a unique reaction outcome. It is important to note that even a small free loop between the G- and T-segments suggested by Kampranis et al. (1999)
as an alternative to the bound loop, cannot exclude other outcomes of the reaction. Indeed, there is no way to restrict the size and conformations of the free loop if its pathway is not controlled by DNA-protein interaction.
|
|
A model that assumes that both the G- and T-segments belong to a continuous DNA stretch interacting with the enzyme cannot explain, however, catalysis of the catenation/decatenation reaction by DNA gyrase. Indeed, in this case the G- and T-segments belong to different DNA molecules and cannot form a continuous stretch of the double helix. Thus, we have to conclude that there are two modes of the enzyme action. In the major mode DNA gyrase interacts with a continuous stretch of DNA and introduces negative supercoiling into closed circular DNA molecules. In the minor mode the second part of the stretch, which corresponds to the T-segment, is substituted by another DNA segment. The later segment can belong to a remote part of the same DNA or another molecule. In this second mode the enzyme is able to knot/unknot DNA molecules and catenate/decatenate DNA circles. The catalytic activity of the major mode is at least 1000 times greater than that of the minor mode. It cannot be excluded, however, that the activity of the minor mode is related with a small fraction of damaged enzymes.
Formation of knots in circular DNA has never been observed in the reaction catalyzed by DNA gyrase. According to the computer simulations performed for the model shown in Fig. 1, trefoils have to make up a large fraction of reaction products if relaxed circular molecules serve as a substrate (see Fig. 4 A). On the other hand, the fraction of trefoils that corresponds to the thermodynamic equilibrium does not exceed 3% for 7-kb DNA molecules under close to physiological ionic conditions (Rybenkov et al., 1993
). Thus, if DNA gyrase works according to this model, it would greatly increase the fraction of knots over the equilibrium level on nicked circular DNA (see Fig. 4 A). This effect would be the opposite of the property of other type II topoisomerases that strongly decrease the steady-state fraction of knots relative to the equilibrium level (Rybenkov et al., 1997a
). This is not surprising because the model shown in Fig. 1 is exactly opposite to one that was suggested to explain topology simplification by the topo II enzymes (Vologodskii et al., 2001
). The computational analysis showed that a topoisomerase reduces the fraction of knots and links below the equilibrium level if the T-segment passes through the G-segment from inside to outside the loop formed by the G-segment upon binding with the enzyme (Vologodskii et al., 2001
). Conversely, the fraction of knots and links will be increased above the equilibrium level if the T-segment passes through the G-segment from outside to inside the loop, as it is shown in Fig. 1. The data presented in Fig. 4 A reflect this feature of the gyrase model. In the model shown in Fig. 7 the T-segment passes through the break from inside to outside the loop. Thus, the minor mode of the enzyme action has to reduce the steady-state fraction of knots below equilibrium level. The fact that formation of knots by DNA gyrase has been never observed experimentally gives additional support to the suggested model. Unfortunately, the very low level of activity of the enzyme in the minor mode makes it difficult to study the issue quantitatively.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
The work was supported by National Institutes of Health grant GM 54215.
Submitted on March 17, 2004; accepted for publication August 26, 2004.
| REFERENCES |
|---|
|
|
|---|
Brian, A. A., H. L. Frisch, and L. S. Lerman. 1981. Thermodynamics and equilibrium sedimentation analysis of the close approach of DNA molecules and a molecular ordering transition. Biopolymers. 20:13051328.[CrossRef][Medline]
Brown, P. O., and N. R. Cozzarelli. 1979. A sign inversion mechanism for enzymatic supercoiling of DNA. Science. 206:10811083.
Corbett, K. D., R. K. Shultzaberger, and J. M. Berger. 2004. The C-terminal domain of DNA gyrase A adopts a DNA-bending beta-pinwheel fold. Proc. Natl. Acad. Sci. USA. 101:72937298.
Fisher, L. M., K. Mizuuchi, M. H. O'Dea, H. Ohmori, and M. Gellert. 1981. Site-specific interaction of DNA gyrase with DNA. Proc. Natl. Acad. Sci. USA. 78:41654169.
Gebe, J. A., J. J. Delrow, P. J. Heath, B. S. Fujimoto, D. W. Stewart, and J. M. Schurr. 1996. Effects of Na+ and Mg2+ on the structures of supercoiled DNAs: comparison of simulations with experiments. J. Mol. Biol. 262:105128.[CrossRef][Medline]
Gellert, M., K. Mizuuchi, M. H. O'Dea, and H. A. Nash. 1976. DNA gyrase: an enzyme that introduces superhelical turns into DNA. Proc. Natl. Acad. Sci. USA. 73:38723876.
Grainge, I., S. Pathania, A. Vologodskii, R. Harshey, and M. Jayaram. 2002. Symmetric DNA sites are functionally asymmetric within Flp and Cre site-specific DNA recombination synapses. J. Mol. Biol. 320:515527.[CrossRef][Medline]
Hagerman, P. J. 1981. Investigation of the flexibility of DNA using transient electric birefringence. Biopolymers. 20:15031535.[CrossRef][Medline]
Hagerman, P. J. 1988. Flexibility of DNA. Annu. Rev. Biophys. Biophys. Chem. 17:265286.[CrossRef][Medline]
Hagerman, P. J. 1990. Sequence-directed curvature of DNA. Annu. Rev. Biochem. 59:755781.[CrossRef][Medline]
Hagerman, P. J., and B. H. Zimm. 1981. Monte Carlo approach to the analysis of the rotational diffusion of wormlike chains. Biopolymers. 20:14811502.[CrossRef]
Hammermann, M., N. Brun, K. V. Klenin, R. May, K. Toth, and J. Langowski. 1998. Salt-dependent DNA superhelix diameter studied by small angle neutron scattering measurements and Monte Carlo simulations. Biophys. J. 75:30573063.
Hammermann, M., C. Stainmaier, H. Merlitz, U. Kapp, W. Waldeck, G. Chirico, and J. Langowski. 1997. Salt effects on the structure and internal dynamics of superhelical DNAs studied by light scattering and Brownian dynamic. Biophys. J. 73:26742687.
Heddle, J. G., S. Mitelheiser, A. Maxwell, and N. H. Thomson. 2004. Nucleotide binding to DNA gyrase causes loss of DNA wrap. J. Mol. Biol. 337:597610.[CrossRef][Medline]
Horowitz, D. S., and J. C. Wang. 1984. Torsional rigidity of DNA and length dependence of the free energy of DNA supercoiling. J. Mol. Biol. 173:7591.[CrossRef][Medline]
Kampranis, S. C., A. D. Bates, and A. Maxwell. 1999. A model for the mechanism of strand passage by DNA gyrase. Proc. Natl. Acad. Sci. USA. 96:84148419.
Kirkegaard, K., and J. C. Wang. 1981. Mapping the topography of DNA wrapped around gyrase by nucleolytic and chemical probing of complexes of unique DNA sequences. Cell. 23:721729.[CrossRef][Medline]
Klenin, K., M. Hammermann, and J. Langowski. 2000. Modeling dynamic light scattering of supercoiled DNA. Macromolecules. 33:14591466.[CrossRef]
Klenin, K., J. Langowski, and A. Vologodskii. 2002. Computational analysis of the chiral action of type II DNA topoisomerases. J. Mol. Biol. 320:359367.[CrossRef][Medline]
Klenin, K. V., A. V. Vologodskii, V. V. Anshelevich, A. M. Dykhne, and M. D. Frank-Kamenetskii. 1988. Effect of excluded volume on topological properties of circular DNA. J. Biomol. Struct. Dyn. 5:11731185.[Medline]
Klenin, K. V., A. V. Vologodskii, V. V. Anshelevich, A. M. Dykhne, and M. D. Frank-Kamenetskii. 1991. Computer simulation of DNA supercoiling. J. Mol. Biol. 217:413419.[CrossRef][Medline]
Klenin, K. V., A. V. Vologodskii, V. V. Anshelevich, V. Y. Klisko, A. M. Dykhne, and M. D. Frank-Kamenetskii. 1989. Variance of writhe for wormlike DNA rings with excluded volume. J. Biomol. Struct. Dyn. 6:707714.[Medline]
Kreuzer, K. N., and N. R. Cozzarelli. 1980. Formation and resolution of DNA catenanes by DNA gyrase. Cell. 20:245254.[CrossRef][Medline]
Levene, S. D., and D. M. Crothers. 1986. Topological distributions and the torsional rigidity of DNA. A Monte Carlo study of DNA circles. J. Mol. Biol. 189:7383.[CrossRef][Medline]
Levine, C., H. Hiasa, and K. J. Marians. 1998. DNA gyrase and topoisomerase. IV. Biochemical activities, physiological roles during chromosome replication, and drug sensitivities. Biochim. Biophys. Acta. 1400:2943.[Medline]
Liu, L. F., and J. C. Wang. 1978a. DNA-DNA gyrase complex: the wrapping of the DNA duplex outside the enzyme. Cell. 15:979984.[CrossRef][Medline]
Liu, L. F., and J. C. Wang. 1978b. Micrococcus luteus DNA gyrase: active components and a model for its supercoiling of DNA. Proc. Natl. Acad. Sci. USA. 75:20982102.
Marians, K. J. 1987. DNA gyrase-catalyzed decatenation of multiply linked DNA dimers. J. Biol. Chem. 262:1036210368.
Mizuuchi, K., L. M. Fisher, M. H. O'Dea, and M. Gellert. 1980. DNA gyrase action involves the introduction of transient double-strand breaks into DNA. Proc. Natl. Acad. Sci. USA. 77:18471851.
Morrison, A., and N. R. Cozzarelli. 1981. Contacts between DNA gyrase and its binding site on DNA: features of symmetry and asymmetry revealed by protection from nucleases. Proc. Natl. Acad. Sci. USA. 78:14161420.
Murasugi, K. 1996. Knot Theory and Its Applications. Birkhauser, Boston, MA.
Orphanides, G., and A. Maxwell. 1994. Evidence for a conformational change in the DNA gyrase-DNA complex from hydroxyl radical footprinting. Nucleic Acids Res. 22:15671575.
Peng, H., and K. J. Marians. 1995. The interaction of Escherichia coli topoisomerase IV with DNA. J. Biol. Chem. 270:2528625290.
Reece, R. J., and A. Maxwell. 1991. DNA gyrase: structure and function. Crit. Rev. Biochem. Mol. Biol. 26:335375.[Medline]
Rybenkov, V. V., N. R. Cozzarelli, and A. V. Vologodskii. 1993. Probability of DNA knotting and the effective diameter of the DNA double helix. Proc. Natl. Acad. Sci. USA. 90:53075311.
Rybenkov, V. V., C. Ullsperger, A. V. Vologodskii, and N. R. Cozzarelli. 1997a. Simplification of DNA topology below equilibrium values by type II topoisomerases. Science. 277:690693.
Rybenkov, V. V., A. V. Vologodskii, and N. R. Cozzarelli. 1997b. The effect of ionic conditions on DNA helical repeat, effective diameter, and free energy of supercoiling. Nucleic Acids Res. 25:14121418.
Rybenkov, V. V., A. V. Vologodskii, and N. R. Cozzarelli. 1997c. The effect of ionic conditions on the conformations of supercoiled DNA. II. Equilibrium catenation. J. Mol. Biol. 267:312323.[CrossRef][Medline]
Rybenkov, V. V., A. V. Vologoskii, and N. R. Cozzarelli. 1997d. The effect of ionic conditions on the conformations of supercoiled DNA. I. Sedimentation analysis. J. Mol. Biol. 267:299311.[CrossRef][Medline]
Shaw, S. Y., and J. C. Wang. 1993. Knotting of a DNA chain during ring closure. Science. 260:533536.
Stigter, D. 1977. Interactions of highly charged colloidal cylinders with applications to double-stranded DNA. Biopolymers. 16:14351448.[CrossRef][Medline]
Taylor, W. H., and P. J. Hagerman. 1990. Application of the method of phage T4 DNA ligase-catalyzed ring-closure to the study of DNA structure. II. NaCl-dependence of DNA flexibility and helical repeat. J. Mol. Biol. 212:363376.[CrossRef][Medline]
Ullsperger, C., and N. R. Cozzarelli. 1996. Contrasting enzymatic activities of topoisomerase IV and DNA gyrase from Escherichia coli. J. Biol. Chem. 271:3154931555.
Vologodskii, A. V. 1994. DNA extension under the action of an external force. Macromolecules. 27:56235625.[CrossRef]
Vologodskii, A. V., and N. R. Cozzarelli. 1993. Monte Carlo analysis of the conformation of DNA catenanes. J. Mol. Biol. 232:11301140.[CrossRef][Medline]
Vologodskii, A. V., and N. R. Cozzarelli. 1995. Modeling of long-range electrostatic interactions in DNA. Biopolymers. 35:289296.[CrossRef][Medline]
Vologodskii, A. V., S. D. Levene, K. V. Klenin, M. D. Frank-Kamenetskii, and N. R. Cozzarelli. 1992. Conformational and thermodynamic properties of supercoiled DNA. J. Mol. Biol. 227:12241243.[CrossRef][Medline]
Vologodskii, A. V., A. V. Lukashin, and M. D. Frank-Kamenetskii. 1975. Topological interaction between polymer chains. Sov. Phys. JETP. 40:932936.
Vologodskii, A. V., A. V. Lukashin, M. D. Frank-Kamenetskii, and V. V. Anshelevich. 1974. Problem of knots in statistical mechanics of polymer chains. Sov. Phys. JETP. 39:10591063.
Vologodskii, A. V., and J. F. Marko. 1997. Extension of torsionally stressed DNA by external force. Biophys. J. 73:123132.
Vologodskaia, M., and A. Vologodskii. 2002. Contribution of the intrinsic curvature to measured DNA persistence length. J. Mol. Biol. 317:205213.[CrossRef][Medline]
Vologodskii, A. V., W. Zhang, V. V. Rybenkov, A. A. Podtelezhnikov, D. Subramanian, J. D. Griffith, and N. R. Cozzarelli. 2001. Mechanism of topology simplification by type II DNA topoisomerases. Proc. Natl. Acad. Sci. USA. 98:30453049.
Wang, J. C. 1996. DNA topoisomerases. Annu. Rev. Biochem. 65:635695.[CrossRef][Medline]
Wang, J. C. 1998. Moving one DNA double helix through another by a type II DNA topoisomerase: the story of a simple molecular machine. Quart. Rev. Biophys. 31:107144.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |