help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Biophys J. BioFAST on August 31, 2004.
doi:10.1529/biophysj.104.042994
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.104.042994v1
87/5/3066    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Vologodskii, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Vologodskii, A.
Biophysical Journal 87:3066-3073 (2004)
© 2004 The Biophysical Society

Computational Analysis of DNA Gyrase Action

Alexander Vologodskii

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
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
DNA gyrase introduces negative supercoiling into circular DNA by catalyzing the passage of one DNA segment through another. The efficiency of the reaction is many times higher than that of other topological transformations. We analyze, by a computer simulation, the reaction selectivity for a model of DNA gyrase action that assumes existence of a free loop between the G- and T- DNA segments participating in the reaction. A popular model of this type assumed that the selectivity can be provided by the conformation of the DNA segment wrapped around the enzyme into the right-handed helix turn (G-segment). We simulated the distribution of the reaction products for this model. Equilibrium sets of DNA conformations with one segment of the double helix wrapped around the enzyme were constructed. From these sets we selected conformations that had a second segment properly juxtaposed with the first one. Assuming that the juxtapositions result in the strand-passing reaction, we calculated the reaction products for all these conformations. The results show that different products have to be formed if the enzyme acts according to the model. This conclusion can be extended for any model with a free loop between the G- and T-segments. An alternative model that is consistent with the major experimental observations and the computational analysis, is suggested.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
DNA gyrase is a bacterial enzyme that introduces negative supercoiling into circular DNA (Gellert et al., 1976Go). It belongs to the type II DNA topoisomerases that catalyze passing one double-stranded DNA segment through another one. The free energy required for supercoiling comes from ATP hydrolysis coupled with the strand-passing reaction. The enzyme is important for maintaining a certain level of DNA supercoiling inside bacterial cells and for DNA replication (Levine et al., 1998Go; Reece and Maxwell, 1991Go; Wang, 1996Go).

The active form of the Escherichia coli enzyme consists of two A and two B subunits (reviewed by Reece and Maxwell, 1991Go). 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., 1981Go; Kirkegaard and Wang, 1981Go; Liu and Wang, 1978aGo). 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, 1991Go; Wang, 1998Go). This model of the G-segment conformation is supported by the results from nuclease protection experiments (Liu and Wang, 1978aGo; Morrison and Cozzarelli, 1981Go) as well as by DNase I (Fisher et al., 1981Go; Kirkegaard and Wang, 1981Go) and hydroxyl radical footprinting (Orphanides and Maxwell, 1994Go). 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., 1999Go; Liu and Wang, 1978bGo; Peng and Marians, 1995Go). Quantitative analysis of these experiments showed that each complex between the enzyme and DNA increases the linking number in the ligated molecules by 0.5–0.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, 1979Go). If the T-segment belongs to second circular DNA, the reaction can link or unlink the molecules.



View larger version (31K):
[in this window]
[in a new window]
 
FIGURE 1  The model of DNA gyrase action. The enzyme wraps a DNA segment, G-segment (red), around itself. Thermal motion brings another segment, T-segment, into the entrance gate of the enzyme and the segment passes through the temporary cleaved G-segment, then the break is resealed (Wang, 1998Go). Depending on DNA conformation that precedes the T-segment entering, the reaction can result in creating two negative supercoils (from conformation in panel A), two positive supercoils (B), or creating positive trefoil knot (C) and even more complex products.

 
There are data, however, that do not agree with this model. The strongest result of this kind, obtained by Bates and Maxwell (1989)Go, is the ability of the enzyme to introduce supercoiling in the very small DNA circles, 174 bp in length. Clearly, in this case, the strand passing cannot follow the model simply because the contour length of the circles is insufficient to form the starting conformation (Fig. 1 A). Indeed, if ~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., 2004Go). For us, however, this is an example of a model in which the G- and T-segments are separated by a free DNA loop. The major conclusion of the work will be essentially the same for any model of this type. In the alternative type of models the G- and T-segments belong to one DNA stretch whose conformation is completely specified by the interaction with DNA gyrase.

Although DNA gyrase can catenate, decatenate, and unknot circular DNA molecules (Kreuzer and Cozzarelli, 1980Go; Mizuuchi et al., 1980Go), the efficiency of these reactions is very low (Marians, 1987Go; Ullsperger and Cozzarelli, 1996Go). 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, 1981Go; Hagerman and Zimm, 1981Go; Rybenkov et al., 1997dGo), DNA cyclization (Hagerman, 1990Go; Levene and Crothers, 1986Go; Taylor and Hagerman, 1990Go; Vologodskaia and Vologodskii, 2002Go), equilibrium distributions of topological states (Klenin et al., 1988Go, 1989Go; Rybenkov et al., 1993Go, 1997cGo; Shaw and Wang, 1993Go; Vologodskii and Cozzarelli, 1993Go), elasticity of single molecules (Vologodskii, 1994Go; Vologodskii and Marko, 1997Go), and light and neutron scattering data on supercoiled DNA (Gebe et al., 1996Go; Hammermann et al., 1998Go, 1997Go; Klenin et al., 2000Go). 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., 2002Go; Klenin et al., 2002Go; Vologodskii et al., 2001Go). 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)Go and recently by Corbett et al. (2004)Go. We discuss and elaborate important details of such a model.


    METHODS OF COMPUTATION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
DNA model and simulation procedure
Circular DNA was modeled as a discrete worm-like chain consisting of rigid cylinders of equal length with harmonic potentials with respect to the angles between adjacent segments and the value of torsional deformation, {Delta}Tw (Klenin et al., 1991Go; Vologodskii et al., 1992Go). Electrostatic interactions between the segments were taken into account through a hard core potential specified by DNA effective diameter (Vologodskii and Cozzarelli, 1995Go; Vologodskii et al., 1992Go). The equilibrium ensemble of chain conformations was simulated by the Metropolis procedure (Vologodskii et al., 1992Go). Each simulation was performed for the chains with a particular topology and {Delta}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., 1974Go). If the value of 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.



View larger version (18K):
[in this window]
[in a new window]
 
FIGURE 2  Analysis of the simulated conformations. (A) Juxtaposition of the G-segment with a potential T-segment. It is assumed that the G-segment is bounded with the enzyme and thus has a rigid conformation. The G-segment consists of five straight segments (red) that form approximately one turn of the right-handed helix. A segment (yellow) is treated as a potential T-segment if the distance {rho} and the angle {phi} are less than and correspondingly. (B) Topological test of the juxtaposition. Two closed contours were formed by closing the ends of the G-segment and by closing the remaining part of the model chain. The closings of both parts of the chain were done by small deformations of the end segments. It was then tested that the contours are not linked. (C) Extension of the G-segment during the analysis of simulated conformations. The conformation of the G-segment was extended by doubling the length of the second and fourth straight segments of the loop and keeping the length and direction of the third straight segment. The extended and original conformations are shown by bright and dim red, correspondingly.

 
The subsets of conformations, having a segment juxtaposed with the G-segment, were selected from the constructed equilibrium sets for the analysis of potential outcome of the strand-passing reaction.

Parameters
All computations were performed for a DNA persistence length of 50 nm (Hagerman, 1988Go), a torsional rigidity, C, of 3.0 x 10–19 erg·cm (Horowitz and Wang, 1984Go; Klenin et al., 1989Go), 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., 1981Go; Rybenkov et al., 1993Go, 1997bGo; Stigter, 1977Go). 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:

  1. The distance between the middle of the testing segment and the middle of the G-segment contour was less than a distance (Fig. 2 A);
  2. The angle {phi} (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);
  3. There is no linkage between small closed contour formed by closing the ends of the G-segment and the contour formed by closing the rest of the chain without the G-segment (Fig. 2 B). To test this we calculated the Alexander polynomial for these two closed contours, for (Vologodskii et al., 1975Go). This test excluded conformations when the tested segment was inside the loop formed by the G-segment.
  4. To determine the product that could be obtained by the strand-passing reaction from a particular conformation of the chain, the G-segment was extended as shown in Fig. 2 C. We tested that two circular contours, created by closing the extended G-segment and by closing the chain with the excluded G-segment, form the simplest topological link. This was done by calculating 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, {delta}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 {delta}Wr (Kleninet al., 2002Go). 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 {delta}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 {delta}Wr (see, e.g., Murasugi, 1996Go; in particular Example 6.4.2 and Theorem 6.4.7).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
We assume that the reaction catalyzed by gyrase is not diffusion limited, and thus can be analyzed in terms of the equilibrium distribution of DNA conformations. This assumption seems valid for the majority of enzymatic reactions. An initial state for the simulation corresponds to the complex between gyrase and DNA in which a segment of circular DNA is wrapped around DNA gyrase. We assume that the complex is ready to absorb another DNA segment, the T-segment. We also assume that the gate for the T-segment is located in the middle of the wrapped G-segment. The conformation of the G-segment is kept unchanged during the simulation. Conformations of the rest of the circular DNA are sampled to simulate the equilibrium distribution. We select from the constructed set of conformations those that have a segment properly juxtaposed with the assumed gate located in the middle of the G-segment (Fig. 3). Each of these conformations has equal probability to become a subject of the strand-passing reaction. This is true because the distribution of DNA conformations is a continuous function and, thus, the probabilities to find a segment in any two very close positions are nearly equal. We imitate the strand-passing reaction for each of the selected conformations and calculate the topological state of the product (see "Methods of Computation"). The topological state of circular DNA is specified by two parameters, the linking number difference of the complementary strands, {Delta}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.



View larger version (25K):
[in this window]
[in a new window]
 
FIGURE 3  Typical simulated conformation of DNA 3.5 kb in length with a segment forming a complex with DNA gyrase. The G-segment, shown by red, is assumed to be wrapped around the enzyme (not shown). The conformation of the segment was kept unchanged during the simulation whereas conformations of the rest of the model chain were sampled according to the equilibrium distribution. There is another segment in the shown conformation, which is juxtaposed with the break (not shown) in the middle of the G-segment. If the juxtaposed segment were a T-segment, the strand-passing reaction for the shown conformation would result in the introduction of two negative supercoils.

 
First, we calculate the distributions for relaxed ({Delta}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.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 4  Computed distributions of the reaction products for the gyrase model. The data correspond to one strand-passing reaction in each DNA molecule. (A) The product distributions for initially relaxed unknotted molecules of different lengths. (B) The product distributions for the molecules whose initial state corresponds to positive trefoil with {Delta}Lk = 2. This state is the second major product obtained from relaxed unknotted molecules.

 
One might suggest that the positive trefoils formed after the first strand-passing reaction will be unknotted in the second round of the process. To test this possibility we simulated the product distribution formed from positive trefoils (Fig. 4 B). The results of this simulation show that the second strand-passing reaction would rarely return the model chains to an unknotted state. In most cases the second reaction results in formation of more complex knots.

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.



View larger version (11K):
[in this window]
[in a new window]
 
FIGURE 5  Computed distributions of the reaction products for the gyrase model. The distributions were obtained for initially unknotted molecules 3.5 kb in length with {Delta}Lk of –2 (A), –4 (B), and –6 (C). In each panel the bars correspond to adding two (–) supercoils and keeping the unknotted state (bar 1), adding two (+) supercoils and keeping the unknotted state (bar 2), forming negative trefoil and decreasing Lk by two (bar 3), and forming positive trefoil and increasing Lk by two (bar 4). The product distributions correspond to one strand-passing reaction in each molecule.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The simulation results clearly show that the model of DNA gyrase action shown in Fig. 1 cannot explain the high selectivity of the strand-passing reaction catalyzed by the enzyme. If DNA gyrase were to act according to this model, it would produce a large fraction of knotted molecules. Because negatively supercoiled molecules overwhelmingly dominate the population of DNA gyrase products, we have to conclude that the model presented in Fig. 1 is incorrect. The failure of the model to provide a unique reaction outcome results from the assumption that there are free loops between the G- and T-segments that do not interact with protein. Conformations of these loops have to depend on the properties of the double helix only. Because DNA is a sufficiently flexible molecule on the scale of thousands of basepairs, the loops accept many different conformations that correspond to different products of the strand-passing reaction. Negative supercoiling can represent the overwhelming reaction outcome only if the G- and T-segments belong to one DNA stretch whose conformation is completely specified by the interaction with DNA gyrase. Models of this kind were suggested by Kampranis et al. (1999)Go and recently by Corbett et al. (2004)Go. Below we elaborate on a version of such a model with some additional features.

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, 1991Go). 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, 1989Go). 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)Go 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.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 6  Introduction of two negative supercoils into a DNA minicircle by DNA gyrase. The model assumes that DNA-protein interaction causes the motion of the upper part of the DNA circle (yellow) through the break introduced in the G-segment (red). The diagram addresses only conformational changes in DNA during the strand-passing reaction.

 


View larger version (14K):
[in this window]
[in a new window]
 
FIGURE 7  The model of DNA gyrase action. A continuous stretch of the double helix interacting with the protein includes both the G- and T-segments, and the DNA-protein interaction specifies the conformational changes in the stretch. Introduction of negative supercoiling is the only possible outcome of the reaction in this model. To keep the symmetry of the complex we assume that the required DNA-protein interaction can occur at either right (A) or left side (B) of the enzyme. The diagram does not address conformational changes in DNA gyrase.

 
The model shown in Fig. 7 A does not fit the requirement of the symmetry in the complex, however. More than this, half of the G-segment does not play an active role in the reaction, and it is not clear why its extensive interaction with the protein is needed. To overcome this problem we suggest that either of two halves of the enzyme can interact with the stretch of G- and T-segments (see Fig. 7). This would make both sides of the bound segment equivalent for footprinting or nuclease digestion. Of course, in any particular DNA stretch bound with the enzyme, the DNA sequence can determine preferential formation of one of two possible loops. Therefore the level of the symmetry observed in footprinting or nuclease digestion experiments should depend on the G-segment sequence, and this could explain certain discrepancy between different nuclease protection and footprinting experiments (Fisher et al., 1981Go; Heddle et al., 2004Go; Kirkegaard and Wang, 1981Go; Liu and Wang, 1978aGo; Morrison and Cozzarelli, 1981Go; Orphanides and Maxwell, 1994Go).

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., 1993Go). 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., 1997aGo). 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., 2001Go). 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., 2001Go). 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
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
I thank Nicholas Cozzarelli, Neville Kallenbach, and Yuri Yarim-Agaev for helpful discussions and comments.

The work was supported by National Institutes of Health grant GM 54215.

Submitted on March 17, 2004; accepted for publication August 26, 2004.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS OF COMPUTATION
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Bates, A. D., and A. Maxwell. 1989. DNA gyrase can supercoil DNA circles as small as 174 base-pairs. EMBO J. 8:1861–1866.[Medline]

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:1305–1328.[CrossRef][Medline]

Brown, P. O., and N. R. Cozzarelli. 1979. A sign inversion mechanism for enzymatic supercoiling of DNA. Science. 206:1081–1083.[Abstract/Free Full Text]

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:7293–7298.[Abstract/Free Full Text]

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:4165–4169.[Abstract/Free Full Text]

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:105–128.[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:3872–3876.[Abstract/Free Full Text]

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:515–527.[CrossRef][Medline]

Hagerman, P. J. 1981. Investigation of the flexibility of DNA using transient electric birefringence. Biopolymers. 20:1503–1535.[CrossRef][Medline]

Hagerman, P. J. 1988. Flexibility of DNA. Annu. Rev. Biophys. Biophys. Chem. 17:265–286.[CrossRef][Medline]

Hagerman, P. J. 1990. Sequence-directed curvature of DNA. Annu. Rev. Biochem. 59:755–781.[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:1481–1502.[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:3057–3063.[Abstract/Free Full Text]

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:2674–2687.[Abstract/Free Full Text]

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:597–610.[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:75–91.[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:8414–8419.[Abstract/Free Full Text]

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:721–729.[CrossRef][Medline]

Klenin, K., M. Hammermann, and J. Langowski. 2000. Modeling dynamic light scattering of supercoiled DNA. Macromolecules. 33:1459–1466.[CrossRef]

Klenin, K., J. Langowski, and A. Vologodskii. 2002. Computational analysis of the chiral action of type II DNA topoisomerases. J. Mol. Biol. 320:359–367.[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:1173–1185.[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:413–419.[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:707–714.[Medline]

Kreuzer, K. N., and N. R. Cozzarelli. 1980. Formation and resolution of DNA catenanes by DNA gyrase. Cell. 20:245–254.[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:73–83.[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:29–43.[Medline]

Liu, L. F., and J. C. Wang. 1978a. DNA-DNA gyrase complex: the wrapping of the DNA duplex outside the enzyme. Cell. 15:979–984.[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:2098–2102.[Abstract/Free Full Text]

Marians, K. J. 1987. DNA gyrase-catalyzed decatenation of multiply linked DNA dimers. J. Biol. Chem. 262:10362–10368.[Abstract/Free Full Text]

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:1847–1851.[Abstract/Free Full Text]

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:1416–1420.[Abstract/Free Full Text]

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:1567–1575.[Abstract/Free Full Text]

Peng, H., and K. J. Marians. 1995. The interaction of Escherichia coli topoisomerase IV with DNA. J. Biol. Chem. 270:25286–25290.[Abstract/Free Full Text]

Reece, R. J., and A. Maxwell. 1991. DNA gyrase: structure and function. Crit. Rev. Biochem. Mol. Biol. 26:335–375.[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:5307–5311.[Abstract/Free Full Text]

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:690–693.[Abstract/Free Full Text]

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:1412–1418.[Abstract/Free Full Text]

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:312–323.[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:299–311.[CrossRef][Medline]

Shaw, S. Y., and J. C. Wang. 1993. Knotting of a DNA chain during ring closure. Science. 260:533–536.[Abstract/Free Full Text]

Stigter, D. 1977. Interactions of highly charged colloidal cylinders with applications to double-stranded DNA. Biopolymers. 16:1435–1448.[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:363–376.[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:31549–31555.[Abstract/Free Full Text]

Vologodskii, A. V. 1994. DNA extension under the action of an external force. Macromolecules. 27:5623–5625.[CrossRef]

Vologodskii, A. V., and N. R. Cozzarelli. 1993. Monte Carlo analysis of the conformation of DNA catenanes. J. Mol. Biol. 232:1130–1140.[CrossRef][Medline]

Vologodskii, A. V., and N. R. Cozzarelli. 1995. Modeling of long-range electrostatic interactions in DNA. Biopolymers. 35:289–296.[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:1224–1243.[CrossRef][Medline]

Vologodskii, A. V., A. V. Lukashin, and M. D. Frank-Kamenetskii. 1975. Topological interaction between polymer chains. Sov. Phys. JETP. 40:932–936.

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:1059–1063.

Vologodskii, A. V., and J. F. Marko. 1997. Extension of torsionally stressed DNA by external force. Biophys. J. 73:123–132.[Abstract/Free Full Text]

Vologodskaia, M., and A. Vologodskii. 2002. Contribution of the intrinsic curvature to measured DNA persistence length. J. Mol. Biol. 317:205–213.[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:3045–3049.[Abstract/Free Full Text]

Wang, J. C. 1996. DNA topoisomerases. Annu. Rev. Biochem. 65:635–695.[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:107–144.[CrossRef][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.104.042994v1
87/5/3066    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Vologodskii, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Vologodskii, A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2004 by the Biophysical Society.